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Foreword 


The following report is an assessment of scientific predictions for the atmospheric effects of a proposed 
fleet of high-speed civil transport (HSCT) aircraft. It was commissioned by the Atmospheric Effects of 
Stratospheric Aircraft (AESA) element of the National Aeronautics and Space Administration (NASA) 
High-Speed Research Program (HSRP), which has technology goals that, if achieved, would promise 
environmental compatibility and economic practicability for the HSCT fleet. When the program was 
conceived in 1989, this latest revival of interest in supersonic transport (SST) was greeted with healthy 
skepticism by those who well recalled the related atmospheric studies of the 1970s. For the technology 
of that era, it was predicted that nitrogen oxides (NO x ) from a fleet of SSTs would do significant 
damage to stratospheric ozone [Johnston, 1971]. 

Our understanding of the atmosphere had grown significantly by the beginning of the AESA program 
[Johnston etal, 1991; Douglass etal., 1991], and today we are on the verge of having predictive tools 
that can be applied with a high degree of confidence. However, today's best informed opinion 
regarding HSCT technology is the same as before; that the SSTs studied in the 1970s would have been 
an environmental problem. But in reading the present report and projecting to where the science may 
be in a few years, it is important to realize that aircraft technology has also progressed and will continue 
to advance. 

Earlier atmospheric predictions assumed that jet engine combustors produced NO x at a rate of about 30 
grams of equivalent nitrogen dioxide (NO 2 ) per kilogram of fuel burned. The original HSRP target 
emission, 5 grams per kilogram of fuel burned [Shaw et al., 1993, Sweetman, 1995, Williams, 1993], 
has now been achieved in tests of practical combustor configurations. At that low level of emissions, 
current atmospheric models suggest that HSCT technology indeed can be environmentally compatible 
[Stolarski and Wesoky, 1993a]. Therefore, it is the parallel development of science and technology that 
offers the promise of success for the HSRP. 

Recognition of the progress made in understanding the atmospheric processes to which an HSCT fleet 
would contribute, and the technological means for controlling resulting emissions, is demonstrated by a 
September 1994 recommendation from the International Coordinating Council of Aerospace Industry 
Associations to the Emissions Working Group (WG3) of the International Civil Aviation Organization's 
(ICAO) Committee on Aviation Environmental Protection (CAEP): 

WG3 should recommend to CAEP that it charter a sub-group, with appropriate Terms of 
Reference, to develop Annex 16 emissions certification requirements for future 
supersonic aeroplanes, consistent with technical practicability and economic 
reasonability, for consideration at CAEP/4. 

The next meeting of CAEP (i.e., CAEP/3) will be in December 1995, with CAEP/4 likely to occur three 
to five years later. 

Although this progress in the regulatory forum is also very promising for HSCT technology, it must be 
recognized that ICAO has indicated [CAEP, 1992] that any future emission standard will likely occur 
only “when the environmental need has been accepted by an international consensus (e.g., by 
UNEP/WMO).” And, although significant “developments [have] changed decidedly the scientific 
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community's judgment with respect to the expected impact of the NO x component of the HSCT 
effluent” [Anderson, 1995], such a consensus does not yet exist. The most current ozone assessment of 
the United Nations Environment Programme (UNEP) and World Meteorological Organization (WMO) 
[WMO, 1995] clearly recognizes the progress of the NASA program and other research, but also states 
that there are “important uncertainties in supersonic assessments.” 

These uncertainties have been highlighted in an evaluation by the National Research Council (NRC) of 
the National Academy of Sciences that NASA requested in 1992. A panel of ten experts, including four 
from Europe, received extensive briefings from NASA managers and investigators before an interim 
assessment report [NASA, 1993] was transmitted to them in June 1993. In Congressional testimony 
[Graedel, 1994] regarding the NRC [1994] report, the panel chairman noted that the then 

current assessment of environmental effects resulting from HSCT fleet emissions reveals 
no paramount concerns, but this assessment is clearly preliminary ... Although many 
AESA studies during the next year or two will increase our understanding of aspects of 
the atmosphere-aircraft interaction, the assessment's overall reliability will not change 
significantly by 1995. 

The present report seeks to incorporate AESA results from the last two years into an updated 
assessment. It will also serve to initiate the second phase of AESA, in direct response to the NRC 
evaluation [Graedel, 1 994] which recommended that 

The extension of AESA through 1998, coupled with strong program emphasis on 
resolving key uncertainties through a combination of airborne missions and computer 
model incorporation of mission results, would markedly increase the confidence one 
could bring to the assessment results. 

Convergence of a confident scientific assessment, the related HSCT technology, and the regulatory 
policy at the end of this decade would be a truly exciting achievement. 


Howard L. Wesoky 

Manager, Atmospheric Effects of Aviation 
Office of Aeronautics 

National Aeronautics and Space Administration 



Executive Summary 


1. What are the concerns about high-speed civil transport operations in the stratosphere that 
motivate this assessment? 

The operation of high-speed civil transport (HSCT) aircraft will introduce exhaust gases and 
particulates directly into the stratosphere. These exhaust components could alter stratospheric ozone 
and the Earth’s climate. 

• Important emissions from HSCTs include water vapor, carbon dioxide (CO2), nitrogen oxides 
(NO x ), sulfur oxides (SO x ), carbon monoxide (CO), nonmethane hydrocarbons (NMHCs), and soot. 
The atmospheric impact of these emissions depends on how they change the steady-state abundance 
over background atmospheric amounts and their role in atmospheric photochemical, dynamical, and 
radiative processes. Atmospheric ozone change is possible because these HSCT emissions directly 
affect photochemical processes which control the ozone abundance. Climate change is also possible 
because these HSCT emissions are directly involved in the radiative balance of the atmosphere. 

2. What processes are important in today's atmosphere and in a future atmosphere perturbed 
by HSCT emissions? 

The amount of stratospheric ozone is determined by photochemical production and loss processes, and 
by the transport of air throughout the atmosphere. Photochemical loss of ozone is dominated by 
catalytic reactions involving NO x , hydrogen oxides (HO x ), and halogen radicals (chlorine oxides 
(C 10 x ) and bromine oxides (BrO x )). Heterogeneous reactions which occur on or in stratospheric 
aerosol particles play an important role by reducing the ozone loss due to NO x and increasing that due 
to HO x and halogens. Atmospheric circulation determines the time spent by air in regions of 
photochemical loss and the distribution of exhaust gases emitted by a fleet of HSCTs. 

• NO x catalytically destroys stratospheric ozone but also inhibits ozone destruction by the HO x and 
halogen radical catalytic cycles through the formation of more stable gases. The effect of adding 
NO x to the atmosphere by HSCTs is sensitive to the balance of these two effects. In the 
stratosphere above ~26 km, NO x catalysis is the most important ozone loss process; increased NO x 
leads to increased local ozone loss. In the lower stratosphere, NO x interference with the dominant 
HO x and halogen catalysis of ozone destruction is more important. For increased NO x (e.g., from 
HSCT exhaust), the sign and magnitude of the local chemical ozone change in the lower 
stratosphere depend on the background amounts of trace gases and aerosol surface area. 

• Reactions on sulfate aerosols composed of sulfuric acid (H2SO4) and water (H2O) play an 
important role in determining the abundance of NO x . Because of these reactions, the concentration 
of NO x depends on the variable input of volcanic SO2 to the stratosphere. Observations show that 
the concentration of NO x falls in response to increases in aerosol surface area. Because increases in 
sulfur, water vapor, and soot from HSCT emissions will change the surface area, composition, and 
number of particles in the lower stratosphere, the emissions indirectly influence ozone loss rates. In 
low-temperature regions near the poles, both liquid and frozen aerosols form that contain nitric acid 
(HNO3) in addition to H2O and H2SO4. Reactions on these aerosols enhance halogen-catalyzed 
loss of ozone and may sediment from the stratosphere, effectively reducing NO x . The increased 
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water vapor and H2SO4 from HSCT exhaust could increase the importance of such processes, 
especially if the stratosphere is cooled by climate change. 

• Transport of air between mid-latitudes and the tropics is restricted, especially between 20 and 28 
km. Transport of air from the primary flight corridors in mid-latitudes into the tropics results in a 
more rapid spread of HSCT exhaust from the lower to upper stratosphere, where the catalytic 
destruction of ozone by NO x is more efficient. 

• For a projected fleet of HSCT aircraft, climate can be affected by changes in the abundance of water 
vapor, soot, or sulfate in the stratosphere. Climate also can be affected by changes in cirrus cloud 
properties or by the radiative consequences of changes in ozone and its vertical distribution. 

3. What has been learned about the character of engine emissions and the chemical and fluid- 
dynamic processing that occur in the aircraft wake? 

In situ aircraft exhaust sampling experiments and modeling of engine processes and the dilution of ex- 
haust gases into the atmosphere are used to define the important chemical and dynamical features of ex- 
haust plumes in the stratosphere. 

• Measurements of reactive nitrogen (NO y ) species, condensation nuclei (CN), and CO 2 were 
recently made in the exhaust of the Concorde supersonic aircraft and the NASA ER-2 subsonic 
aircraft while in flight. The emission indices (Els) were calculated independently of plume 
dynamics because measurements of CO 2 provide a measure of fuel burned. The measured El of 
NO x , expressed as grams of equivalent nitrogen dioxide (NO 2 ) per kilogram of fuel burned, is in 
excellent agreement with earlier results deduced from altitude chamber measurements of exhaust 
from the Concorde engine as part of the Climatic Impact Assessment Program (CIAP). These 
results suggest that the methodology developed to calculate the El of NO x at cruise from altitude 
chamber data will be appropriate to evaluate new HSCT engines. 

• Observations show that NO x is the most abundant NOy species in both supersonic and subsonic air- 
craft plumes. The presence of only a small percentage of HNO3 and nitrous acid (HONO) is in 
agreement with observations of HO x in the plume and theoretical model calculations. 

• The large number of detectable particles in the Concorde plume conflicts with the results of plume 
models that predict the formation of a large number of much smaller sulfate particles and a smaller 
fractional conversion of SO x to condensed sulfate. These results indicate that current understanding 
of plume chemistry and particle formation processes is not complete. If new HSCT aircraft produce 
particles at a rate comparable to the Concorde, significant increases in particle number and surface 
area could occur in the lower stratosphere of the Northern Hemisphere. Including these increases in 
predictive models could change calculated ozone and climate impacts from an HSCT fleet. 

4. What are the predicted atmospheric changes associated with HSCTs? 

Five two-dimensional (2-D) photochemical models were used to calculate the impact of an HSCT fleet 
for a variety of cases including sensitivity tests for the El of NO x , cruise altitude, background 
atmospheric chlorine amount, and fleet size. Individual models were used to test sensitivity to a variety 
of effects not included in these basic predictive calculations. 
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• The range of model results is shown in Table ES-1. Conditions include a baseline fleet of 500 
Mach 2.4 aircraft burning 8.2 x 10 10 kg of fuel per year in an atmosphere with a background 
stratospheric aerosol amount equivalent to a relatively clean period between major volcanic 
eruptions. The models do not include the formation or reactivity of polar stratospheric cloud (PSC) 
particles. The calculated changes in ozone in Table ES-1 result from a decrease in the middle and 
upper stratosphere and an increase in the lower stratosphere and troposphere. The time required for 
photochemistry and transport to effect these ozone changes precludes large ozone changes near the 
assumed HSCT flight corridors. The EIno x values in Table ES-1 reflect those expected for future 
HSCT combustors (5 g NC> 2 /kg fuel) and those closer to current engines (15 g NO 2 /kg fuel). 

• Laboratory observations continue to refine photolysis and reaction rate coefficients of importance to 
the lower stratosphere and to identify important reactions on surfaces of stratospheric aerosol 
particles. The formulation of the chemistry in the 2-D models is derived from the latest 
recommendation of a NASA panel that conducts a biennial review [DeMore et al . , 1994]. 

• Four of the five models used in the calculations were tested against two benchmarks, one for 
photolysis rates and one for chemistry. These indicate mostly minor differences, implying that the 
range in model results follows, at least in part, from differences in the 2-D representation of 
dynamics. The chemistry benchmark tested the model chemistry integrators against a 
photochemical steady-state model that has been used in detailed comparisons with atmospheric 
measurements of key radicals in volved in the chemistry of stratospheric ozone. 

• When the background aerosol surface area used in the models is increased by a factor of 4 (roughly 
equivalent to the median amount in the stratosphere over the last 2 decades), the calculated 
perturbation to total ozone is less negative by a few tenths of a percent. Variations in future 


Table ES-1. Calculated steady-state total column ozone change (%) in the Northern Hemisphere 
averaged over a year for several HSCT fleet scenarios. 


Scenario 

2-D model results 

Mach number 

Cly 

(ppbv) 

EInOx 

(g NC> 2 /kg fuel) 

Number 

of 

aircraft 

Ozone change (%) a 

2.4 

3 

0 

500 

-0.3 to -0.1 

2.4 

3 

5 

500 

-0.3 to +0.1 

2.4 

3 

5 

1000 

-0.7 to +0.03 

2.4 

3 

10 

500 

-0.5 to 0.0 

2.4 

3 

15 

500 

-1.0 to -0.02 

2.4 

3 

15 

1000 

-2.7 to -0.6 

2.4 Cruise +1 km 

3 

5 

500 

-0.5 to +0.02 

2.4 Cruise -2 km 

3 

5 

500 

-0.06 to +0.1 

2.4 

2 b 

5 

500 

-0.4 to +0.02 


a Range of average values obtained from the five 2-D models used in this assessment. 
b Expected in a 2050 atmosphere. 
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stratospheric aerosol levels, primarily due to the frequency of large volcanic eruptions, are a source 
of uncertainty in HSCT ozone impact predictions. 

• Sulfur emitted by a fleet of HSCTs can increase the sulfate aerosol surface area throughout the 
lower stratosphere. Although the increase in surface area decreases the impact of NO x emissions, it 
also in-creases ozone loss due to the existing chlorine and hydrogen radicals in the stratosphere. If 
the sulfur emitted by HSCTs immediately forms small particles in the plume, then the global 
increase in sulfate surface area is maximized. 

• A model that included aerosol processes was used to calculate ozone changes when 100% of the 
sulfur (EIso 2 = 0-4 g sulfur/kg fuel) from an HSCT fleet was emitted as small particles into a 
volcanically unperturbed stratosphere. The change in total column ozone from the injection of 
small particles depended in a complex way on EIno x and background chlorine amount. For 
EInOx = 5 g/kg fuel, the calculated ozone depletion increased when small particles were formed; the 
overall perturbation became as large or larger than that for the EInox = 15 g/kg fuel case. Inclusion 
of the reaction of BrONC >2 with water on aerosols increased the effect of the small particles. 

• The model results for ozone change in the basic scenarios vary more widely when heterogeneous 
reactions on PSCs are included. This results from uncertainty in PSC microphysical processes and 
from ambiguity in how to best represent such processes in a 2-D model. 

• A general circulation model (GCM) was used to calculate the global surface air temperature and 
stratospheric temperature responses to predicted changes in ozone and water vapor from a fleet of 
HSCTs with EIno x = 15 g NC> 2 /kg fuel. Since the resulting changes in temperature are smaller than 
the climatic fluctuations in the model, they likely would not be detectable in the atmosphere. 

• Climate effects due to sulfur and soot emissions, while likely to be small, have yet to be evaluated. 
Formation of persistent contrails from HSCTs is expected to be infrequent. 

5. How consistent are atmospheric and laboratory measurements with the current representa- 
tion of processes in the 2-D models used to predict the HSCT effects on ozone? 

Atmospheric observations can provide tests of our understanding of the photochemical and transport 
processes responsible for controlling the distribution of ozone and the potential impact of HSCTs. 
Laboratory measurements define the formulation of chemistry in models. Subsets of this formulation 
are tested by comparing the model predictions against a sufficiently complete set of atmospheric 
observations. Although constraints on transport formulations in 2-D models are less direct, 
observations of the distribution of and the correlations between long-lived gases provide a test of these 
formulations. An extensive series of 2-D model tests, conducted in 1992 [Prather and Remsberg, 
1993], led to improvements in the models and an understanding of the differences among them. 
Comparisons to observations made since that time are ongoing, as are further model improvements. 

• Models reproduce the basic features of the ozone column and vertical profile. Consistent with 
observations, the models suggest that ozone concentrations should have decreased in the past 
decade. However, the predicted changes are smaller than those measured by a factor of 1.5 to 2, 
particularly in the lower stratosphere. Inclusion of a realistic representation of the temporal changes 
in stratospheric sulfate aerosol in a model improves the agreement with the observed trends, but the 
predicted changes are still smaller than observed. These differences suggest that important 
uncertainties remain in the ozone changes calculated for an HSCT fleet. 
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• Model simulations reproduce some, but not all, of the basic features of the distribution of long-lived 
trace gases in the stratosphere. Since trace gases are the sources of the free radicals that catalyze 
ozone removal, model-measurement differences reflect inaccuracies in the model transport 
formulation and affect comparisons of reactive species calculations with measurements. 

• Important regions for simulation of atmospheric transport are those with high spatial gradients in 
long-lived species such as N 2 O, NO y , or ozone. Examples are the regions in the stratosphere 
between the tropics and subtropics and at the edge of the polar vortices in winter. Aircraft and 
satellite observations have increased the level of understanding of the gradients and the barriers to 
transport that create them. In general, the 2-D models used to predict HSCT effects exhibit weak 
gradients, suggesting that they have too much mixing across these barriers. 

• Observations of HO x , NO x , C10 x , and BrO x in the stratosphere allow direct evaluation of the photo- 
chemical destruction cycles of ozone. Above 25 km altitude, observations of HO x and NO x free 
radicals agree with photochemical models to within -40% when the models are constrained by 
observed abundances of long-lived chemical species, temperature, aerosol surface area, and 
overhead ozone column. Observations of chlorine monoxide (CIO) at these altitudes are up to a 
factor of 2 lower than model predictions. Below 20 km, midday observations of HO x , NO x , C10 x , 
and BrO x also generally agree with similarly constrained models to within -40%. However, HO x 
observations exceed model calculations by more than a factor of 3 near sunrise and sunset, times 
when HO x concentrations are small. Incorporation of the reaction of bromine nitrate (BrONO^ 
with H 2 O on sulfate aerosols reduces this discrepancy but increases the discrepancy in other 
species. Between about 20 and 25 km, observations of HO x and C10 x radicals are not available to 
provide a rigorous test of models. 

• The observations and modeling of NO x , HO x , and C10 x free radicals in the lower stratosphere agree 
in showing that the concentrations of HO x and C10 x depend on NO x radicals. When the NO x 
concentration is low, such as observed in the lower stratosphere of the Northern Hemisphere during 
the fall of 1992 and spring of 1993, the HO x and C10 x /BrO x reactions dominate the local ozone loss 
rates. Thus, according to model calculations, when NO x from HSCTs is added to these regions, 
reductions in HO x and C10 x /BrO x abundances generally result in reductions in local ozone loss 
rates. 

• Observations show that PSCs occur during the Antarctic winter and spring and Arctic winter. 
Heterogeneous reactions on these PSC surfaces activate large fractions of the available chlorine that 
can then destroy ozone in the presence of sunlight. Widespread ozone destruction does occur in the 
Antarctic because sedimenting PSC particles also denitrify the stratosphere and alter the buffering 
of C10 x by NO x . The proper inclusion of PSC processes in the 2-D predictive models would reduce 
an important uncertainty in the HSCT impact results. 

6. What are the uncertainties in the atmospheric impacts of HSCT operations? 

Several classes of uncertainty are associated with predictions of HSCT impacts made with 2-D models 

as shown in Table ES-1. Uncertainties in some model parameters can be readily quantified and evalu- 
ated directly with the models. Other sources of uncertainty have been identified but not yet quantified. 

Finally, uncertainty from the possibility of missing processes or major errors in understanding cannot 

be quantified, but only discussed in terms of the overall confidence in the assessment. 



• The only quantified uncertainty in this assessment is that due to uncertainties in gas-phase reaction 
rate coefficients. These have been propagated through a model to give an estimated uncertainty of 
~1% (l-o) in the calculated ozone column perturbation in the Northern Hemisphere. Uncertainties 
in heterogeneous chemical reactions and photolysis rates will increase the overall chemical 
uncertainty. Constraining these results with atmospheric observations will reduce the uncertainty. 

• A number of parametric sensitivity tests were performed for processes not included in the basic as- 
sessment scenarios. If the combined range of these tests is taken as an estimate of uncertainty, the 
overall uncertainty due to the processes not included in the basic scenarios is equal to or greater 
than the reaction-rate uncertainty. 

• Uncertainties in microphysics include both the formation of numerous small particles in the HSCT 
exhaust plume and the possibility of increased denitrification in polar regions due to increases in 
condensible gases from the HSCT exhaust. Although the numerous small particles in the HSCT 
plume can indirectly affect ozone loss rates, the particle formation process is poorly understood and 
the uncertainties are large. Although the likelihood of a significant increase in Arctic denitrification 
and in polar ozone loss through chlorine and bromine chemistry is small, the implications of such 
processes for this assessment are large. 

• Uncertainties in model transport remain important and difficult to quantify. Observations of season- 
ally varying and short- and long-lived tracers provide diagnoses of the degree of mixing and rates of 
transport between tropical and mid-latitude regions of the stratosphere, as well as providing 
information on residence times for stratospheric trace species. 

• The impact of ozone and water vapor from HSCTs on climate was calculated to be within the 
climatic fluctuations of the model used. Thus, the expected perturbation and its uncertainty are 
difficult to quantify. 

• The possibility of missing processes in the models creates an unquantifiable uncertainty in the 
HSCT assessment. However, the probability that important processes have been overlooked 
decreases as we continue to explore HSCT perturbations with models and continue to challenge 
these models with atmospheric measurements. 

7. How can uncertainties in the predicted atmospheric effects of HSCTs be reduced? 

Future efforts to reduce the overall assessment uncertainty are expected to focus on the major 
unquantified uncertainties and attempts to quantify them. The primary methodology for improving our 
predictive capability continues to involve testing model representations of the atmosphere through 
comparison with observational data. Substantial progress in that direction will be realized in the near 
term by systematic comparisons with recent data from the Upper Atmospheric Research Satellite 
(UARS), the Atmospheric Laboratory for Applications and Science (ATLAS) on the Space Shuttle, and 
ER-2 aircraft, as well as a host of balloon platforms. Specific additional strategies for addressing the 
major uncertainties identified in this assessment are listed below in no specific order of priority: 

• Extension of the observational database for long-lived tracers and radicals to include a number of 
seasonal cycles and undersampled regions of the atmosphere, notably the summer high latitudes and 
tropics. 



Systematic development of 3-D chemical-transport models for testing transport formulations and 
zonal-symmetry assumptions in 2-D predictive models. 

Inclusion of standardized descriptions of particle formation and evolution in 2-D predictive models 
and comparison with global aerosol distributions derived from existing satellite and aircraft data. 

Further characterization, through laboratory and atmospheric observations, of the composition and 
temporal evolution of aircraft exhaust aerosols in the wake/plume regime, and of the composition 
and microphysics of PSCs. 
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INTRODUCTION 

The High-Speed Civil Transport (HSCT) is a pro- 
posed aircraft that would dramatically shorten inter- 
continental commercial flights. For example, 
transit times from Los Angeles to Tokyo would be 
shortened from over 10 hours to less than 4 1/2 and 
flights from Los Angeles to Sydney from 14 hours 
to 7 1/2 hours, including a refueling stop. Such an 
aircraft might fly at Mach 2.4 (1600 mph) at an alti- 
tude of 20 km (65,000 ft), well into the stratosphere 
where most atmospheric ozone resides. 

The ozone abundance in the stratosphere above us 
is important for three reasons. One is that the 
amount of ultraviolet radiation from the sun 
reaching the surface of the Earth depends upon the 
amount of stratospheric ozone, particularly for the 
wavelengths that cause sunburn and skin cancer. 
The concern has been that reductions in the amount 
of stratospheric ozone might occur as a result of 
emissions from supersonic airliners and industrial 
production of long-lived halogenated compounds, 
thus leading to involuntary exposure of the human 
population to increased risk of skin cancer. The 
second reason is that the enhanced ultraviolet 
radiation at the surface could damage the biosphere, 
particularly the oceanic plankton which are so vital 
to the marine food chains. The third reason is that 
ozone is one of the three atmospheric trace gases 
(water (H2O) and carbon dioxide (CO2) are the 
others) that largely determine the transmission of 
ultraviolet, visible, and infrared radiation through 
the atmosphere, and hence affect climate. 
Chlorofluorocarbons (CFCs), nitrous oxide (N2O), 
and methane (CH4) are other gases which can affect 
climate and whose abundances are being affected 
by human activity. 

Figure 1 is a simple schematic of the atmosphere 
between the ground and 30 km. The sunlight ab- 
sorbed at the surface produces a rapidly 
overturning, convective thermal structure in the 
lower atmosphere, called the troposphere, in which 
temperature decreases with height. However, 
above approximately 12-16 km the temperature 
increases with altitude due to the absorption of 


ultraviolet sunlight by ozone. This region of 
permanent temperature inversion is called the 
stratosphere. Because of this inversion, the 
stratosphere is characterized by slower removal of 
injected pollutants. The boundary between the tro- 
posphere and the stratosphere is marked by a tem- 
perature minimum called the tropopause. Air 
crosses it in an upward direction in the tropics at an 
altitude of -16 km (10 miles), and is dried by the 
very low temperatures there. Although clouds and 
rain can remove soluble molecules from the 
troposphere, they cannot do so in the dry 
stratosphere except for the special case of the very 
low-temperature regions occurring during winter 
over the poles. Air returns to the troposphere at 
mid- to high latitudes, where the height of the 
tropopause is 12 km or less. 

By the mid-1960s, it was realized that the 
distribution of ozone in the stratosphere could not 
be quantitatively understood using only simple 
photochemistry. Balancing the production of ozone 
(O3) from the photolysis of molecular oxygen (O2) 
by sunlight with loss via the recombination of 
atomic oxygen and ozone led to significant overes- 


Ozone Concentration (10 12 molecules' 3 ) 



Figure 1 . A simple schematic of the atmosphere 
between the ground and 30 km. Temperature is 
shown by the solid curve using the bottom scale 
and concentration of ozone molecules is shown as 
a dashed curve using the top scale. Both are for 
mid-latitude conditions from the U.S. Standard 
Atmosphere [1976]. Typical cruise altitudes for 
HSCTs and subsonic aircraft are indicated. 
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timates of the total amount of ozone. In addition, 
this model could not explain quantitatively why 
overhead column abundances of ozone observed at 
high latitudes are larger than in the tropics. 

These calculations showed two things. One was the 
importance of winds in moving ozone from where 
it is produced in the tropics to the polar regions 
where it more easily accumulates. The other was 
that some additional chemical processes must be 
destroying ozone. It is these processes that make it 
possible for industrial chemicals to affect ozone, 
and that together make it a complicated task to 
produce an assessment of the effects of a specific 
activity, such as the operation of a fleet of HSCT 
aircraft. 

Missing from the earlier models was the knowledge 
that certain molecules existing naturally in the 
atmosphere - such as H2O and N2O (produced 
biologically) - could react with energetic oxygen 
atoms in the stratosphere to form reactive molecular 
fragments called free radicals. These radicals, 
although present in very small amounts, can destroy 
many ozone molecules by series of reactions in 
which ozone is consumed but the free radicals are 
not, a process called catalytic destruction. These 
series of reactions can repeat many times as long as 
the free radicals remain in the stratosphere. The 
free radicals of concern were the hydroxyl radical 
(OH) and the hydroperoxy radical (HO2) from 
H2O, and nitric oxide (NO) and nitrogen dioxide 
(N0 2 ) from N 2 0. In 1970 and 1971, it was 
suggested that the emission of water and NO from a 
proposed supersonic transport fleet could 
significantly reduce stratospheric ozone via these 
reactions. 

We will now turn to the question of why the pro- 
posed supersonic aircraft might have a much larger 
effect than subsonic ones, and why the 
redistribution of chemicals in the stratosphere by 
the winds is both important and difficult to estimate 
accurately. The most efficient altitude at which an 
aircraft can fly is determined by its speed. For a 
Mach 2.4 aircraft, the optimum cruise altitude is 
~20 km (13 miles), well into the stratosphere. 
Transport of the exhaust products to the 
troposphere is slow and occurs episodically, so 
material injected at 20 km stays in the atmosphere 
for months to years. These long residence times 


give the emitted exhaust material time to affect the 
ozone as it is redistributed over the globe by winds. 

An important attribute of the chemical and 
dynamical processes linking aircraft exhaust to 
ozone concentrations is nonlinearity. In a 
nonlinear system, relatively small changes may 
result in complex and sometimes amplified 
responses. When a system consists, as the 
stratosphere does, of nonlinear chemistry coupled 
with nonlinear meteorology, it is possible that the 
response of the system to perturbations such as 
aircraft exhausts will be complicated. A good 
example of a nonlinear response of a system to a 
small change is that of water just above its freezing 
point. If the temperature of the system is lowered 
only slightly, the water will freeze. However, if the 
system is at room temperature, the same change in 
temperature does not cause any appreciable effect. 

During the 1970s and early 1980s, the focus moved 
away from nitrogen oxides (NO x ) emitted by 
aircraft to the effects of chlorine from CFC 
molecules released from aerosol spray cans and 
refrigeration systems and also to bromine contained 
in agricultural fumigants and fire extinguishers. 
The primary radicals of concern in these cases are 
chlorine monoxide (CIO) and bromine monoxide 
(BrO). During this period, it was appreciated that 
different free radicals are coupled by reactions 
between members of these families to produce 
species which were mixed, such as chlorine nitrate 
(CIONO2), bromine nitrate (BrON02), and nitric 
acid (HNO3). This meant that the response of the 
stratosphere to aircraft would be complex and 
would depend on the amount of chlorine in the 
stratosphere. 

In 1985, a vivid demonstration of the nonlinearity 
of the system became apparent - the ozone hole 
over Antarctica. It was seen in the historical ozone 
record only when the stratospheric chlorine content 
reached a certain value, after which the magnitude 
of the ozone loss grew rapidly. The explanation 
was provided by a series of ground-based and 
airborne observations, later confirmed and put into 
a global context by satellite observations. The 
ozone hole chemistry demonstrated that reactions 
occurring on the surfaces of aerosol and cloud 
particles consisting of mixtures of sulfuric acid 
(H2SO4), HNO3, and H2O are important in the 
stratosphere. The abundances of these species are 
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all affected directly by aircraft emissions. A second 
important feature was that the rate of the primary 
gas-phase reaction sequence destroying the ozone 
increased nonlinearly with increasing chlorine 
abundance. 

These events led investigators to examine whether 
chemistry occurring on aerosol particles could have 
important effects, not just at the low temperatures 
over Antarctica and the Arctic during winter, but 
also in mid-latitudes. In the last five years it has 
been shown, both in the laboratory and in the 
atmosphere, that reactions occurring on H 2 SO 4 
aerosols play an important role throughout the 
lower stratosphere. The effect of this chemistry is 
to reduce the calculated amount of NO x , which 
destroys ozone catalytically, and to increase HNO 3 , 
which does not. Including this chemistry in models 
made the calculated impact of aircraft exhaust upon 
the ozone layer considerably smaller than in earlier 
calculations. 

The effect of the heterogeneous chemistry is not 
limited to changes in NO x alone; the coupling 
between NO x and the chlorine and bromine 
chemistry on one hand and between the NO x and 
hydrogen oxides (HO x ) chemistry on the other is 
such that changes in NO x also induce changes in 
the free radical abundances in these other chemical 
families. The close coupling between the chemical 
families and the large influence of heterogeneous 
chemistry in the lower stratosphere suggests that 
emissions from supersonic aircraft engines can 
affect stratospheric ozone two ways. First, 
emissions of NO x can alter all the free radical 
abundances, thus directly affecting ozone loss. 
Second, emissions of water vapor and particles can 
alter the characteristics of the existing aerosols and 
cloud particles, thus affecting the heterogeneous 
chemistry and indirectly influencing ozone loss. 
These effects further underscore the complexity of 
the stratospheric response to the emissions from su- 
personic aircraft. 

Model calculations suggest that ozone should have 
decreased from the early 1 980s at least in part as a 
result of the increase in chlorine emissions from 
1945 to the present. Such a decrease is evident in 
the observed mean record of total ozone in northern 
mid-latitudes in the overhead column over the last 
30 years, which is shown in Figure 2 [Bojkov and 
Fioletov, 1995]. The ozone record in Figure 2 



Year 


Figure 2. Deviations from the average of the 
1964-1980 column ozone level for the global 
average and for middle and polar latitudes (35" - 
90'S plus 35"-90’N). Both data sets show a 
dramatic decline (about 5% for the global average) 
over the last 15 years. Solid line represents 
ground-based data (smoothed by a 1-2-1 weight 
function); dotted line represents TOMS data. Note 
that the data is plotted against two separate scales 
on the left axis. [Adapted from Bojkov and 
Fioletov, 1995] 


shows a marked reduction in the last 15 years, the 
interpretation of which is both a complicated 
exercise and an important goal. This figure shows 
many variations that occur naturally and which are 
the result of the nonlinear chemical and 
meteorological systems interacting with each other. 
The disentangling of the causes and effects in this 
time history is not easy, nor can it alone tell us 
everything that is needed to understand the system. 

This report focuses on calculating the pattern of 
ozone concentration changes arising from a 
projected fleet of HSCTs. It pays less attention to 
climatic effects and none to effects on the 
biosphere. Such factors have been dealt with at 
length in the United Nations Environmental 
Programme/World Meteorological Organization 
(UNEP/WMO) and Intergovernmental Panel on 
Climate Change (IPCC) assessments [WMO 1986, 
1989, 1990, 1992, 1995; IPCC, 1992, 1995], which 
also make it clear that the assessment of the effects 
of HSCTs must be set in a context of parallel 
effects on the ozone layer arising from CFCs, 
subsonic aircraft, and increasing amounts of 
greenhouse gases. 
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The basic path in the assessment methodology leads 
from a determination of HSCT fleet operation and 
emissions, to calculation of the atmospheric 
distribution of HSCT emissions, and finally to the 
calculation of steady-state changes in column ozone 
throughout the atmosphere. The computational tool 
primarily employed is the two-dimensional (2-D) 
model, in which the time-consuming complexity of 
the real three-dimensional (3-D) atmosphere is re- 
duced to more manageable proportions by 
averaging around latitude circles. 

An important source of information used to place 
the 2-D model results in a realistic context is the set 
of observations gathered from aircraft, balloons, 
and satellites; it is an important objective to link 
the observational perspective with the uncertainties 
on the numerical estimates. It is important to 
appreciate that because of this simplification, 2-D 
models do not adequately simulate all of the 
features of the atmosphere, either with or without 
HSCT operations. They are an advance on the one- 
dimensional (1 -D) models used in the assessments 
made 20 years ago [CIAP, 1975] in that they 
include seasonal and pole-to-pole variations rather 
than a representation of a single atmospheric 
column. The imperfections of the simulations are 
accompanied within the Atmospheric Effects of 
Stratospheric Aircraft (AESA) component of the 
National Aeronautics and Space Administration 


(NASA) High-Speed Research Program (HSRP) by 
efforts to provide a description of uncertainties, 
with the objective of providing a balanced 
assessment of the current state of predictive ability. 
A significant portion of text in the main body of 
this report is devoted to a consideration of the un- 
certainties in 2-D calculations. 

This scientific assessment seeks to describe 
research performed in the last few years and to help 
direct efforts that would further improve the ability 
to quantify and reduce uncertainties. Advances in 
the knowledge of stratospheric transport, chemistry, 
and aerosol microphysics are discussed in Chapter 
1 along with a brief discussion of observed changes 
in ozone over the last 20 years. The implication of 
these advances in understanding are discussed 
within the context of predicting the effect of 
stratospheric aircraft, and areas of uncertainty are 
underlined. Chapter 2 presents and discusses the 
basic assessment scenarios and the calculations, 
mostly by two-dimensional models, of HSCT 
impact. Chapter 3 presents an evaluation of the 
assessment based on an understanding of how the 
stratosphere works and how this knowledge is 
represented in models. Studies of sensitivity to 
effects not included in the basic assessment are 
shown. Finally, Chapter 4 summarizes the current 
knowledge of the potential impacts of HSCTs and 
discusses the implications of this knowledge, 
including uncertainties for AESA Phase II. 
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1 ADVANCES IN ATMOSPHERIC 
CONCEPTS AND SCIENTIFIC 
ASSESSMENT APPROACH 

1 .1 Observations of Stratospheric Ozone 

Observations of stratospheric ozone establish a con- 
text for conducting the assessment of the effects of 
HSCTs. Particularly important are the mid-latitude 
observations taken over the last three decades, as 
shown in Figure 2 of the Introduction. The total 
ozone column shows clear variations that are related 
to physical and chemical processes in the strato- 
sphere. Most important is the overall decline over 
the last 15 years. Observations of the altitude distri- 
bution of ozone shows that this change is greatest 
below 25 km in the lower stratosphere (Figure 3). 
This region contains a large fraction of the ozone 
column. It is a region where gas-phase chemistry, 
transport, and heterogeneous chemistry on small 
particles are all important. This decreasing trend, the 
altitudes of greatest change, and our current under- 
standing of stratospheric chemistry all suggest that at 
least some of this ozone change is caused by in- 
creases in stratospheric chlorine and bromine chemi- 
cals. 

Shortly after the large injection of sulfur into the 
stratosphere by Mt. Pinatubo in June 1991, global 
stratospheric ozone concentrations decreased signifi- 



Figure 3. Comparison of trends in the vertical distri- 
bution of ozone during the 1980s. Ozonesonde and 
Umkehr trends are those from Miller et al. [1995]. 
Trends are given in percent per decade relative to 
1980. 95% confidence limits are shown. [WMO, 
1995] 


SBUV/SBUV2 GLOBAL OZONE 50N-50S 



Figure 4. Total ozone measured by SBUV and 
SBUV(2) since January 1992 compared with the 
1980s range and average: 50*N-50*S; global ozone. 
[WMO, 1995] 


cantly. For several years following the eruption, 
values were found to be outside the range of obser- 
vations made during the 1979-1990 period (Figure 
4). Studies suggest that these ozone changes resulted 
from changes in transport, radiation, and chemistry 
caused by the enhancement of sulfate aerosol in the 
stratosphere. Because aerosols are gradually re- 
moved from the stratosphere, the large decrease in 
ozone was temporary, and within two to three years 
the ozone trends returned to values present before the 
eruption (Figure 2). 

The most pronounced ozone changes, i.e., the ozone 
hole over Antarctica, have been observed at the high 
latitudes (Figure 5). Ozone losses similar in struc- 
ture but smaller in magnitude also have been ob- 
served over the Arctic in March. These rapid 
changes result directly from chlorine and bromine 
catalytic destruction of ozone, but the reduction of 
the NO x in these regions, initiated by the presence of 
particles, is essential for such rapid ozone loss. 

Because HSCTs would fly in the lower stratosphere, 
their exhaust would interact with the processes that 
caused the observed ozone change. The exhaust of 
proposed HSCTs will contain CO 2 , H 2 O, NO x , car- 
bon monoxide (CO), sulfur dioxide (SO 2 ), and soot, 
just as the exhaust from other aircraft do. 
Understanding the effects of this exhaust on the 
stratosphere, particularly on stratospheric ozone, re- 
quires a detailed knowledge of the aircraft exhaust 
and the processes that control ozone. Advances in 
this understanding have come from atmospheric ob- 
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SOUTH POLE STATION 



Figure 5a. Comparison of the South Pole pre-de- 
pletion ozone profile in 1993 (average of four 
soundings) with the profile observed when total 
ozone reached a minimum in 1992 and 1993. 
[Adapted from Hofmann etal., 1994] 



Figure 5b. Average ozone profiles found from 
ozonesonde measurements of Edmonton in spring 
(January to April) for 1980-1982 (46 sondes), 1988- 
1991 (42 sondes), and 1993 (13 sondes). [Adapted 
from Kerr etal., 1993]. 


servations, laboratory studies, and model calculations 
supported by the HSRP. Other research programs 
participate in this effort, including NASA’s Upper 
Atmosphere Research Program (UARP) and 
Atmospheric Chemistry Modeling and Analysis 
Program (ACMAP), and programs at the National 
Oceanic and Atmospheric Administration (NOAA) 
and the National Science Foundation (NSF). 

1 2 Current Concepts of Stratospheric 
Processes 

A simple view of stratospheric circulation is that air, 
while circling the globe, enters the stratosphere from 
the tropical troposphere, slowly ascends, moves from 
the tropics to mid-latitudes, and reenters the tropo- 
sphere at middle and high latitudes a few years later 
(Figure 6). During this transit, sunlight produces 
ozone from O2 and also produces, from source gases, 
those reactive species that destroy ozone. The 
amount of ozone results from its photochemical pro- 
duction and loss, and transport from regions with dif- 
ferent ozone concentrations. In the middle strato- 
sphere above 30 km, production and loss are much 
faster than transport, so that ozone reaches a steady - 
state amount for which its production equals its loss. 
In the lower stratosphere, both chemistry and trans- 
port affect the ozone amount as transport processes 
occur on time scales comparable to those of photo- 
chemical ozone production or loss. Photochemical 
loss of ozone is much smaller than production in the 
tropics, while at high latitudes photochemical pro- 
duction is much smaller than loss. Thus, ozone 



Figure 6. Streamlines associated with Lagrangian- 
mean velocities for equinoctial conditions. 
[Dunkerton, 1978] 
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changes due to HSCTs at any location can result both 
from local photochemical changes and from changes 
in the ozone transported to that region. 

1.2.1 Advances in Understanding 

STRATOSPHERIC CHEMISTRY 

Ultraviolet sunlight and the ensuing chemistry de- 
compose some gases that enter the stratosphere into 
chemically active gases. These chemically active 
gases can be grouped into chemical families. 
Oxygen atoms and ozone together are called odd 
oxygen (O x ). The OH and HO 2 reactive species to- 
gether are called HO x . The sum of NO, NO 2 , and 
nitrogen trioxide (NO 3 ) is called NO x . The sum of 
NO x and the more stable gases such as HNO 3 , dini- 
trogen pentoxide (N 2 O 5 ), and CIONO 2 is called odd 
nitrogen (NO y ). The sum of stable hydrogen chlo- 
ride (HC 1 ), CIONO 2 , reactive CIO, and other chlo- 
rine species is called inorganic chlorine (Cly). The 
sum of reservoir Br 0 N 02 , reactive BrO and other 
bromine species is called inorganic bromine (Br y ). 
The partitioning, which refers to the fraction of the 
family found in each of the chemical species, is sen- 
sitive to chemical processes. The total amounts of 
NO y , Cl y , and Br y are sensitive mainly to transport. 

Stratospheric ozone is produced primarily through 
the decomposition of O 2 by solar ultraviolet light. 
Chemical removal results from chemical destruction 
by NO x , HO x , halogen (CIO and BrO), and atomic 
oxygen (O) free radicals. Prior to 1990, photochemi- 
cal models predicted that throughout the middle and 
lower stratosphere, the NO x catalytic cycle: 

NO + O 3 -> NO 2 + O 2 
O + NO 2 ->02 + NO 


Net: O + O 3 — > 20 2 

and other N0 X catalytic cycles would dominate the 
photochemical destruction of ozone from below 20 
km to ~40 km altitude. If this was the case, then an 
increase in NO x would always result in a decrease in 
ozone. 

Recent free radical measurements in all chemical 
families demonstrate a different picture. Above 26 
km, where somewhat less than half the ozone column 
resides, the catalytic destruction of ozone can result 
mainly from the N0 X catalytic cycle [Jucks et al., 
1995]. However, in the lower stratosphere below 20 


km, the chemistry becomes more complex, the 
overall ozone removal rate decreases, and the frac- 
tion of the catalytic removal by NO x becomes 
smaller. In the air masses observed from NASA’s 
ER-2 high-altitude research aircraft in May 1993 
during the Stratospheric Photochemistry, Aerosols, 
and Dynamics Expedition (SPADE), ozone destruc- 
tion was dominated by reactions involving HO x , 
CIO, and BrO (Figure 7) [Wennberg et al., 1994]. 

Absolute ozone removal rates depend on chemical 
reactions that couple NO x with HO x , C10 x , and 
BrO x . Observations show that the concentrations of 
HO 2 and CIO are inversely correlated with the con- 
centration of NO x [Cohen et al., 1994; Stimpfle et 
al., 1994; Froidevaux et al., 1994]. The inverse cor- 
relation is consistent with the understanding of the 
chemical mechanisms based on laboratory studies. 
Reaction with NO converts HO 2 to OH; reaction 
with NO 2 converts CIO and BrO to CIONO 2 and 
B 1 ONO 2 . In each case, the NO x reactions suppress 


• NO x a CIO & BrO * HO x 



Figure 7. Photochemical removal rates for ozone. 
Measurements of the concentrations of OH, HO 2 , 
CIO, NO, and ozone obtained on 1 May 1 993 are 
used to infer the ozone removal rates. The HO x rad- 
icals were responsible for more than 40% of the 
photochemical loss in the lower stratosphere. The 
catalytic action of the halogen radicals accounted for 
nearly one-third of the total. [Adapted from 
Wennberg etal, 1994] 
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the rate of ozone catalysis; thus, NO x is critical in 
determining free radical concentrations and the re- 
moval rate of ozone. 

Increases in NO x , such as from aircraft exhaust, can 
cause either increases or decreases in the local ozone 
chemical loss rate (Figure 8 ) [Wennberg et al., 
1994], In regions where HO x and halogen free radi- 
cals dominate ozone loss, increasing concentrations 
of NO x decreases the concentrations of HO x and 
halogen radicals, and thus the ozone loss rates by 
HO x and halogen radicals. For regions of suffi- 
ciently large NO x concentrations, NO x catalysis 
more than compensates for the decreases in the HO x 
and halogen radical catalysis. In these regions, the 
local chemical ozone loss increases with increasing 
NO x . In the limit where NO x is driven toward zero, 
such as can happen in the winter polar stratosphere, 
ozone loss rates increase by one to two orders of 
magnitude because nearly all the available chlorine 
exists as CIO. Thus the preexisting concentrations of 
CIO, BrO, HO x , and, in particular, NO x determine 
the sign and magnitude of changes in local ozone 
loss rates that will occur when exhaust NO x is added. 

Photochemical reactions couple NO x to other mem- 
bers of the nitrogen chemical family, NO y , such as 



Increasing NO x 


Figure 8. Schematic of the ozone removal rate vs. 
NOx- Because of the coupling that exists between 
the radical families, the response of the total ozone 
removal rate to changes in NO x is highly nonlinear. 
At sufficiently low NO x such as that observed during 
the SPADE campaign, the removal rates are in- 
versely correlated with NO x . [Adapted from 
Wennberg eta!., 1994] 
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Figure 9. Partitioning of nitrogen oxide species 
between 18 and 45 km. Symbols represent mea- 
surements made by the Atmospheric Trace Molecule 
Spectroscopy (ATMOS) instrument during the 
ATLAS-3 Space Shuttle mission. The measure- 
ments were made at northern mid-latitudes (40° - 
48°N) during sunset in the first days of the mission 
(4-7 November 1994). The lines represent the re- 
sults of a zero-dimensional photochemical model 
constrained by simultaneous measurements of HCI 
+ CINO3, HNO3 + NO + NO2 + HNO4 + 2N2O5 + 
CINO3, CH4, O3, H2O, C2H6, CO, and temperature 
from ATMOS and aerosol surface area measure- 
ments from SAGE II, assuming a balance between 
production and loss for each calculated species over 
a 24-hour period. Model C (solid lines) uses rates 
for the reactions of Cl + CH4 and OH + HCI from a 
reanalysis of laboratory data and allows for a 7% 
yield of HCI from the reaction of CIO + OH. 
[Adapted from Michelsen etal., 1995] 


HNO 3 , N 2 O 5 , and CIONO 2 within hours to weeks. 
The fraction of N O v that exists in reactive form 
(NO x /NO y ) depends on the amount of solar ultravio- 
let light (and thus altitude, season, and latitude) and 
on the surface area of the aerosol particles that pro- 
vide the sites for heterogeneous chemical reactions. 
In the sunlit mid-latitudes, the partitioning of NO y 
has two distinct regimes (Figure 9). Above -30 km, 
NO x dominates the budget of NO y ; while below, 
NO x becomes a much smaller fraction [Michelsen et 
al., 1995], Figure 9 also shows comparisons to a 
photochemical model which are discussed in Section 
3.2.3. 1. 

Over the last few years, laboratory studies and atmo- 
spheric observations have led to identification of new 
mechanisms that partition NO y . Prior to 1990, the 


8 


Approx, altitude (km) 



1 

n 2 o 5 

C10N0 2 

BrON0 2 

HN0 4 

hno 3 

i 

® 

i 

V 





Sedimentation of 
HN0 3 -containing 

particles 


Chemical ozone destruction at higher altitudes allows 
ultraviolet sunlight to penetrate further into the 
stratosphere. This increased amount of ultraviolet 
sunlight produces more ozone at these lower 
altitudes. This effect, called “self-healing,” partially 
compensates for the chemical destruction of ozone, 
but the net impact still may be ozone loss. This 
coupling between ozone at different levels means 
that changes in the local ozone production and loss 
rates depend on ozone changes at higher altitudes. 

These results illustrate several features of the poten- 
tial effects of HSCT exhaust on stratospheric chem- 
istry: 


Figure 10. Processes connecting the various NOy 
species and NO x . a) The picture circa 1990, the 
easily-photolyzed reservoirs of NO x , such as N2O5, 
CIONO2, BrONC>2, and HNO4, are connected to the 
much more stable HNO3 only via gas-phase reac- 
tions involving NO x ; b) Recent laboratory and atmo- 
spheric measurements demonstrate that heteroge- 
neous reactions directly convert N2O5, CIONO2, 
and Br0N02 to HNO3 [Mozurkewich and Calvert, 
1988; Van Doren et al., 1991, Tolbert et al., 1988; 
Hanson and Ravishankara, 1993, 1994, 1995; 
Molina etal., 1987; Williams etal., 1994], 


view was that NO x radicals interconverted with the 
easily photolyzed NO y reservoirs such as N 2 O 5 , 
CIONO 2 , B 1 ONO 2 , peroxynitric acid (HNO 4 ), and 
with the much more stable HNO 3 , but that no reac- 
tions converted one reservoir species to another 
(Figure 10a). It is now recognized that heteroge- 
neous reactions on particle surfaces existing in the 
lower stratosphere convert N 2 O 5 , CIONO 2 , and 
B 1 ONO 2 directly to HNO 3 (Figure 10b). These 
reactions indirectly reduce NO x concentrations by 
reducing the abundance of easily photolyzed 
reservoirs and accumulating NO y as HNO 3 . 
Stratospheric observations in the mid-latitudes 
provide evidence supporting the conclusions of 
laboratory studies for N 2 O 5 to HNO 3 conversion 
[McElroy et al., 1992; Fahey et al., 1993; Koike et 
al., 1994]. An important sink of NO x is 
sedimentation of particles containing HNO 3 . This 
process apparently occurs only in polar regions dur- 
ing winter [Fahey et al., 1990; Santee et al., 1995], 
The resulting decrease of NO x affects the chemistry 
of the polar vortex, and may have some global con- 
sequences when the affected polar air mixes with 
mid-latitude stratospheric air. 


• NO x catalytically destroys stratospheric ozone 
but also inhibits catalytic ozone destruction by 
the HO x and halogen radical catalytic cycles 
through the formation of more stable gases. In 
the stratosphere above 26 km, NO x catalysis is 
the most important ozone loss process; increased 
NO x leads to increased local ozone loss. In the 
lower stratosphere, NO x interference with the 
dominant HO x and halogen catalysis of ozone 
destruction is more important. For increased 
NO x (e.g., from HSCT exhaust), the sign and 
magnitude of the local chemical ozone change in 
the lower stratosphere depend on the meteorolog- 
ical conditions and the background amounts of 
trace gases and aerosol surface area. 

• HSCT exhaust would increase the water vapor 
and sulfur in the lower stratosphere and indirectly 
influence ozone loss by increasing the amount of 
sulfate aerosol. An increase in the amount of sul- 
fate aerosol will lower NO x and thus reduce the 
direct impact of aircraft-generated NO x on ozone. 
The aerosol increase also will lead to increased 
catalysis of ozone destruction by HO x and halo- 
gens. The net effect of the HSCT exhaust will 
depend on magnitude of the aerosol surface area 
changes and the impact of these changes on each 
of the chemical destruction cycles for ozone. 

Local chemical production and loss is not the only 
means by which ozone is changed. The ozone con- 
centration in a region of the atmosphere is also sen- 
sitive to transport of air from adjoining regions. 
Changes in chemical loss rates by catalysis do not 
necessarily result in proportional changes in ozone. 
The ozone concentration in a region is sensitive, not 
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only to the chemical losses that can occur in that re- 
gion, but also to the magnitude of the perturbations 
which have occurred elsewhere in the stratosphere 
and then propagated to that region by transport. 

1.2.2 Advances in U nderstanding 
Stratospheric Transport 

A calculation of zonal-mean mass transport in the 
stratosphere is shown in Figure 6. The exchange of 
mass between the tropics and mid-latitudes appears 
to be restricted at certain altitudes, and may include 
some transport from the mid-latitudes back into the 
tropics. In the fall, the rapid cooling and formation 
of a strong polar wind jet cause polar air to subside 
many kilometers and leads to the development of a 
transport barrier between the polar regions and mid- 
latitudes. Transport of stratospheric air across the 
tropopause, which determines in part the mean resi- 
dence time of air in the stratosphere, occurs sporadi- 
cally when stratospheric air gets folded into the tro- 
posphere. The paths of the air parcels and their rate 
of mixing with other air parcels are important for un- 
derstanding the response of the stratosphere to HSCT 
exhaust. 

Transport and the dispersal of HSCT exhaust are key 
issues for understanding the potential effects on 
stratospheric ozone. An important aspect is that 
fraction of HSCT exhaust which would be trans- 
ported from the heavily traveled flight corridors in 
mid-latitudes to the tropics where it can be lofted to 
higher altitudes. Subsequent transport back to mid- 
latitudes would be at altitudes where the NO x from 
the exhaust causes increased ozone loss. The degree 
of isolation between the winter polar vortices and 
mid-latitudes determines the extent to which chem- 
istry occurring in the vortices, and any change in that 
chemistry resulting from HSCT emissions, will af- 
fect global ozone levels. Thus, the deposition of 
HSCT exhaust in the lower stratosphere becomes a 
global issue because of transport. 


satellites. Mixing ratios of stratospheric constituents 
fluctuate at a fixed location due to reversible 
horizontal and vertical displacements, but examining 
changes in one tracer relative to another effectively 
removes this variability [Ehhalt et al., 1983]. 
Observations of compact, simple relationships be- 
tween tracers in the lower stratosphere [Fahey et al., 
1989; Woodbridge et al., 1995] (Figure 11) suggest 
that the surfaces of constant mixing ratio, or 
“isopleths,” are similar for these species if quasi-hor- 
izontal mixing is faster than their local chemical life- 
times [Holton, 1986; Mahlman et al., 1986; Plumb 
and Ko, 1992]. 

If the observed relationships hold from pole to pole, 
then the global stratospheric tracer transport would 
be well represented by the “global diffuser” model 
(Figure 12a). In this model, transport occurs by the 


O CFC-12 
+ CCI, 


VP* 0 




*9 




4 » 1 1 «m ^ i 





0^*0 o ♦ ♦ 

aTP V f 

■ 

+ 

8 

m o 0 

a o 

o o °0 0 


eg a ,ii imuiif***** 


150 200 250 

N 2 0 (ppbv) 


Recent observations have led to new insights into 
stratospheric dynamics, and a clearer conceptual 
picture of transport is emerging. The most powerful 
of the observations involve simultaneous measure- 
ments of suites of long-lived gases, called tracers, 
with contrasting sources, sinks, and lifetimes in the 
stratosphere (e.g., N2O, NOy, CFCs, CO 2 , H 2 O, 
CH4, CO, ozone, and aerosols). These measure- 
ments have been obtained from aircraft, balloon, and 


Figure 11. Concentrations of halocarbons in the 
lower stratosphere from Whole Air Sampler vs. N 2 O 
from the Airborne Tunable Laser Absorption System. 
These relationships have been discussed for NOy 
and N 2 O by Fahey et al. [1990]. Their theoretical 
basis is discussed by Holton [1986], Mahlman ef al. 
[1986], and Plumb and Ko [1992]. [Adapted from 
Woodbridge eta!., 1995] 
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Figure 12. Idealized models for global-scale 
transport of tracers, a) Global diffuser model, with 
rapid exchange between mid-latitudes and tropics 
due to quasi-isentropic mixing in the winter mid- 
latitude surf zone, b) Tropical pipe model, with 
vertical advective transport in the tropics, 
entrainment of tropical air into mid-latitudes, and 
quasi-isentropic mixing restricted to subtropical 
latitudes and higher. [Adapted from Plumb, 1995] 


simple, mean motions depicted in Figure 6 and rapid 
mixing from pole to pole along quasi-horizontal sur- 
faces. The result is globally uniform mixing surfaces 
and therefore globally uniform tracer-tracer correla- 
tions. 

However, observations showed confinement of Mt. 
Pinatubo and background aerosols in the tropics 
[Trepte and Hitchman, 1992; Trepte et al., 1993] 
(Figure 13); tracer gradients across the subtropics 
[Randel et al ., 1993]; and gradients in ozone, NOy, 
and the ozone/NO y ratio between the tropics and 
mid-latitudes in the lower stratosphere [Murphy et 
al ., 1993]. These observations led to the idea of a 
subtropical “barrier” to transport of air from mid- lati- 
tudes into the tropics and to the development of an- 
other conceptual picture - the “tropical pipe” model 
(Figure 12b). This model is similar to the global 
diffuser model except that air can travel across the 
“subtropical barrier” only from the tropics to mid- 
latitudes, resulting in tracer-tracer correlations that 
are significantly different in the tropics than in mid- 
latitudes. 


The atmospheric effects of HSCTs are different in 
the two concepts. If only the global diffuser model is 
applicable, HSCT exhaust would freely enter the 
tropics and be lofted to high altitudes where it would 
lead to enhanced ozone loss. If only the tropical pipe 
model is applicable and all upwelling is confined to 
the tropics, then the middle stratosphere would be 
largely protected from exhaust products released 
outside the tropics; changes to atmospheric composi- 
tion would be largely confined to low altitudes and 
mid- to high latitudes. 

Consistent with aspects of both models, aircraft and 
satellite observations indicate that air enters the 
stratosphere in the tropics and indeed is transported 
rapidly poleward into mid-latitudes up to -20 km 
[McCormick et al ., 1993; Trepte et al 1993; 
Boering et al., 1994, 1995]. The observed propaga- 
tion of seasonal variations in CO 2 and H 2 O implies a 
vertical ascent rate in the tropics of -20 meters per 
day, comparable to that expected from calculated 
heating rates, and transport to mid-latitudes on time 
scales as short as one month. 

Likewise, the isolation of the tropics from mid-lati- 
tude air is not complete and may have seasonal and 
interannual variations. Analyses of tropical tracer 
correlations demonstrate that 25 to 40% of air in the 
tropics at 20 km was transported from mid-latitudes 
[Avallone and Prather, 1995; Minschwaner et al, 
1995]. Thus, the tropical pipe is “leaky,” though in- 
dications are that the pipe becomes less leaky and the 
tropics become more isolated at altitudes above 20 
km [Michelsen et al, 1995; Mote et al . , 1995]. 

Another critical issue is the residence time of the air- 
craft emissions, which is determined by the rate at 
which the emissions are flushed out of the lower 
stratosphere into the troposphere. The final stage of 
transport across the tropopause is effected through 
“tropopause folding” events associated with active 
upper tropospheric synoptic weather systems. Until 
a few years ago, it was thought that these processes 
controlled the overall rate of stratosphere-tropo- 
sphere exchange, but it is now recognized that the 
rate of mass exchange is controlled by the large-scale 
stratospheric circulation [Haynes et al, 1991; Holton 
et al. , 1995]. However, synoptic-scale folding events 
may be a controlling influence on the residence time 
of exhaust gases injected into the very low mid-lati- 
tude stratosphere (below -16 km). 
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SAGE II Aerosol Observations 
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Figure 13. Latitude-altitude cross sections of the SAGE I1 1- pm extinction ratio measurements that show the 
effect of the Mt. Pinatubo eruption in June 1991 on aerosol abundance in the lower atmosphere. The specific 
observation dates are indicated with crosses below each panel for the periods: a) 15 April to 25 May 1991 
(pre-eruption); b) 14 June to 26 July 1991 (early austral winter); d) 20 August to 30 September 1991 (late 
austral winter); and, i) 29 March to 9 May 1992 (full dispersal). No data were used 2 km below the tropopause 
(blacked out). Small triangles indicate truncation altitude for the SAGE II data. Lidar data were used below 
this altitude. Isentropes (constant potential temperature in K) appear as white contour lines. Similar concepts 
come from observations of UARS data [Randel et al., 1993] and NO y and ozone by ER-2 instruments in the 
lower stratosphere [Murphy et al., 1993], [Adapted from Trepte eta!., 1993] 


Understanding the transport of air into and out of the 
polar vortex is also important, since chemical and 
microphysical processes which occur solely in the 
vortices may be affected by HSCT exhaust. Recent 
observations from balloons, aircraft, and satellites 
suggest that above -16 km the polar vortices are 
largely isolated from exchange with the mid-latitudes 
[Vomel et al., 1995; Mote, 1995; Waugh et al.. 


1995]. Concentrations of trace species at mid-lati- 
tudes are not dramatically affected during winter by 
transport out of the vortex. However, transport of 
material out of the vortex in spring and exchange of 
material near the vortex edge may be important. 

Transport controls the dispersal of the HSCT exhaust 
and thus the concentration of its constituents and its 
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distribution throughout the stratosphere. Because the 
abundance of the exhaust constituents and the meteo- 
rological conditions of their location dictate the 
chemical reactions that will occur, transport and 
chemistry both determine the potential effects of 
HSCT exhaust. 

1.2.3 Advances in Understanding Aerosol 
Microphysics 

HSCT exhaust emits condensible gases and particles 
(aerosols) into the atmosphere. Sulfate aerosols re- 
sult from the conversion of fuel sulfur to sulfur ox- 
ides (SO x ) in the engine combustor, followed by 
conversion to H2SO4 in the expanding exhaust 
plume or in the free atmosphere and, in turn, growth 
of preexisting particles or nucleation of new parti- 
cles. Soot aerosols are a product of incomplete fuel 
combustion. 

Increases in atmospheric sulfate and soot aerosols 
could affect both climate and chemistry. Sulfate 
aerosols are primarily scatterers of incoming solar 
radiation, while soot aerosols are primarily ab- 
sorbers. Therefore, increases in sulfate and soot may 
affect the radiative forcing of the atmosphere-Earth 
system. These aerosols may also impact the forma- 
tion, and thus the radiative effects, of cirrus clouds. 
From a chemistry perspective, increased surface area 
from HSCT exhaust aerosols and condensible gases 
could change the rates of heterogeneous reactions 
that alter the stratospheric NO y partitioning, and con- 
sequently the HO x , Cly, and Br y partitioning. 

Advances have occurred in the quantification of air- 
craft-induced aerosol effects. Measurements show 
that both soot and sulfate from HSCT exhaust could 
increase the abundance of small particles in the lower 
stratosphere. Soot emissions could double or triple 
the amount of natural stratospheric soot. Because the 
natural soot loading is only 1% that of sulfate, the 
potential impact of this increase is unclear. Other 
observations have shown that significant numbers of 
persistent small sulfate particles form in the exhaust 
plume of supersonic aircraft, with appreciable total 
surface area available for catalyzing heterogeneous 
reactions [Fahey et al., 1995 a]. Sulfur-soot interac- 
tions in the near-field region behind the aircraft are 
possible, such as the activation of soot particles as 
nuclei for condensation of H2O and H2SO4. 


There is a greater understanding of the complex mi- 
crophysics that occur in the low-temperature atmo- 
spheric regions such as the Arctic and Antarctic. 
Observations of polar stratospheric clouds (PSCs) 
show considerable variability that is consistent with 
the complex ternary phase diagram of HNO3, H2O, 
and H2SO4 rather than the simpler model of equilib- 
rium nitric acid trihydrate (NAT) thermodynamics. 
For temperatures above 195 K, the behavior of 
stratospheric aerosols is consistent with liquid sulfate 
composition in the range of 50 to 80 % H2SO4 by 
weight. At temperatures of -192 K, both HNO3 and 
H2O can condense to form ternary solution droplets 
whose composition depends on the ambient 
HNO3/H2SO4 ratio. The existence of supercooled 
liquid particles at these temperatures is consistent 
with earlier in situ and remote sensing observations 
and with laboratory and theoretical studies that pre- 
dict significant thermodynamic barriers to freezing of 
both binary H2SO4-H2O and ternary H2SO4-H2O- 
HNO3 solutions. Other laboratory studies show 
multiple metastable solid phases are possible when 
the solution droplets freeze. Understanding the de- 
tailed temperature-dependent behavior of these parti- 
cles is a key in determining their efficiency in cat- 
alyzing heterogeneous reactions and in triggering 
denitrification, the irreversible removal of NO y in 
polar regions. The mechanisms underlying the steep 
temperature dependence of both microphysical and 
reaction kinetic processes must be accurately coupled 
in order to model high-latitude, cold-temperature 
photochemical cycles. 

1.3 Advances in HSCT Emissions 

Characterization: Measurements in the 

Plume of the Concorde 

During the Airborne Southern Hemisphere Ozone 
Experiment/Measurements for Assessing the Effects 
of Stratospheric Aircraft (ASHOE/MAESA) ER -2 
campaign in 1994 , measurements of NO, NO y , CO2, 
H2O, condensation nuclei (CN), N2O, OH, and HO2 
were made in the exhaust plume of a Concorde su- 
personic aircraft while in flight (Figure 14 ) [Fahey et 
al, 1995 a], The emission indices (Els) were calcu- 
lated independently of plume dynamics because 
measurements of CO2 provide a measure of fuel 
burned. The observations show: 

• The measured El of NO x , 23 ±5 grams as equiva- 
lent NO2 per kilogram of fuel burned, is in 
agreement with the earlier results deduced from 
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Time (UT s) 

Figure 14. Emissions from an Air France Concorde 
measured in situ from the ER-2 on 8 October 1 994. 
The bottom axis indicates Universal Time (UT). 
[Adapted from Fahey etal., 1995a] 


altitude chamber measurements as part of the 
Climatic Impact Assessment Program (CIAP). 
These results suggest that the methodology de- 
veloped to extrapolate the El of NO x at cruise 
from altitude chamber data will be appropriate to 
evaluate new HSCT engines. 

• No enhancements in CO above stratospheric 
background levels were observed in the plume. 
This limits the El of CO to be less than 3.5 g 
CO/kg fuel. Provided the combustion efficiency 


of the future HSCT fleet is no lower than that of 
the Concorde, CO emissions will be negligible. 

• The measurements of OH in the plume suggest 
that hydroxyl emissions are approximately 3 to 
5% of the NO x emissions and that the OH reacts 
quickly with NO to form nitrous acid (HONO) 
[Arnold et al., 1992; Hanisco et al . , 1995]. 
During daylight, the regeneration of OH through 
HONO photolysis increases OH above back- 
ground values. The measured OH values suggest 
that only a small fraction of NO x (5%) and SO2 
(1%) are oxidized by OH in the plume. 
Heterogeneous chemistry involving HONO may 
be important, particularly at night [Zhang et al., 
1995]. 

• The emission index for total particles as mea- 
sured with a CN counter exceeded 10 17 (kg 
fuel) -1 . Approximately 70 to 90% of these parti- 
cles were volatile when heated to 192°C. 
Laboratory studies indicate that the volatile par- 
ticles are likely sulfuric acid and water mixtures 
and not soot. The lower size-limit of the CN 
counter (0.009 nm diameter) implies that between 
12 and 45% of the fuel sulfur is oxidized to form 
H2SO4 [Fahey et al . , 1995a]. The OH radical is 
generally assumed to be responsible for sulfur 
oxidation; however, the level of OH measured in 
the Concorde plume is consistent with no more 
than 1% oxidation. These results suggest that 
sulfur oxidation by OH is not the dominant pro- 
cess in the plume. A suitable mechanism ratio- 
nalizing the observed degree of oxidation has yet 
to be identified. 

The Concorde plume measurements highlight the is- 
sue of particulate (i.e., soot and sulfur) emissions 
from HSCTs. The potential for the aircraft emissions 
of H2O, sulfur, and particulates to change the 
NO x /NOy ratio in both the polar and mid-latitudes 
may be more significant than the emissions of NO x . 
These results heighten the requirement that 
characterization and understanding of the potential 
changes in aerosol chemistry be appreciated. 

1.4 HSCT Perturbations: The Assessment 
Approach 

The AES A component of NASA’s HSRP is designed 
to assess the potential atmospheric effects of a fleet 
of HSCT aircraft operating in the stratosphere. As in 
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Figure 15 . Schematic of the HSCT Scientific Assessment Process. 


previous assessments, the current effort emphasizes 
changes in ozone that could result from HSCT opera- 
tions, although preliminary studies of changes in 
climate are also included. As outlined in previous 
AESA reports [Prather and Wesoky, 1992; Stolarski 
and Wesoky, 1993a, b, 1995; Prather and Remsberg, 
1993; NASA, 1993], significant resources have been 
devoted to carefully handling the details of the as- 
sessment process. 

This document describes predictions of the steady- 
state ozone changes that might result from HSCT op- 
eration (Figure 15). The predictions cannot be ob- 
tained directly from observations or the concepts 
outlined in this chapter; they must come from photo- 
chemical transport model calculations. 2-D photo- 
chemical transport models have been chosen for this 
task. 

The end result of the assessment process is not just 
the predicted column ozone change or changes in the 
climate-related chemical species, temperature, and 
atmospheric dynamics. As important is the statement 
of the uncertainty in the prediction. The discussion 
of uncertainties provides the greatest information 
about confidence in the prediction of ozone change 


and about advancement in the understanding of 
stratospheric processes and ozone. 

The assessment is based on a large array of tools, of 
which the 2-D models are one component. Other 
tools used in the assessment include atmospheric 
observations, laboratory studies, HSCT emissions 
inventories, emissions characterization, and a variety 
of photochemical transport and process models. 
Atmospheric observations supported under AESA 
allow many key processes and constituents to be 
observed directly in the lower stratosphere using 
specialized instruments on airborne platforms. 
Laboratory studies have expanded the number of re- 
actions, particularly heterogeneous reactions, that 
must be considered in the models and have provided 
measured rate coefficients. Ground-based data sets 
indicate significant trends in long-lived source gases 
that are produced in the troposphere and transported 
to the stratosphere. 

An essential part of this process is the comparison of 
models and observations. These provide the basic 
understanding of stratospheric processes, the devel- 
opment of conceptual models, and means to test the 
quality of model performance. Observations have 
been most usefully compared with process models 
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that can mimic the meteorological and chemical con- 
ditions of the individual observations. However, 2-D 
and 3-D model results have been compared only in 
limited cases. Of the advances in understanding out- 
lined above, only some have been included in the 
2-D predictive models. Sensitivity tests have been 
performed for the aerosol effects and transport, but 
the advances in transport and aerosol effects have not 
yet been fully implemented. 

Particles formed in HSCT exhaust plumes may have 
an important role in changing global stratospheric 
aerosol abundances and in increasing the rates of 
heterogeneous reactive processes. Climatic effects 
by an HSCT fleet have been recognized as a poten- 
tial concern, but appear to be small compared to 
other anthropogenic climate perturbations. Climatic 


assessments are preliminary and depend on future 
improvements of the climate models and their appli- 
cation to this issue. The direction and emphasis of 
the AESA Phase II program will be influenced by the 
understanding of the relative importance of the ul- 
traviolet and climate effects caused by ozone 
changes. 

The atmosphere is changing. Climate is changing as 
radiatively active gases are being added to the atmo- 
sphere. Stratospheric chemistry is changing as CFCs 
and thus stratospheric chlorine are reduced. The 
subsonic aircraft fleet is expected to approximately 
double in size by the year 2015. The framework that 
is being developed for this assessment must be robust 
enough to make confident predictions in the face of 
atmospheric change. 
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2 MODEL PREDICTIONS OF HSCT 
EFFECTS FOR BASIC SCENARIOS 

An initial inventory of emissions from an HSCT fleet 
was given in [Miake-Lye et al., 1992] in the first 
AES A Program Report [Prather and Wesoky, 1992], 
which concluded that the nonexhaust emissions are 
unlikely to affect the chemical content and radiative 
balance of the atmosphere. The primary emissions 
(in term of mass emitted) from aircraft are H 2 O and 
CO 2 , which are the end products of jet fuel 
combustion. NO x (= NO + NO 2 ), CO, hydrocarbons 
(HCs), and soot also are produced in aircraft engine 
combustors and vary in quantity according to the 
combustor temperature, pressure, and design. Sulfur 
is present as an impurity in the fuel. For a given 
engine, the Els of NO x , CO, and HCs will vary with 
the power setting and, thus, will be different for 
different phases of the flight [ICAO, 1995]. For 
example, NO x levels are highest at high thrust 
(takeoff and climb) and lowest during descent. Table 
1 gives typical Els of the HSCT engines based on the 
NASA Interim AESA Assessment Report [NASA, 
1993] and on a recent Boeing report [Baughcum et 
al., 1994]. The emission characteristics of Concorde, 
as measured in an altitude chamber simulating 
supersonic cruise conditions [CIAP, 1975], are 
included for comparison. As described in the 
previous section, recent measurements in the wake of 
the Concorde flying at supersonic cruise yielded an 
El of 23 ±5 g NO x /kg fuel, which is in agreement 
with the CIAP results [Fahey etal., 1995a]. 

The Els given in Table 1 can be combined with the 
expected fuel use of an HSCT fleet to give the total 
emission into the stratosphere. The increase in 
stratospheric burden of a species can be obtained by 
using the emitted amount and an estimated residence 
time of the species in the stratosphere. Table 2 
summarizes the estimated range of perturbations 
from an HSCT fleet with supersonic cruise fuel use 
of ~50 x 10 9 kg yr 1 . The effect from the emitted 
sulfur may depend on the assumed gas-to-particle 
conversion rates in the plume (see Sections 3.2.4.3 
and 3.2.4.5). 

2.1 Description of Fleet Scenarios 

In order to assess the potential impact of a fleet of 
HSCTs, an inventory of the HSCT emissions as a 
function of latitude and altitude is needed as model 
input. Two combustor concepts are being evaluated: 


Table 1. Typical Els (in units of grams emission per 
kilogram fuel) for the Concorde and the projected 
HSCT engines at stratospheric cruise altitudes. 


Species 

Concorde 

Projected 

HSCT 

Engines 

Carbon Dioxide (CO 2 ) 

3155 

3155 

Water (H 2 0) 

1237 

1237 

CO 

3.5 

2.9 

NO x (as N0 2 ) 

18 

5-15 

Hydrocarbons (HCs) 

0.2 

0.3 

Soot (as carbon) 

0.03 

0-0.02 

Sulfur Dioxide (SO 2 ) 

0.8 

0.4 


lean, premixed, pre vaporized and two-stage rich- 
burn, quick-quench [see Shaw et al., 1993]. 
Although full-scale combustors utilizing these 
concepts have not been built, preliminary 
measurements suggest that low-NO x levels can be 
achieved. However, because of uncertainty in what 
NO x Els can be achieved in future engine designs, 
the HSCT assessment has been conducted treating 
EInOx as a parameter over the range of 5-15 g 
NO x /kg fuel. The Els for HSCT engines at 
stratospheric cruise conditions used in the 
calculations are given in Table 1. 

To calculate a 3-D inventory of HSCT emissions, it 
is necessary to determine the ground rules under 
which an HSCT would operate (e.g., no supersonic 
flight over land), identify which cities would be 
connected via HSCT flights, and calculate the 
number of flights between each city. These 
assumptions were based on a forecast of passenger 
demand in the year 2015. Then, based on the 
preliminary design of the aircraft and the engine, the 
fuel bum and emissions were calculated for each 
flight, placed onto a 3-D grid, and summed to 
produce the emission inventory for the fleet. For the 
HSCT assessments, emission inventories were 
developed on a 1° longitude by 1 ° latitude by 1 km 
altitude grid. 

For the 1993 AESA Interim Assessment [NASA, 
1993], emission inventories were developed for vari- 
ous HSCT fleets (Mach 1.6, 2.0, and 2.4), assuming 
a passenger demand corresponding to that of a fleet 
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Table 2. Catalog of potential stratospheric perturbations by HSCT fleet of 500 aircraft. 


Emittant 

Projected 

Perturbation 

Potential Atmospheric Interactions 

Ozone Effects 

NO x 

Peak 50-100% in- 
crease a 

Ozone depletion by NO x catalysis, interference with 
C10 x , HO x , and BrO x catalysis 



Conversion to HNO 3 , NAT condensation, increased 
CIO x catalysis in polar vortex 

h 2 o 

Peak 10-20% in- 
crease 

HO x formation and ozone depletion, interference with 
NO x catalysis, enhanced C10 x catalysis 



NAT/ice condensation, denitrification in polar vortex, 
increased C10 x catalysis 

Sulfur 

10 - 200 % increase b 
in surface area 

Increased aerosol surface area, enhanced ozone deple- 
tion by C10 x , decreased ozone depletion by NO x 

Soot 

Highly uncertain: 
0-300% increase 

Additional nucleation sites for aerosols, increased sur- 
face area-changes in chemical reactivity of aerosols 

Hydrocarbons 

~ 0 . 1 % increase 
compared to CH 4 

Source of CO, HO x , H 2 0 

CO 

5-20% increase 

Modification of catalysis by HO x and NO x 

Radiative Forcing Effects 

C0 2 

Current subsonic 
-3% of C0 2 from 
fossil fuel - HSCT 
- 1 % 

Direct change in IR radiative forcing 

h 2 o 

Peak 10-20% in- 
crease 

NAT/ice condensation, cirrus cloud formation, change 
in radiative forcing 



Direct change in IR radiative forcing 

Sulfur 

10 - 200 % increase b 
in surface area 

Increased aerosol mass loading, change in radiative 
forcing 

Soot 

Highly uncertain: 
0-300% increase 

Additional nucleation sites for aerosols, increased sur- 
face-area change in radiative forcing 

NO x 

Peak 50-100% in- 
crease 

Ozone depletion, change in radiative forcing 


a Estimated for EIno x = 5 

b Depends on assumptions concerning gas to particle conversion in the plume 


of 500 Mach 2.4 aircraft. Simple ground rules were 
used which assumed that the HSCT would capture a 
constant fraction of the passenger flow between the 
candidate cities [Wuebbles etal., 1993; Baughcum et 
al., 1994; Landau et al . , 1994], Subsonic air traffic 
emissions were also evaluated and displacement of 


subsonic aircraft by HSCTs was considered 
explicitly. 

Since some routes are more profitable for airlines 
than others and since the biggest benefit of the HSCT 
for passengers is the shorter travel time, an HSCT 
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marketing analysis was completed to take these 
factors into account. Updated emission scenarios for 
Mach 2.0 and 2.4 HSCTs were calculated using a 
Boeing projection of year 2015 passenger demand 
[Baughcum and Henderson, 1995]. HSCTs were 
scheduled assuming a universal (i.e., single, world- 
wide) airline. The fleet sizes calculated were for a 
very active, high-utilization HSCT fleet and should 
represent an upper limit of the amount of emissions 
generated by that number of aircraft. Revised 
subsonic aircraft scenarios were also developed to be 
consistent with this HSCT fleet. The new scenarios 
predict more emissions in the Southern Hemisphere 
and in the tropics than did the 1993 Interim 
Assessment scenarios. The total estimated fuel usage 
for a 500-aircraft, Mach 2.4 fleet in these scenarios is 
82 x 10 9 kg per year. Figure 16 shows the 
percentage of fuel burned at altitudes above 1 3 km as 
a function of latitude for a fleet of 500 Mach 2.4 
HSCTs. Figure 17 shows the fuel burned as a 
function of altitude for an all-subsonic fleet and for a 
modified subsonic fleet plus 500 Mach 2.4 HSCTs in 
2015. The NO x emissions at supersonic cruise alti- 
tudes are assumed to be proportional to fuel burned 
for a constant supersonic cruise EInOx- 

For the following assessment model runs, the Mach 
2.4 HSCT with EInOx = 5 at stratospheric cruise has 
been chosen as the standard scenario. Such an 
aircraft would cruise in the 18- to 21 -km altitude 
range. For an assumption of constant HSCT 
passenger demand, the studies have shown the 
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Figure 16. Latitudinal distribution of fuel used above 
13 km altitude by a Mach 2.4 HSCT fleet. 



Figure 17. Total estimated fuel burn as a function of 
altitude for an all-subsonic 2015 fleet and for a 2015 
fleet which includes a modified subsonic fleet plus 
500 Mach 2.4 HSCTs. 


emission levels for a Mach 2.0 fleet to be comparable 
to the Mach 2.4 fleet shifted to lower altitude by 2 
km [Baughcum et al., 1994; Baughcum and 
Henderson, 1995]. Because of this, the altitude 
sensitivity has been evaluated in the model scenarios 
by shifting the Mach 2.4 HSCT emissions by +1 km 
and -2 km. 

In order to evaluate how growth in the HSCT fleet 
size would affect the geographical distribution of 
engine exhaust, emission scenarios were calculated 
by both Boeing [Baughcum and Henderson, 1995] 
and McDonnell Douglas [Metwally, 1995] for fleets 
of -1000 HSCTs. (As a reference point, -1000 
Boeing 747s have been produced over the last 25 
years.) The Boeing study was based on a universal 
airline fleet of 991 active HSCTs. The McDonnell 
Douglas study used four regional airlines (European, 
Far East, Americas-east, and Americas-west) to 
schedule departures of 1059 HSCTs. Both studies 
used the same baseline aircraft and emission 
characteristics, but different assumptions about 
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market penetration and utilization. Thus, while they 
generated similar amounts of global fuel burned and 
emissions, some differences in the geographical 
distribution were predicted. The market penetration 
model determines the fraction of the HSCT 
passenger demand between two cities; hence, 
differences in geographical distribution for a given 
total passenger demand are due largely to the 
different market penetration assumptions. For the 
same model, changes in airline scheduling (e.g., a 
universal airline vs. four regional airlines) would be 
expected to modify somewhat the number of HSCTs 
required to meet the passenger demand but not the 
total emissions deposited within the stratosphere. 

To evaluate the effects of doubling the HSCT fleet 
size, two model runs (EIno x = 5 and 15) were made 
using the Boeing study. The emissions from the 
McDonnell Douglas study are sufficiently similar 
that it was deemed unnecessary to perform a 
different set of scenarios. The total fuel burned for 
the HSCT component of this scenario was 157 x 109 
kg per year, with ~63% between 17 and 21 km. 

The 2-D models used the 3-D emission data by 
summing the emissions within a latitude belt and in- 
troducing them into the model grid box representing 
that latitude belt. This procedure assumes that the 
emitted materials are instantly mixed zonally. 
Previous model studies showed that the effects on the 
calculated ozone changes from adjusting the source 
function for subsidence of the plume [Rodriguez et 
al., 1994] and for the nonzonal excursion of the 
plume before it is zonally mixed [Plumb et al, 1995] 
are small. The annual amount was put into the 
model at a uniform rate through the year without 
accounting for possible seasonal variations in HSCT 
traffic. 

2.2 Description of 2-D Models Used in the 
Assessment 

For the practical purpose of accomplishing long-term 
simulations for an assessment such as this, approxi- 
mations are made to reduce computing time. One 
approach that has proven to be practical is the use of 
2-D models that solve systems of equations for 
zonal-mean quantities. These 2-D models take 
advantage of the fact that atmospheric transport 
occurs rapidly around a latitude circle on a time scale 
of one to two weeks. Thus, they consider only 
longitudinally averaged or zonal-mean quantities. 


They include a relatively complete suite of chemicals 
and reactions as a function of latitude and altitude. 

The behavior (as a function of latitude, altitude, and 
time) of the zonal-mean concentrations of the key 
chemical constituents are obtained by solving a set of 
continuity equations. The time rate of change of the 
concentration of a species at a particular location is 
given in terms of the changes due to photochemical 
interactions and transport to and from neighboring 
locations. A number of inputs are required to specify 
these quantities. These include boundary conditions 
at the ground and the top of the atmosphere to 
specify the exchange rates, chemical reaction rate 
coefficients and absorption cross sections to calculate 
the photochemical rates, and the transport fields to 
specify the transport rates. 

The dynamical foundation for describing the 
transport of the zonal-mean quantities in 2-D models 
was developed in the 1960s [Murgatroyd and 
Singleton, 1961; Reed and German, 1965]. These 
ideas were refined and incorporated into numerical 
models with chemistry during the 1970s and early 
1980s [see review in Tung, 1984]. While all the 
models that participated in this assessment use some 
form of residual mean circulation calculated from 
diabatic heating rates as an approximation to the 
advective transport circulation for constituents, they 
differ in the choice of the circulation and eddy 
diffusion coefficients [see Prather and Remsberg, 
1993]. The diabatic heating rates are calculated off- 
line or taken from published results such as those of 
Dopplick [1972, 1979], Shine [1987], Rosenfield et 
al. [1987], and Briegleb [1992a, b]. The calculated 
diabatic heating rates differ from one another 
because of differences in the assumptions in the 
radiative code and differences in the background 
temperatures assumed in the calculation. The 2-D 
models also differ in their assumptions concerning 
mixing processes by large-scale eddies (planetary 
waves) and small-scale processes. Formulations are 
being developed following Garcia et al. [1992] that 
self-consistently represent eddy processes, but these 
have not yet been applied to the aircraft assessment. 

The set of equations for the different species are cou- 
pled because the local photochemical rate for a 
species depends on the concentration of other species 
at that location. Even in the 2-D approximation, 
these integrations take significant computer time. 
Thus, some of the models employ grouping schemes 
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Table 3. Models participating in the assessment. 


Model 

Investigators/Institution 

AER 

M. Ko and D. Weisenstein, 
Atmospheric and Environmental 
Research, Inc. (USA) 

CAMED 

R. Harwood, J. Kinnersley, Y. Roberts, 
J. Pyle, University of Edinburgh/ 
University of Cambridge (UK) 

CSIRO 

K Ryan, I. Plumb, P. Vohralik, and 
L. Randeniya, Commonwealth Science 
Commonwealth Scientific and 
Industrial Research Organization 
(Australia) 

GSFC 

C. Jackman, D. Considine, and 
E. Fleming, NASA Goddard Space 
Flight Center (USA) 

LLNL/UI 

D. Kinnison and P. Connell, Lawrence 
Livermore National Laboratory/ 

D. Wuebbles and K. Patten, University 
of Illinois (USA) 


for constituents to reduce the number of equations. 
Several constituents such as those in the NO x family 
are grouped together and the continuity equation 
with transport terms is solved for the family. The 
chemical partitioning occurring within the family is 
calculated under the assumption that a steady state is 
reached. 

The 2-D models that provided results for this assess- 
ment are identified in Table 3. With the exception of 
the Commonwealth Scientific and Industrial 
Research Organization (CSIRO), the models 
provided results for the 1993 Interim Assessment 
Report. Descriptions of these four models were 
given in that report [NASA, 1993]. However, these 
models have been modified since then, not only in 
the chemistry used but also in their treatment of 
atmospheric processes. These models will continue 
to evolve as our knowledge of atmospheric processes 
improves. It is inappropriate to ascribe the 
differences in ozone predictions in this and previous 
reports solely to changes in the models, since 
significant changes have been made in the boundary 
conditions chosen for the assessment and in the 
emission scenarios. Brief descriptions of the model 
changes, together with a full description of the 
CSIRO model, are given in Appendix A. 


Heterogeneous processes on background sulfate 
aerosols have been shown to play an important role 
in controlling the chemistry of NO x , HO x , C10 x , and 
BiO x which, in turn, control the concentration of 
ozone in the stratosphere [WMO, 1995]. The effect 
of such reactions is to reduce the role of NO x relative 
to the other ozone catalysts, implying a smaller 
potential effect of NO x from HSCT exhaust on 
stratospheric ozone [e.g., Weisenstein etal., 1991]. 
All the 2-D predictive models include heterogeneous 
chemistry on sulfate aerosols in the stratosphere 
using specified sulfate aerosol distributions. For the 
assessment calculations, PSCs were not included in 
the 2-D model runs. With the reduction in the 
importance of NO x as a catalyst of ozone depletion, 
HO x from exhaust water vapor emissions takes on a 
larger role. Thus, the 1993 AES A Interim 
Assessment Report [NASA, 1993] and this 
assessment emphasize model evaluations of the 
combined effects of NO x and H 2 O from HSCT 
exhaust. The current report, however, also discusses 
other factors that previously were given less 
emphasis, such as PSC formation (see Section 
3.2.4.6) and the effect of sulfur and soot emissions 
on aerosol loading in the stratosphere (see Sections 
3.2.4.5 and 3.2.4.7). 

2.3 Model Intercomparisons 

Variations among the predictive model results arise 
from combinations of differences in the model com- 
ponents of transport, radiation, and photochemistry. 
Two exercises were carried out recently to examine 
the differences in radiation and photochemistry. 

2.3.1 PHOTOLYSIS BENCHMARK 

Radiation subroutines within 2-D models usually 
make approximations to speed up the calculation. 
During the Models and Measurements exercise 
[Prather and Remsberg, 1993] and in previous model 
intercomparisons, differences were discovered 
among model photolysis calculations even when the 
cross sections and overhead ozone were specified. 
For this assessment, a benchmark photolysis 
calculation has been developed for which several 
detailed radiation models [Anderson and Lloyd, 
1990; Dahlback and Stamnes, 1991; M. Prather, 
personal communication] have achieved the same 
answers to within 2%. The benchmark gives us the 
ability to check directly the precision of the radiation 
modules used in the 2-D predictive codes. The 
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Figure 18. Photolysis benchmark intercomparison for J(N02) as a function of altitude for 13 models for a 
clear-sky calculation at O' solar zenith angle and ground albedo of 0.3. The solid line is the UCI benchmark. 
The models with similar values (>0.015 sec* 1 ) in the upper troposphere are multi-stream calculations, while 
the lower group are two-stream calculations. The AER model used in the runs now includes the UCI code. 


initial intercomparisons yielded differences between 
some of the 2-D model radiation routines and the 
benchmark. These variations were deemed unac- 
ceptable by the 2-D modeling groups, who then up- 
graded their radiation calculations or converted to a 
parameterized look-up table derived from a more ac- 
curate radiation code. One example of the results of 
this intercomparison, the photolysis rate for NO 2 as a 
function of altitude, is shown in Figure 18. The 
spread of the calculations for this benchmark case is 
~±5% in the upper stratosphere. In the lower strato- 
sphere and troposphere, the results fall into two 
groups: multi-stream calculations with photolysis 
coefficients peaking at a little over 0.015 sec- 1 , and 
two-stream calculations with peak photolysis coeffi- 
cients under 0.015 sec- 1 . This is significantly better 
agreement than in previous intercomparisons. 

2.3.2 Chemistry Solver Benchmark 

After considering chemical kinetic theory and 
practice, numerical accuracy, and computational 
speed, a variety of photochemistry solution 


techniques have been chosen by the 2-D predictive 
model groups. In the same manner as the photolysis 
benchmark discussed above, a chemistry benchmark 
calculation was used to identify differences in the 
steady-state solutions calculated by the different 
chemistry solvers. 

This initial benchmark was defined using a Gear- 
type integration code with a consensus set of 
chemical species and reactions common to the 
majority of the 2-D predictive models using rate 
coefficients from JPL-94-26 [DeMore et al., 1994]. 
The chemical code in each 2-D model was used to 
calculate the constrained-equilibrium concentrations 
of the radical species with fixed values for ozone; 
CH 4 ; H 2 O; CO; molecular hydrogen (H 2 ); and fixed 
family totals of inorganic NO y , Cly, and Br y . In the 
calculation, all models used the same photolysis 
rates. 

Three cases or regions were chosen to represent 
situations in which different families dominate ozone 
loss: A) 40 km/2.7 hPa, 65°N, April 11, with fixed 
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species and families measured or inferred from the 
Atmospheric Trace Molecule Spectroscopy 
(ATMOS) measurements outside the polar vortex 
between April 8 and 14, 1993; B) 20 km/67 hPa, 
38°N, May 1 1 with concentrations from the SPADE 
aircraft campaign diurnal flights in 1993; and C) 20 
km/64 hPa, 59°N, November 4 with concentrations 
from the ASHOE/MAESA aircraft mission northern 
survey flight in 1994. 

The benchmark diumally averaged volume mixing 
ratios for several important radical species were 
compared to 2-D model results. Models (e.g., 
CSIRO and the Lawrence Livermore National 
Laboratory/ University of Illinois (LLNL/UI)), which 
calculate diurnal coefficients directly using short 
time steps, agree with a variety of key radicals to 
within better than 1%. The Atmospheric and 
Environmental Research, Inc. (AER) model, which 
also calculates diurnal coefficients directly but with a 
longer time step to speed up the calculations, agrees 
with the benchmark to within a few percent for key 
radicals. The Goddard Space Flight Center (GSFC) 
model gets somewhat larger differences, from a few 
percent up to 10 to 15%. This arises from their more 
approximate approach to calculating diurnal effects 
so as to obtain sufficient model speed to make many 
runs. 

The benchmark mechanism results have also been 
compared for region C to the solver module in the 
photostationary state (PSS) model used by Salawitch 
and coworkers [e.g., Salawitch et al., 1994a, b] in 
their studies of aircraft observations. (The photolysis 
rates used in the benchmark were provided by 
Salawitch.) Diurnal-average radical concentrations 
are almost uniformly consistent to within 1% 
between the benchmark result and PSS solver results 
with a closely similar mechanism. This begins to 
provide a connection between the predictive 2-D 
models and the data now available to evaluate their 
accuracy. 

2.4 Model Predictions of Ozone Response 
to HSCT Emissions 

2.4.1 The Scenarios 

A scenario is defined by specifying the boundary 
conditions for the trace gases in the atmosphere and 
the emissions from a HSCT fleet. To obtain the 
ozone change for each scenario, the 2-D model 


background case is run using the prescribed 
boundary conditions without any HSCT emissions 
(but with projected subsonic emissions for the year 
2015) until a repeating seasonal cycle is achieved. 
The perturbed case includes emissions due to the 
HSCT fleet and subsonic emissions modified to 
account for subsonic flights displaced by the HSCT 
fleet. The model calculated change in ozone for the 
scenario is given by the difference in the calculated 
ozone distributions between the background and the 
perturbed cases. 

Aircraft emissions of NO x , H2O, HCs (as CH4), and 
CO were included in the calculations. Some of the 
assessment scenarios were developed to evaluate the 
sensitivity of the effect of HSCTs to changes in the 
background atmosphere (e.g., background chlorine 
amount), and changes in the HSCT fleet (e.g., 
changes in altitude of injection, El no*, and fleet 
size). In most scenarios, the trace gas concentrations 
in the background atmosphere are those predicted to 
be appropriate for 2015. In particular, the ground - 
level concentrations of CFCs were specified to give a 
stratospheric mixing ratio for Cl y = 3 parts per billion 
by volume (ppbv), consistent with the projected 
concentration given in WMO [1995]. The bromine 
concentration is assumed to be at 19 parts per trillion 
by volume (pptv) [WMO, 1995]. Additional 
calculations were performed for Cl y = 2.0 ppbv 
corresponding to the chlorine level for 2050, and 0.6 
ppbv corresponding to the natural background. 
Variations in the HSCT fleet include different values 
of El NO* (0. 5, 10 and 15), different cruise altitudes, 
and different fleet size. 

Background aerosol surface areas, assuming 
relatively little volcanic activity, were as given in the 
lower limit case specified in WMO [1992] and were 
assumed not to be perturbed by the aircraft operation. 
The sensitivity of the HSCT ozone perturbations to 
ambient sulfate aerosol loading is discussed in 
Section 3. 2. 4.4. Effects from emitted sulfur 

compounds (due to fuel impurities) are discussed in 
Section 3.2.4.5. 

For the calculations, all models used the chemical re- 
action rate coefficients and photolysis cross sections 
recommended in the most recent Jet Propulsion 
Laboratory (JPL) review [DeMore et al . , 1994]. The 
runs did not include the heterogeneous reaction of 
B1ONO2 + H2O or any consideration of PSCs. 
Calculations on the effects of B1ONO2 hydrolysis 
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Table 4a. Calculated steady-state total column ozone change (%) in the Northern Hemisphere averaged over 
a year for the basic assessment scenarios. 


Scenario 

Model 

■ 

Mach Number 

HS9| 

EInOx 

(g N0 2 /kg fuel) 

Number of 
Aircraft 

AER 

GSFC 

CSIRO 

LLNLV 

UI 

CAMED 

■ 

2.4 

3 

0 

500 

-0.3 


-0.2 

-0.1 

— 

n. 

2.4 

3 

5 

500 

-0.3 


-0.2 

-0.3 

+0.08 

m. 

2.4 

3 

10 

500 

-0.4 


- 0.3 

-0.5 

-0.4 

IV. 

2.4 

3 

15 

500 

-0.7 


-0.5 

-0.9 

-1.0 

■a 

2.4 Cruise +1 km 

3 

5 

500 

-0.4 


-0.3 

-0.5 

+0.02 

VL 

2.4 Cruise +1 km 

3 

15 

500 

-1.2 

m 

-0.9 

-1.4 

-1.4 

vn. 

2.4 Cruise -2 km 

3 

5 

500 

-0.04 


+0.0 

-0.06 

+0.1 

vra. 

2.4 Cruise -2 km 

3 

15 

500 

-0.03 


+0.08 

-0.1 

-0.3 

DC. 

’”**2.4 

2 

5 

500 

-0.4 

El 

- 0.3 

-0.4 

+0.02 

X. 

2.4 

2 

15 

500 

-1.0 

Ksl 

-0.7 

-1.2 

-1.2 

XI. 

2.4 

0.6 

15 

500 

-1.5 

bub 

— 

— 

— 

xn. 

"" 2.4 

3 

5 

1000 

-0.6 

ppj 

-0.4 

-0.7 

+0.03 

xm. 

2.4 

3 

15 

1000 

-1.9 

E9 

-1.4 

-2.1 

EH 


and PSCs are discussed in Sections 3. 2. 3. 2 and 
3.2.4.6, respectively. 

2.4.2 CALCULATED CHANGES IN COLUMN 

Amount of Ozone 

Analysis of the results from these updated models 
confirm the previous finding that the effect on ozone 
from the emitted CO and HCs (as CH 4 ) is small. 
Sensitivity calculations from the National Center for 
Atmospheric Research (NCAR) model [Tie et al., 
1994] show that the effect in the model of the 
emitted HC, when treated as nonmethane 
hydrocarbons (NMHCs), is also small. 

Table 4a gives model-calculated steady-state ozone 
changes for a number of basic assessment scenarios 
relative to the all-subsonic aircraft cases, averaged 
over a year for the Northern Hemisphere, while 
Table 4b gives calculated changes for a latitude band 
between 40°N and 50°N. The basic HSCT fleet was 
the Mach 2.4, 500-aircraft fleet with Els of 0, 5, 10, 
and 15 g NO x /kg fuel burned (scenarios I through 
IV). The EIno x = 0 case represents the effects from 


HSCT H 2 O emissions in the absence of interference 
from emitted NO x . For the EInox = 0. 5, and 10 
cases, all of the models calculated ozone depletions 
in the Northern Hemisphere of <0.5%. The 
calculated ozone depletions generally increase with 
larger EIno x - All the models, except GSFC, 
calculated a Northern Hemisphere total ozone 
decrease of between -0.5 and -1% for EIno x = 15. 
Figure 19 shows the calculated ozone column change 
in the 40 to 50° latitude band as a function of EInox- 
Note that while most models show a similar response 
to changes in EIno x , the gradient calculated by the 
University of Cambridge/University of Edinburgh 
(CAMED) model is much steeper. 

The effect of changing flight altitudes is shown in 
Figure 20. Scenarios V and VI are cases where the 
Mach 2.4 fleet was assumed to be cruising 1 km 
higher in altitude than in the basic scenarios. The re- 
sults for EInox = 5 were more negative than the 
basic scenario. For EInOx = 15, most models 
calculated depletions of greater than 1%. The results 
shown for scenarios VII and VIII are for calculations 
in which the Mach 2.4 fleet was assumed to cruise 2 
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Table 4b. Calculated steady-state total column ozone change (%) between 40°N and 50°N averaged over a 
year for the basic assessment scenarios. 



Scenar 

io 

Model 

■ 

Mach Number 

Cly 

(ppbv) 

ElNOx 

(g NC>2/kg fuel) 

Number of 
Aircraft 

AER 

GSFC 

CSIRO 

LLNL/ 

UI 

CAMED 

I. 

2.4 

3 

0 

500 

-o^r 

-0.3 

-0.3 

-0.1 

— 

n. 

2.4 

3 

5 

500 

-0.3 

-0.05 

-0.2 

-0.3 

+0.2 

m. 

2.4 

3 

10 

500 

-0.5 

+0.07 

-0.3 

-0.6 

-0.4 

IV. 

2.4 

3 

15 

500 

-0.8 

+0.06 

-0.5 

-1.0 

-1.2 

m 

2.4 Cruise +1 km 

3 

5 

500 

™-oJ" 

-0.2 

-0.4 

-0.6 

+0.09 

VI. 

2.4 Cruise +1 km 

3 

15 

500 

-1.3 

-0.2 

-1.1 

-1.6 

-1.5 

VII. 

2.4 Cruise -2 km 

3 

5 

500 

-0.03 

+0.07 

+0.04 

-0.06 

+0.2 

VIII. 

2.4 Cruise -2 km 

3 

15 

500 

-0.01 

+0.3 

+0.2 

-0.2 

-0.2 

IX. 

2.4 

2 

5 

500 

-0.4 

-0.2 

-0.3 

-0.4 

+0.07 

X. 

2.4 

2 

15 

500 

-1.1 

-0.3 

-0.8 

-1.4 

-1.4 

XI. 

2.4 

0.6 

15 

500 

-1.7 

-0.7 

— 

— 

— 

XII. 

2.4 

3 

5 

1000 

-0.7 

-0.2 

-0.5 

-0.7 

+0.1 

XIII. 

2.4 

3 

15 

1000 

-2.2 

-0.6 

-1.6 

-2.3 

— 



El (NOx) 

Figure 19. Predicted column ozone change (%) in the 40° to 50°N latitude band as a function of EInOv for a 
fleet of 500 Mach 2.4 HSCTs. 
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Models: 



AER 

GSFC 

CSIRO 

LLNL 

CAMED 


500 HSCTs 
EI(NOx) =5 


Altitude shift relative to Mach 2.4 (km) 


Figure 20. Predicted column ozone change (%) in the 40' to 50’N latitude band as a function of altitude 
displacement from the Mach 2.4 base case, assuming 500 HSCTs with EInq x = 5. 


km lower in altitude than the basic scenario, 
simulating a Mach 2.0 fleet. In these cases, all 
models calculated very small (<0.3%) depletions for 
both EInOx = 5 ar *d 15. This demonstrates the strong 
sensitivity of the effects of NO x to altitude. All the 
models show a similar dependence on flight altitude. 

The effect of decreasing chlorine was tested by run- 
ning the Mach 2.4 fleet with EInOx = 5 and 15 and a 
background chlorine amount of 2 ppbv (expected to 
be reached in -2050), and the EIno x = 15 case with 
Cly = 0.6 ppbv (scenarios IX through XI). As 
expected, the models calculate a more negative 
perturbation as the chlorine background decreases 
because of reduced interference between the C10 x 
and NO x cycles. The calculated globally and 
annually averaged ozone de-pletion varies linearly 
with the Cly concentration and is more negative for 
smaller Cly. 

Scenarios XII and XIII are for the double Boeing 
fleet described in Section 2.1. The increased 
perturbation over the basic Mach 2.4 fleet 
calculations for EInOx = 5 is nearly linear, except for 
the GSFC and CAMED models which started with a 
result near zero. The perturbation increase is greater 
than linear for EInox = 15. Note that the calculated 


depletion for the double fleet at EIno x = 5 is greater 
than for the EIno x = 10 case with the smaller fleet. 
This is primarily because the mature fleet is not just a 
doubling of the flights in the same corridors, but 
includes a higher percentage of flights crossing the 
tropics where the exhaust is more readily lofted to 
high altitudes. Another difference is that the double 
fleet at El no,, = 5 emits twice the water vapor as a 
fleet of 500 HSCTs with EIno x = 10. 

Figure 21 shows the calculated depletion in column 
ozone as a function of latitude for each model for the 
month of March. All of the models show a greater 
calculated ozone change near the poles than at low 
latitudes. There is significant variation in the latitude 
dependence because the various models have a 
different crossover point in altitude from ozone 
increase in the lower stratosphere to ozone decrease 
in the middle and upper stratosphere (see Section 
2.4.3). At 60°N, the model-calculated ozone changes 
for the EInOx = 15 case range from near zero for the 
GSFC model to -1.3% for the LLNL/UI model. For 
EInOx = 5, the range is from -0.1% to nearly -0 .5%. 
The model-to-model variation is less in the Southern 
Hemisphere. 
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EL,,* = 5 gN0 2 /kg fuel EI N0 , = 15 gN0 2 /kg fuel 



Figure 21. Calculated column ozone change (%) during March as a function of latitude for a Mach 2.4 HSCT 
fleet from the several models; a) EIno x = 5 g N02/kg fuel, and b) EIno x = 15 g N 02 /kg fuel. The changes 
reflect the difference between a calculation of the projected 2015 supersonic and subsonic fleets and the 
projected 2015 subsonic fleet in the absence of a supersonic fleet. 


2.4.3 Calculated local changes in ozone 
Concentration 

The above discussions represent a brief overview of 
the column integrals of ozone change. The change in 
ozone column is made up of a calculated decrease in 
the middle to upper stratosphere and a calculated in- 
crease in the lower stratosphere and troposphere. 
The balance between these increases and decreases is 
dependent on the dispersion of the aircraft effluent 
by atmospheric transport processes as well as the 
transport of long-lived species, including ozone, in 
the background atmosphere. An example of the 
dispersion of aircraft-generated NO x (converted to 
NO y ), calculated by the AER model, is shown in 
Figure 22. All of the models used in this study 
obtained similar total amounts of perturbed NO y , but 
showed differences in its distribution with altitude 
and latitude. The calculated added amount of NO y in 
the HSCT-perturbed atmosphere corresponded to a 
residence time of 13 to 16 months. Figure 23 shows 
the corresponding calculation of ozone change as a 
function of latitude and season from the AER model 


for EIjmo x = 5 and 15 g NO 2 /kg fuel. The residence 
time for exhaust gases and the photochemical 
response time of ozone in the lower stratosphere 
imply that the perturbation would decrease to less 
than half its value in a few years if HSCT operations 
in the stratosphere were to cease. 

Figure 24 compares the calculated change in ozone 
concentration as a function of altitude from 5 models 
at 45°N for the Mach 2.4 El = 5 case. All of the 
models predict ozone depletion at altitudes above 24 
km, while most of the models predict ozone increase 
between 14 and 20 km. With the exception of the 
CAMED model, the models predict relatively small 
contributions to the integrated ozone column due to 
their treatment of tropospheric processes. To isolate 
the changes that are the result of changes in the sub- 
sonic fleet, an extra model scenario was run in which 
the emissions above 1 3 km altitude were assumed to 
be zero. The models obtained ozone changes of 0 to 
-0.03% relative to the subsonic-only case. This 
demonstrates that the tropospheric emissions from 
supersonic aircraft and the change in subsonic emis- 
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AER - NOy June (El=5 - Subsonic) % Diff. 




AER - NOy June (El=1 5 - Subsonic) % Diff. 



Figure 22. Calculated change in the total odd 
nitrogen (NO») concentration during June (in 
percent) as a function of latitude and altitude for a 
Mach 2.4 HSCT fleet from the AER model; a) EIno x 
= 5 g N 02 /kg fuel, and b) EIno x = 15 9 N 02 /kg fuel. 
The changes reflect the difference between a 
calculation of the projected 2015 supersonic and 
subsonic fleets and the projected 2015 subsonic 
fleet in the absence of a supersonic fleet. 

sions in the HSCT scenario are small effects on 
column ozone. 

Ozone profiles as a function of altitudes calculated 
by the models are very similar at altitudes above ~26 
km (see Figure 25). Between approximately 16 and 
26 km, the calculated ozone profiles show a wider 
spread of results. This is a critical region for 
understanding the HSCT perturbation and is the 
altitude band where the variations between models of 
HSCT impact at 45 °N latitude were greatest. These 
differences in ozone concentrations in the lower 
stratosphere are, at least in part, due to transport 
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Figure 23. Calculated change in the ozone 
concentration during June (in percent) as a function 
of latitude and altitude for a Mach 2.4 HSCT fleet 
from the AER model; a) EIno x = 59 N 02 /kg fuel, 
and b) EIno x = 15 g N 02 /kg fuel. The changes 
reflect the difference between a calculation of the 
projected 2015 supersonic and subsonic fleets and 
the projected 2015 subsonic fleet in the absence of a 
supersonic fleet. 


effects. Partitioning of the chemical families is a 
function of the ozone concentration (see Section 
3.2.3.2). The differing background ozone amounts in 
the models will thus lead to different chemical 
sensitivities to added NO x . 

Figure 26 shows the percentage of ozone loss in the 
background atmosphere due to each catalytic cycle 
as a function of altitude and latitude calculated in a 
representative 2-D model (LLNL/UI) for March. 
Catalytic destruction due to NO x is one of the 
significant background atmospheric ozone loss 
mechanisms in the middle and upper stratosphere. In 
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Figure 24. Calculated change in ozone 
concentration as a function of altitude for a fleet of 
500 Mach 2.4 HSCTs (EIno x = 5 )- The results are 
shown at 45’N latitude for the month of March. 



Ozone (10 12 mol cm* 3 ) 


Figure 25. Calculated ozone concentrations as a 
function of altitude for the 2015 all-subsonic scenario 
(Cly = 3 ppbv). The results are shown at 45'N 
latitude for the month of March. 
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Figure 26. Percentage of ozone loss rate due to each catalytic cycle as a function of latitude and altitude for 
March in the 2015 atmosphere without HSCTs from the LLNL model: NO x , HO x , CIO x , and BrO x . BrO x 
chemistry does not include hydrolysis of BrONC >2 on aerosols. Reactions involving bromine and chlorine 
molecules are classified as BrO x . Reactions involving chlorine and nitrogen are classified as CIO x . Reactions 
involving nitrogen and hydrogen are classified as NO x . 
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the lower stratosphere, NO x catalysis is a minor 
mechanism and loss is dominated by HO x catalysis 
with significant contributions from C10 x and BrO x . 
Similar calculations for the other models yield 
qualitatively similar results. They do, however, 
differ in detail. Some of this can be attributed to 
variations in the dynamics that yield differing 
background concentrations, while part can be 
attributed to variations in the counting schemes used 
to assign O x loss to each chemical family. When 
rate-determining steps for catalytic cycles are used, 
several differing assignments for O x loss can be 
made which yield equivalent results. 

One analysis that can be made is the change in cat- 
alytic cycles as NO y and H 2 O are added from the air- 
craft exhaust. Figure 27 shows the CSIRO model- 
calculated change in the loss term from each catalytic 
cycle when the aircraft perturbation is added for 
EInOx = 15 during March. The quantity plotted is 


-100 • (LNOx(pert) - Lno x (base))/L to tal (base) for the 
NO x term in the upper left panel and the equivalent 
terms for the HO x , C10 x , and BrO x cycles in the 
other panels. In the region below 20 km, the added 
NO x and H 2 O each interfere with both the HO x and 
the NO x cycles, causing a partial cancellation of the 
negative effects of the NO x cycle. Nearer to the 
pole, the BrO x cycle interferes with the negative 
NO x perturbation. 

The model results show that there is no simple 
spatial correspondence between changes in NO y and 
changes in production and loss. The portion of the 
added NO y that is transported to the tropics above 
-25 to 30 km altitude leads to a large change in 
ozone removal rate. This is because the sulfate 
aerosol surface area decreases at higher altitudes, and 
photolysis of HNO 3 is faster in the tropics and at 
high altitudes, leading to an increase in the NO x /NO y 
ratio and, therefore, an increase in the sensitivity to 
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Figure 27. Change in loss term for each catalytic cycle divided by total chemical loss term (in percent) for 
March from CSIRO model. Changes are calculated for EIno x = 15 9 N02/kg fuel: change in NO x loss term, 
change in HO x loss term, change in CIO x loss term, and change in BrO x loss term (not including Br0N02 
hydrolysis). Reactions involving bromine and chlorine molecules are classified as BrO x . Reactions involving 
chlorine and nitrogen are classified as CIO x . Reactions involving nitrogen and hydrogen are classified as 
NO x . 
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NO x perturbations. The NO x catalysis reactions 
proceed in this region with less interference from the 
HO x , C10 x , and BrO x reactions. The locations of 
large ozone depletions are not necessarily coincident 
with locations of large changes in the local 
photochemical terms. In the lower mid-latitude 
stratosphere, the models generally calculate a 
decrease in ozone (in percent) in spite of a small de- 
crease in the local ozone removal rate. The ozone 
depletion that the models calculate in the lower mid- 
latitude stratosphere is the result of chemical losses 
which occurred at lower latitudes and higher 
altitudes and transport. The ozone-depleted air from 
the tropical middle stratosphere has then been 
transported downward and poleward, yielding the 
calculated ozone decreases in the area of the flight 
corridors. 

2.5 Model Predictions of the Effects of 
HSCT Emissions on Climate Forcing 

Potential HSCT climate effects result from the direct 
injection of CO 2 , H 2 O, soot, and sulfur. Indirect ef- 
fects on climate are derived from the exhaust effects 
on ozone, cirrus formation, and background sulfate 
aerosol concentrations. The direct effect of CO 2 
from HSCTs is expected to be small relative to that 
from fossil fuels because, at the present time, the 
projected HSCT fleet emissions of CO 2 are less than 
-0.3% of the emissions from all fossil fuel 
combustion. 

Regional climate changes that may follow from a in- 
creased climate forcing cannot be predicted at this 
time. A common approach to estimating the relative 
importance of various atmospheric changes to 
climate is to determine the difference in stratospheric 
radiative forcing at the tropopause. Lacis et al. 
[1990] showed that the radiative forcing due to 
changes in ozone is sensitive to the vertical 
distribution of the projected ozone changes, with a 
maximum forcing sensitivity per unit ozone 
occurring in the lower stratosphere and upper 
troposphere. Thus, the calculated changes in ozone 
shown in Section 2.4 have potential for affecting 
climate. 

Rind and Lonergan [1995] used the low-resolution 
Goddard Institute for Space Studies (GISS) middle- 
atmosphere 3-D general circulation model (GCM) to 
predict the climate effects of the projected ozone 


changes calculated by the GSFC 2-D model. The 
GSFC model predicted ozone increases of up to 1 % 
between 1 5 and 30 km in the tropics and decreases 
between 30 and 50 km, with a maximum decrease of 
2.5% varying from 25 km at the poles to 35 km in 
the tropics for a Mach 2.4 fleet with EIno x = 15 
(similar to that shown in Figure 23). This profile 
change leads to a net positive forcing and, in the 
GISS model, produced a predicted global-average 
surface-temperature change of +0.025 °C. Because 
the standard deviation of annual average 
temperatures in the GISS control run was 0.25°C to 
0.34°C, depending on the length of integration, this 
predicted change is not statistically significant. The 
predicted surface temperature change was small, 
partly because the effect of the tropospheric ozone 
increase offset the effect of the stratospheric ozone 
decrease. 

However, at high latitudes the projected 2015 ozone 
change that was used had smaller tropospheric in- 
creases. Thus the thermal forcing was more consis- 
tently negative and, combined with amplified climate 
feedbacks due to sea ice and dynamical changes, a 
relatively consistent reduction of a few tenths of a 
degree occurred for the surface air temperature. This 
value is still smaller than the natural interannual vari- 
ability in either the model or surface observations; its 
significance will depend on model sensitivity and the 
ozone forcing used. 

In the stratosphere, the ozone reductions predicted by 
most models, on the order of a few percent, resulted 
in in situ temperature reductions of a few tenths of a 
degree, while dynamical changes produced 
additional effects of up to ±1°C. The dynamical 
changes resulted from a decrease in vertical stability 
as the stratosphere cooled. Once again, these 
changes are smaller than the model or observed 
interannual variability [Rind and Lonergan, 1995; 
Shah et al., 1995], although they were fairly 
consistent. 

Another calculation by Rind and Lonergan [1995] 
was an idealized perturbation of H 2 O (i.e., an 
increase of 0.2 parts per million by volume (ppmv) 
everywhere in the stratosphere). This calculation 
produced results that were inconsistent during a 50- 
year simulation. In one five-year average a warming 
of 0.06°C was found, while in a second a cooling of 
-0.07°C was found. Within the statistical 
fluctuations of their model, the resultant surface 
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temperature perturbation is not significant. They 
also tested the model sensitivity to changes in 
stratospheric H 2 O by assuming a uniform increase 
from 3 to 6 ppmv. This calculation resulted in a 
consistent surface temperature warming of 0.24°C. It 
is not clear at this time whether longer simulation 
times would yield consistent perturbation predictions 
for the smaller perturbation in water vapor. 

In the stratosphere, the increase in water vapor 
generally produced cooling of a few tenths of a 
degree, although dynamical changes associated with 
the decreased vertical stability and an increase in 
residual circulation of <5% produced some regions 
of greater effects. These changes are smaller than 
model or observed interannual variability [Rind and 
Lonergan, 1995; Shah et al., 1995], but once again 
were somewhat persistent. 

One interesting aspect of these perturbations, as 
highlighted by Rind and Lonergan [1995], is the fact 
that the climate sensitivity due to changes in the 
altitude distribution of radiatively important species 
can be different from the climate sensitivity derived 
from CO 2 experiments. The climate sensitivity is 
defined as the change in global-average surface 
temperature for a given radiative forcing at the 
tropopause. They found a reduced sensitivity 
relative to that calculated for a CO 2 forcing in both 
the ozone and H 2 O perturbations. They traced the 
result to the altered response of high clouds in their 
model to the changes in dynamical stability in the 
case of ozone, and to the radiative forcing gradient in 
the case of H 2 O perturbation. High clouds increase 
in CO 2 perturbations, thereby increasing their 
contribution to the greenhouse effect. In both the 
ozone and H 2 O perturbations, high clouds either 
decreased (H 2 O) or increased less than anticipated 
for a similar forcing by CO 2 (ozone), creating a 
smaller sensitivity to the same forcing. This result 
should be viewed with caution, because different 
models will respond differently to the various 
radiative perturbations (e.g., for typical GCMs, 
climate sensitivities due to a doubling of CO 2 vary 
between 1 .5 and 4.5°C). However, it introduces the 
possibility that results may depend on the specific 
scenario introduced and that the results of different 
types of radiative perturbations may not be simply 
additive [compare Taylor and Penner, 1994]. Based 
on the results of Taylor and Penner [1994], it also 
may be important to sort out the effects of 


hemispheric differences in radiative forcing due to 
aircraft emissions. 

The direct radiative and climatic effects of increased 
concentrations of sulfate aerosols and soot have not 
been assessed. However, the simulations of Weaver 
et al. [1995] for NO y can be used to estimate their 
possible radiative effects. Weaver et al. [1995] used 
the 3-D GSFC chemical transport model (CTM) to 
simulate the distribution of NO y for a Mach 2.4 
scenario with El of 5 g N02/kg fuel. The 
stratospheric residence time for both sulfate aerosols 
and soot should be similar to that calculated for NO y , 
although small differences will occur because of the 
fall velocities of the aerosols. The concentration of 
sulfate or soot aerosol is then given by the ratio of 
Els scaled by the appropriate molecular weights. For 
Els of 0.4 g SO 2 /kg fuel and 0.02 g soot/kg fuel, this 
scaling leads to maximum vertical column 
concentrations of 3 x 10- 4 g m-2 for sulfate and 10-5 
g m-2 for soot extending over a broad latitude belt 
from 40°N to 75°N in January. The numbers 
obtained by this simple scaling for sulfate are a 
factor of ~4 larger than those calculated by 
Weisenstein et al. [1995] in their 2-D model. 
Assuming specific extinction cross sections for 
visible light of 10 m2 g-1 gives total optical depths 
for sulfate and soot of only 0.003 and 0.0001, 
respectively. These perturbations to aerosol optical 
depth are less than one-tenth and one one-hundredth, 
respectively, of the global-average perturbations due 
to tropospheric anthropogenic sulfate aerosols, and 
are much smaller than the peak optical depths for the 
Mt. Pinatubo volcanic aerosol which was +0.3. 
Thus, the radiative perturbations from sulfate and 
soot emissions are expected to be small if these 
estimates are accurate. Additional work is needed to 
define more precisely the expected sulfur and soot 
concentration changes and their calculated climatic 
impact, as well as their potential to act as cloud CN 
in the lower stratosphere/upper troposphere. 

The final aircraft climate issue concerns the effect of 
H 2 O emissions on possible contrail formation in the 
stratosphere. 2-D models predict an increase of ap- 
proximately 0.8 ppmv H 2 O at latitudes of maximum 
flight frequencies. Such loading might increase 
stratospheric water vapor to ~4-6 ppmv near 53-90 
mbar. Using average temperatures of -56 to -58°C at 
these heights, saturation mixing ratios with respect to 
ice are 100 to 400 ppmv. Thus, in an average case 
the dry stratosphere would ensure that widespread 
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formation of contrails would not occur as a result of 
HSCT emissions. However, there remain questions 
about contrail formation, especially in winter when 
much lower temperatures can occur. 

In summary, the climate changes which have been 
calculated from HSCT ozone and water vapor 
emissions tend to be small and less than model or 
observed interannual variability, but the in situ 
stratospheric responses tend to be consistent during 
the course of the simulation. Other effects (e.g., 
those due to sulfur and soot) are also likely to be 
small, but need to be quantified and tested. 

The temperature changes shown above for aircraft 
effects on ozone and water vapor are equilibrium 
simulations, using estimated 2015 releases. Due to 
the heat capacity of the ocean, this equilibrium 
response would not be felt for some 30 years (by 
which time the aircraft perturbations would likely 
have changed). To provide a complete perspective 
of these results, aircraft-induced climate changes 
should be compared with those anticipated to occur 


by 2015 from other perturbations, particularly 
anthropogenic-induced trace-gas effects (primarily 
CO 2 ). According to IPCC [1995], with a doubled 
CO 2 sensitivity of some 1.5° to 4.5°C by 2015, 
surface temperature warming of a few tenths of a 
degree or more is likely and stratospheric cooling of 
similar or somewhat greater magnitude is possible 
[Rind et al . , 1990]. Hence the aircraft results which 
would have occurred by 2015 would likely be 
considerably smaller in comparison. 

Note that neither the direct temperature effects of al- 
tered ozone nor the dynamically induced temperature 
changes were incorporated in assessment model pho- 
tochemical calculations; neither were the model- 
calculated residual circulation changes. In addition, 
prospective changes in stratospheric temperatures or 
dynamics due to other anthropogenic perturbations 
(e.g., CO 2 ) by 2015 were not included. Sensitivity 
tests of the impact of these in situ changes should be 
made to incorporate in an iterative fashion potential 
climatic impacts and photochemical predictions. 



3 ASSESSMENT UNCERTAINTIES 

2-D model predictions for ozone changes in response 
to HSCT fleet operation are given in Chapter 2. 
Confidence in the construction quality of the models 
has increased since the NASA 1993 Interim 
Assessment Report [NASA 1993] by applying com- 
putational tests or “benchmarks” to the models, as 
described in Section 2.3. However, the degree to 
which the 2-D models adequately represent all of the 
qualitative and quantitative features of the atmo- 
sphere continues to be limited. This chapter defines 
the major areas of uncertainty associated with the 
quantitative assessment of HSCT impacts. These 
uncertainties are quantified where possible, either 
through formal propagation of errors in the model 
calculations or through model sensitivity studies. On 
account of the nonlinear nature of the atmospheric 
system, many of the uncertainties to be discussed in 
this chapter cannot be quantified at this time. 
Strategies to quantify those areas of uncertainty in 
the future are identified. 

For the purposes of this chapter, the areas of uncer- 
tainty are organized around two major topics, the air- 
craft emissions which serve as input to the predictive 
models (Section 3.1) and the model representations 
of atmospheric processes (Section 3.2). The mod- 
eled processes include photochemistry, transport, and 
aerosol microphysics. The treatment of these pro- 
cesses in the models included in Chapter 2 ranges 
from very detailed (photochemistry) to highly pa- 
rameterized (microphysics). A summary of the ma- 
jor uncertainties and their relationship to the overall 
confidence in the results of Chapter 2 is discussed in 
Section 3.3. 

3.1 Aircraft Emissions 

Future HSCT emissions cannot be forecast with great 
accuracy; hence, they contribute to the uncertainty of 
the overall assessment. Model results presented in 
Chapter 2 include inventories for aircraft-generated 
NO x , H 2 O, CO, and HCs (as CH 4 ). Measurable 
quantities of soot and SO x are also emitted by air- 
craft but were not treated in those calculations. Input 
of the inventories into the models is achieved by 
simple dilution of the aircraft plume at the altitude of 
injection, with no chemical changes taking place, 
into a 1 ° longitude by 1 ° latitude by 1 km altitude 
grid. Further scaling of the exhaust composition is 
required when models employ larger grid sizes (e.g., 


zonal-average 2-D models), again with the same 
“frozen” chemistry assumption. Since HSCT en- 
gines have not been built, there are no direct obser- 
vational constraints that can be placed on the esti- 
mated Els of the important exhaust species, on 
potential chemical changes occurring on the model 
subgrid scale (i.e., in the aircraft plume/wake), or on 
the nature of the exhaust dilution process from 
subgrid to grid scale (e.g., degree of anisotropy, 
subsidence, etc.). An assessment of the model input 
quality requires that the uncertainties associated with 
each of these estimates or assumptions relating to 
aircraft exhaust and plume evolution be evaluated. 

Uncertainties in the HSCT emission inventories de- 
rive from several sources. Inventory estimates for 
any of the exhaust species are affected by uncertain- 
ties in the total amounts of HSCT fuel burn and its 
geographic partitioning. As stated in Section 2.1, the 
flight routing assumed in this assessment represents 
an upper limit to the amount of emissions derived 
from a fleet of the specified size. As HSCT technol- 
ogy matures, better estimates of fuel consumption 
will be possible. 

Emission inventory uncertainties associated with the 
individual exhaust species vary widely. NO x emis- 
sions are subject to engineering tradeoffs in combus- 
tor design. Since initial assessments of HSCTs iden- 
tified NO x as having the largest potential impact on 
calculated ozone changes, the NO x emission inven- 
tory has been treated parametrically in this assess- 
ment, with an El range of a factor of 3 chosen to en- 
compass both the HSCT design goal and the present 
state-of-the-art subsonic engine performance. 
Emission inventories for CO 2 and H 2 O are directly 
related to the stoichiometric conversion of the fuel 
amounts. Because combustor efficiencies are greater 
than 99.9%, the uncertainty in CO 2 and H 2 O due to 
this factor is < 0 . 1 %, far less than that associated with 
the total fuel consumption estimate used for the 
HSCT fleet. Thus, the uncertainty in CO 2 and H 2 O 
emissions is approximately equal to the uncertainty 
in fuel consumption. 

The effects of CO 2 emissions from an HSCT fleet 
are not considered here. Although emitted primarily 
into the stratosphere, CO 2 will become well mixed 
throughout the atmosphere, contributing to the over- 
all greenhouse gas forcing of the climate. In contrast 
to NO x , the CO emissions represent only a small 
(-10%) perturbation to ambient levels. Moreover, 
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the effect of additional CO on ozone levels is indi- 
rect, occurring only through repartitioning of HO x . 
The degree of repartitioning resulting from increased 
CO is small relative to that arising from the NO x 
perturbation. Hence, the uncertainty in CO emis- 
sions has a negligible effect on uncertainty in ozone 
response. The same reasoning applies to HC exhaust 
in the form of CH4. The potential increase in CH4 
levels is essentially negligible relative to ambient 
levels, even assuming stoichiometric conversion of 
the HC fraction to CH4. 

Chemical composition studies of commercial and 
military subsonic engine exhaust [Spicer etal., 1992, 
1994] find that a significant fraction (>50%) of the 
HC exhaust at cruise conditions exists in forms other 
than CH4. The NMHC fraction, if partitioned be- 
tween a small number of species, might constitute 
significant perturbations to the ambient levels of the 
produced species. However, the experimental tests 
indicate that a myriad of NMHC species are pro- 
duced in highly variable amounts, depending on en- 
gine type. Propagation of the uncertainties in 
NMHC emissions to calculated ozone response is 
difficult, owing to photochemical reactivity differ- 
ences between the plethora of NMHC species. 
Based on the current understanding of HC chemistry, 
all NMHCs serve to produce ozone, to sequester 
chlorine into HC1, and to sequester nitrogen in the 
form of nitrites and nitrates. The latter pathway will 
reduce the aircraft-induced change in NO x levels; 
however, the effect is negligible since aircraft NO x 
emissions are, at a minimum, 50 times greater than 
the NMHC emissions. The effect of the first two 
pathways is more difficult to gauge since it involves 
a comparison of cycle efficiencies. A model sensi- 
tivity study to aircraft ethane emissions [Tie, 1994] 
indicates that changes in ozone are less than 0.1% for 
values of ethane slightly higher than levels consid- 
ered possible by this assessment. The predicted 
change in ozone due to ethane emission increases 
under volcanically perturbed aerosol conditions to 
around 0.4%. 

Soot and sulfur emissions from aircraft are, along 
with water vapor, the primary ingredients in particle 
formation; however, they are not included in the cal- 
culations described in Chapter 2. The impact of this 
omission on the assessment uncertainty is discussed 
in Section 3.3. 


Omission of NO x emission plume effects from the 
2-D model input is supported by plume/wake model 
predictions of species in the near-field, which show 
no significant impact in that region. The near-field 
plume and wake has been modeled with 
NO x /SO x /HO x chemistry coupled to fluid mechani- 
cal calculations of the exhaust mixing and dynamics 
[Miake-Lye etal., 1993, 1994; Quackenbush etal., 
1993]. Such calculations indicate that relatively 
small amounts of NO y deposited at the point of wake 
break-up are in the form of HONO or HNO3, the 
bulk remaining in the form of NO x as emitted. 
Measurements in the wake of the ER-2 [Fahey et al., 
1995b], the Concorde [Fahey et al. , 1995a], and sev- 
eral commercial subsonic airliners [Arnold et al. , 
1992] provide field experimental support for the 
modest conversion of emitted NO x to other NO y 
species. 

Another omission from 2-D models is the descent of 
exhaust captured in the aircraft vortex wake, both 
through downward transport due to vortex motion 
and through radiative cooling due primarily to in- 
frared emission of the exhaust water vapor. 
However, the degree of descent is expected to be less 
than a few hundred meters due to vortex descent and 
less than 1 km even for the maximum cooling rates 
occurring at mid- to high latitudes in the winter 
hemisphere [Rodriguez et al., 1994]. In situ obser- 
vations of the Concorde [Barrilleaux, 1995] and ER- 
2 [Fahey etal., 1995b] aircraft exhaust are consistent 
with a small (<1 km) plume subsidence. As de- 
scribed in Section 2.4.2, lowering the injection alti- 
tude by 2 km results in calculated ozone changes that 
are -0.2% and -0.5% less negative for fleet Els of 5 
and 15, respectively. Since the subsidence is likely 
much less than 2 km, an upper limit on the uncer- 
tainty in predicted column ozone changes due to un- 
certainty in the injection altitude is, by linear ex- 
trapolation of the model results, reasonably set at 
0 . 1 %. 

3.2 Model Representation of the 
Atmosphere 

A primary goal of the global stratospheric models 
used in this assessment is to make accurate simula- 
tions of ozone, for both the historical, observed at- 
mosphere and for future atmospheres. This goal can 
be achieved only by accurately representing the 
physical processes controlling atmospheric ozone in 
the models. The ability of the models to meet this 
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goal can be evaluated by comparison with existing 
and future observations and by attempts to reproduce 
measurements of past atmospheric changes. 
However, the evaluation of the 2-D predictive mod- 
els using observations is not straightforward. A ma- 
jor difficulty is that a number of mechanisms in the 
model combine to determine ozone at any particular 
model location. Therefore, the possibility exists that 
the models predict ozone well, even though the rela- 
tive contributions to the ozone value from the various 
model mechanisms are in error. In such a situation, 
the response of a model to a perturbation could differ 
significantly from the atmospheric response. A sec- 
ond difficulty in evaluation is that models are for an 
average, (i.e., typical) atmosphere, while measure- 
ments are taken under specific and perhaps different 
atmospheric conditions. The reasonable performance 
of a model in comparison with a number of different 
observations, each designed to test a specific model 
mechanism, will produce confidence that the models 
will be better suited to predict aircraft effects. 

A systematic evaluation of 2-D models with mea- 
surements in connection with the AESA program has 
been carried out [Prather and Remsberg, 1993]. The 
models have since been revised and improved; cur- 
rent versions and their more complete variants (e.g., 
those including PSCs) need to be reexamined. There 
is now a much more extensive set of aircraft obser- 
vations from the SPADE and ASHOE/MAESA cam- 
paigns, and a set of satellite observations from in- 
struments on the Upper Atmosphere Research 
Satellite (UARS) and the Atmospheric Laboratory 
for Application and Science (ATLAS) payload on 
the Space Shuttle. This chapter will discuss some of 
the valuable comparisons that have been made in the 
past and describe some new comparisons for evalua- 
tion of model response to HSCT perturbations. An 
evaluation of the overall capabilities of the models 
by comparison with ozone observations follows. 
Subsequent discussion will treat the specific model 
components of transport, photochemistry, and parti- 
cle formation. 

3.2.1 Model Comparisons with Ozone 
Observations 

3.2.1. 1 Ozone Climatology 

The Total Ozone Mapping Spectrometer (TOMS) 
provides a global picture of column ozone distribu- 
tion. Figure 28, reproduced from WMO [1992], 


compares TOMS observations with several of the 
2-D models included in this assessment report. 
Although the models have evolved since this com- 
parison was made, it provides a useful illustration 
here. TOMS total ozone reaches a maximum in 
spring in both the Northern and Southern 
Hemispheres. The maximum is on the pole in the 
Northern Hemisphere and at ~60°S in the Southern 
Hemisphere. The magnitude of the Southern 
Hemisphere maximum is smaller than the Northern 
Hemisphere value by some 100 Dobson units (DU). 
The timing of the column ozone maxima varies by 
about a month amongst the models. The models lo- 
cate the maxima at high latitudes in both hemi- 
spheres but do not consistently place the Northern 
Hemisphere maximum at the pole and the Southern 
Hemisphere maximum off the pole. Summertime 
mid- to high- latitude values of column ozone are 
relatively consistent among the models and agree 
fairly well with the TOMS observations. The models 
consistently underestimate the mid-latitude seasonal 
variability seen in TOMS column ozone. Equatorial 
values of column ozone in the models are typically 
within 10% of the TOMS values. However, the lati- 
tudinal extent of the low tropical values in the mod- 
els is typically narrower than indicated by TOMS, 
leading to relatively weak latitudinal gradients in 
model column ozone compared to TOMS observa- 
tions. 

The general behavior of the calculated column ozone 
as a function of latitude suggests that the gross trans- 
port characteristics of the models are correct. The 
differences between the models and the observations 
suggest that there are significant transport uncertain- 
ties in the lower stratosphere. It should be noted that 
the usefulness of column ozone as a simple diagnos- 
tic of transport processes is somewhat limited be- 
cause photochemistry still plays a significant role in 
determining column ozone. 

Comparisons of observed and modeled profile ozone 
test different model features at various latitudes and 
altitudes. In the upper stratosphere, the comparison 
is a test of the overall effect of model photochem- 
istry. In the lower stratosphere, ozone is sensitive to 
lower stratospheric transport but it is difficult to de- 
termine the specific cause of any disagreement. 
Comparisons of model profiles with Solar 
Backscatter Ultraviolet Spectrometer (SBUV) data 
indicate that models produce the correct shape of the 
mixing ratio profile in equatorial regions and at high 
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Figure 28 . Column ozone abundances in Dobson units circa 1980 from TOMS observations as a function of 
latitude and season. Also shown are the results from three models used during the 1991 UNEP/WMO 
Scientific Assessment of Ozone Depletion; these models are not identical to the current versions of the same 
model. [WMO, 1992] 


latitudes in the summer. The ozone peak in the 
models is a few kilometers too low in comparison to 
either SBUV or Stratospheric Aerosol and Gas 
Experiment (SAGE II) profiles. At higher latitudes 
in other seasons, the altitude of the modeled ozone 
peak tends to be ~5 km lower than observations sug- 
gest. Finally, models consistently underpredict 
ozone in the 35 to 50 km regions and overpredict 
ozone at low altitudes, especially in the equatorial 
regions. Recent work shows that the upper strato- 
spheric ozone production and loss rates balance 
when measured values of the ozone and radical con- 
centrations are used [Minschwaner et al., 1993; 
Crutzen et al., 1995; Dessler et al., 1995a]. 

3.2.1.2 Ozone Trends: 1980-1990 

Atmospheric ozone has been observed globally since 
the late 1970s. During the same time, stratospheric 
chlorine levels have been increasing and much of the 
change in ozone is thought to result from these chlo- 


rine increases. The observed ozone response to a 
well-defined perturbation can be compared to the 
calculated ozone response to test the predictive capa- 
bility of the model, provided that increases in halo- 
gen levels are the sole cause. Although the HSCT 
fleet represents a NO x and H 2 O perturbation, the 
close connection between catalytic ozone loss cycles 
makes it difficult to argue that a model which fails to 
predict the atmospheric response to a chlorine per- 
turbation can correctly predict the response to a NO x 
perturbation. 

Trends in column ozone have been estimated from 
TOMS observations and compared to model calcula- 
tions [Stolarski et al., 1991; Bojkov et al., 1990, 
1994; Randel and Cobb, 1994; Reinsel et al., 1994]. 
Gas-phase models substantially underpredict the ob- 
served trends, while models including heterogeneous 
reactions on the sulfate aerosols do much better 
[Rodriguez et al., 1991; Brasseur and Granier, 1992], 
Figure 29, from Chapter 6 of WMO [1995] compares 
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Figure 29. Model-calculated changes in column 
ozone between 1980 and 1990 compared to trends 
derived from observations. [WMO, 1995] 


TOMS, SBUV, and Dobson station total ozone 
trends as a function of latitude with a number of 2-D 
models. The models are in reasonable agreement 
with the measurements in the tropics and lower mid- 
latitudes. At higher latitudes, the models tend to un- 
derestimate the observed trends significantly. In 
particular, the large Southern Hemisphere trend ob- 
served by the satellite measurements throughout the 
year is not captured by the models. Some of this dis- 
crepancy is due to the neglect of PSC processes in 
some of the models and a relatively poor model rep- 
resentation of the polar winter vortex in both hemi- 
spheres. There are also large discrepancies between 
the models and the observations in the Northern 
Hemisphere mid-latitudes, particularly in the late 
winter and spring. It is not known whether the dis- 
crepancy at these latitudes is due to neglected in situ 
chemistry or incorrect representations of lower 
stratospheric transport processes. 

The above model calculations neglected the 
variability of the sulfate aerosol layer over the past 
decade, using a satellite-derived climatology instead. 
A recent study by Solomon et al., [1995] finds that, 
if the temporal variations in the sulfate aerosol layer 
are included, then many of the features of the 
observed trend can be reproduced. However, the 
magnitude of the changes in this study was still 
underestimated by a factor of ~1.5 in comparison to 
TOMS trends. This could be due to an incorrect 
model formulation. The trends calculated from 


TOMS also have some uncertainty. For example, 
trends calculated from Dobson measurements are not 
as strongly negative as those from TOMS [WMO, 
1992], 

The models also can be evaluated by comparison 
with observed profile ozone trends such as those 
shown in Figure 3. This figure illustrates the March 
trends in profile ozone as a function of latitude and 
altitude for 1980 to 1990. Models predict peak frac- 
tional decreases to occur in the 40- to 45-km regions, 
with values similar to the satellite-derived trends. At 
-25 km, SAGE II observations indicate a broad re- 
gion of statistically insignificant trends extending 
from 50°S to 50°N. In the models, the width of this 
region generally extends from 30°S to 30^, with 
trends increasing more rapidly toward higher lati- 
tudes in the models than in the SAGE II derived 
trends. The models generally fail to reproduce the 
very large decreases in low- to mid-latitude lower 
stratospheric ozone derived from the SAGE data. 
Trends in the mid-latitude lower stratosphere derived 
from ozonesonde measurements are less than half of 
the SAGE values at similar latitudes [Logan, 1994; 
Miller et al, 1995; WMO, 1995]. However, the 
model-calculated trends at 15 km in the mid-latitudes 
are generally about half again as small as the sonde- 
derived trends. This underestimate of the lower 
stratospheric mid-latitude trends by the models is 
consistent with the results obtained from comparison 
to column ozone trends mentioned above. Thus, the 
overall pattern and location of ozone change are 
predicted, but the remaining discrepancies in the 
lower stratosphere are troubling since this region is 
critical to the HSCT problem. The reasons for the 
discrepancy between the models and the observations 
are not known. 

3.2.2 ATMOSPHERIC TRANSPORT 

The net impact of aircraft emissions is related to how 
fast the mid-latitude lower-stratosphere emissions are 
dispersed to other regions of the atmosphere. Much 
of the current emphasis is on how material crosses 
from mid-latitudes into the tropics where the pollu- 
tant can then be efficiently transported to the middle 
and upper stratosphere. Just as important, however, 
is transport from the lower stratosphere into the tro- 
posphere, which effectively cleanses the strato- 
sphere. In addition, there is the issue of how aircraft 
emissions are transported into the wintertime polar 
vortex. Finally, there is the possibility of aircraft 
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emissions remaining in the corridors during times of 
weak winds, especially during the summer. 

Transport of stratospheric air between mid-latitudes 
and either the tropics or the polar vortex is driven by 
episodic “wave” events. Similarly, the actual trans- 
port of lower stratospheric emissions into the tropo- 
sphere is linked to large-scale “synoptic” wave activ- 
ity. All of these wave processes are intrinsically 3-D 
and must be parameterized in the 2-D models. The 

2- D wave parameterizations have focused primarily 
on stratospheric planetary waves. 

3- D models explicitly resolve planetary and synoptic 
waves, and are adequately advanced to provide a 
useful tool for evaluating the transport characteristics 
of 2-D models. Generally, 3-D models require a pa- 
rameterization of gravity waves on a smaller spatial 
scale than the model grids, and hence do not provide 
a direct evaluation of gravity-wave transport. 
Currently the most useful integrations are for a small 
number of reactive species or, more commonly, 
long-lived tracers with parameterized loss terms. 
The usual approach to the 3-D chemistry transport 
application is to first determine the 3-D wind and 
temperature fields and then to use these off-line to 
calculate the transport and chemical production and 
loss of constituents. 

There are two approaches to obtaining 3-D wind 
fields. One is to use a self-deterministic GCM. The 
resulting winds are internally consistent solutions to 
the specified equations of motion. Even if the model 
were perfect, the chaotic nature of the weather would 
lead GCM integrations to depart from the behavior of 
the atmosphere and, at best, calculate a climatologi- 
cal-average state. Nevertheless, these models have 
been able to describe many of the averaged observa- 
tions of long-lived chemicals in the stratosphere such 
as N 2 O, ozone, and CO 2 [see e.g., Hall and Prather, 
1993, 1995; Strahan et al., 1994], In addition to 
these chemical tracer tests, one needs to evaluate the 
GCM winds and temperatures with climatological 
averages and statistics. One problem with many 
GCMs is that the computed climatological-average 
state contains known but uncorrectable biases such 
as excessively cold temperatures at the winter pole 
that cause systematically biased winds. 
Unfortunately, chemical studies using GCMs are es- 
pecially sensitive to the impact of temperature biases 
on aerosol processes. 


The second approach is to use winds derived from a 
data assimilation system. In data assimilation, a 
GCM is run and observations are objectively melded 
into the model during the integration [Daley, 1991]. 
Data assimilation techniques have been used widely 
in numerical weather forecasts and are being actively 
extended to more general applications [e.g., Schubert 
et al., 1993]. Through data assimilation vast 
amounts of routine meteorological observations are 
brought to bear on the transport simulation problem, 
and there is systematic comparison of winds and 
temperatures from the assimilation with radiosonde 
and aircraft observations. Data assimilation provides 
high-quality estimates of temperatures, extratropical 
horizontal winds, and other derived quantities that 
can be used to reduce the transport-related uncer- 
tainties in aircraft assessments made using off-line 
CTM calculations. While this approach cannot be 
used to evaluate feedbacks between changes in 
chemistry and climate for an assessment, when com- 
pared with self-deterministic GCMs, its strengths 
include relatively unbiased estimates of present-day 
temperature (no “cold pole”), proper representation 
of planetary- and synoptic-scale variability, proper 
location of climatological features (e.g., vortex, 
Aleutian anticyclone, storm tracks, etc.), timely rep- 
resentation of seasonal transitions, and representation 
of some important characteristics of the quasi-bien- 
nial oscillation (QBO). There are also weaknesses in 
this approach. The residual circulation is especially 
difficult to calculate because it is the difference of 
two large competing processes. In addition, the 
residual circulation is influenced by the data inser- 
tion process, which adds non-negligible forcing 
terms to the momentum and thermodynamic equa- 
tions which impact the residual circulation. Winds in 
the tropics and subtropics are also highly uncertain 
because the temperature observations that dominate 
the input data for the assimilation do not effectively 
constrain the wind. Nevertheless, in 3-D models us- 
ing assimilation winds agree well with ozone obser- 
vations when integrated over yearly time scales 
[Douglass et al., 1995]. Although the impact of er- 
rors in the residual circulation becomes more impor- 
tant with longer integrations, 3-D model integrations 
provide constraining information on all of the trans- 
port uncertainties discussed below. Such constraints 
have not been systematically incorporated into the 
2-D models used in this assessment. 


40 



3.2.2.1 Model Comparisons with Tracer 
Observations 

When results from the 1992 Models and 
Measurements Workshop were summarized [Prather 
and Remsberg, 1993], there were limited observa- 
tions of stratospheric tracers capable of providing 
constraints on how well global-scale transport is rep- 
resented in the 2-D models. Comparison of tracer 
distributions showed a range of agreement between 
observations and the various models. One consistent 
and significant area of disagreement was the tropical 
lower stratosphere. Simultaneous measurements of 
NOy and ozone spanning the tropics and mid-lati- 
tudes in the lower stratosphere were available from 
the tropical 1987 Stratosphere-Troposphere 
Exchange Project (STEP) and Airborne Antarctic 
Ozone Experiment (AAOE) missions [Murphy et al., 

1993] . The difference in the NO y /ozone ratio be- 
tween the tropics and mid-latitudes was simulated by 
the 2-D models [Prather and Remsberg, 1993] and 
provides an indicator of the transport of mid-latitude 
air into the tropics both in the real atmosphere and in 
the models. The models tended to smooth out the 
tropical/extratropical difference in this ratio, suggest- 
ing that the subtropical transport barrier was not well 
represented in the models. Recently, a simulation of 
the effects of a stronger subtropical barrier on the 
NOy/ozone ratio was calculated using a reduced 
horizontal eddy coefficient (K yy ) in the subtropics in 
the AER model [Donnelly et al., 1995]; this steep- 
ened the gradient in NOy/ozone across the subtropics 
and more closely simulated the observations. 

Measurements in the 1960s and 1970s of carbon- 14 
( I 4 C) and strontium-90 that resulted from atmo- 
spheric weapons testing during the 1950s and 1960s 
[summarized by Johnston, 1989; Kinnison et al., 

1994] presented one of the few opportunities to eval- 
uate transport in the models independent of chem- 
istry. Due to the limited data set, however, definitive 
tests of model transport schemes cannot be made 
with 14 C. The comparisons that have been made 
serve to highlight potential model deficiencies. For 
instance, there have been several studies of bomb 
14 C [e.g., Shia et al., 1989, 1993; Jackman et al., 
1991; Kinnison et al., 1994; Rasch et al., 1994]. The 
3-D calculation by Rasch shows the best agreement 
with observations. In particular, there is less pollu- 
tant in the upper stratosphere than in the 2-D calcu- 
lations, which better agrees with observations. The 
physical reasons for this difference have not been 


determined. Although Plumb and Mahlman [1987] 
found good agreement between simulations with a 
GCM and a 2-D model derived from the same GCM, 
this approach has not yet been utilized in 14 C simu- 
lations. The studies of Jackman et al. [1991] and 
Rasch et al. [1994] seem to indicate that the 2-D 
models were, if anything, too slow in their overall 
transport when compared to 3-D models. On the 
other hand, a comparison by Kinnison et al. [1994] 
of 2-D calculations and observed 14 C years after the 
last major bomb tests indicated that the models 
tended to get a profile with altitude that was constant 
above -25-30 km, while the data suggested a strong 
peak and falloff. 

More recent data from the 1992 airborne SPADE 
mission has also been used to characterize the nature 
of atmospheric transport. Using a 1-D global dif- 
fuser model and observed CO 2 and N 2 O correlations 
from the SPADE mission, Wofsy et al. [1994] de- 
duced a more rapid transport of pollutants to the 
middle and upper stratosphere than is found in cur- 
rent 2-D models. Much of the character of these ob- 
servations of CO 2 and N 2 O were matched in a 3-D 
model [Hall and Prather, 1995], and these tests need 
to be examined with the current 2-D predictive 
models used here. On the other hand, as noted 
above, the 2-D models appear to exaggerate the up- 
ward transport of bomb 14 C. These comparisons, 
along with the NOy/ozone ratios discussed above, 
suggest that the 2-D models used in this assessment 
are not global diffuser models. They may underes- 
timate but do not neglect tropical isolation. The 
global diffuser model does not provide a good repre- 
sentation of stratospheric transport. 

In the last several years, new observations of a vari- 
ety of long-lived species have been obtained by in- 
struments on aircraft, satellites, and balloons which 
hold the promise of providing important new tests of 
the representation of transport in models used for 
predicting HSCT effects. The global and seasonal 
coverage of CH 4 , N 2 O, CFOs, hydrofluoric acid 
(HF), and H 2 O measurements from UARS will allow 
seasonal changes in stratospheric circulation to be 
deduced and will provide constraints for models. 
Correlations between the tracers N 2 O, CFC-1 1, NO y , 
ozone, CH 4 , and H 2 O measured by ATMOS and the 
ER-2 aircraft now span the tropics and extratropics 
and will allow a more thorough examination of 
model transport between the tropics and mid-lati- 
tudes. Initial analyses of the data suggest that trans- 
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port in the real atmosphere falls somewhere between 
the two conceptual models of the global diffuser and 
the tropical pipe and may be a strong function of alti- 
tude. Both the in situ and satellite data are consistent 
with estimates that approximately 20 to 45% of the 
air at 20 km in the tropics (equator ±7°) was trans- 
ported from mid-latitudes while air above 20 km in 
the tropics is significantly more isolated (see Section 
1.2.2) [Avallone and Prather, 1995; Minschwaner et 
al., 1995]. However, a more thorough analysis of the 
data and systematic comparisons with model predic- 
tions are required before the transport barriers can be 
completely characterized. In addition, in situ obser- 
vations of seasonally varying correlations between 
CO 2 , H 2 O, and N 2 O will be needed to quantify verti- 
cal and quasi-horizontal transport rates in the lower 
stratosphere. 
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3.2.2.2 Dynamical Barriers: Tropical/Mid- 

Latitude Transport and Polar Vortices 

The global diffuser and the tropical pipe are the ex- 
treme descriptions of mid-latitude transport into the 
tropics, and hence, transport of pollutants to the 
middle and upper stratosphere. Calculations with the 
AER model were made to test the sensitivity of the 
HSCT ozone perturbation to the inclusion of a pipe. 
They simulated the effects of a pipe on transport by 
reducing the K yy in the subtropics at altitudes be- 
tween 14 and 24 km to restrict transport between 
mid-latitudes and the tropics. The addition of the 
pipe resulted in a more positive calculated ozone re- 
sponse to the HSCT fleet, in the range of a few tenths 
of a percent to nearly 0.5%. Figure 30 shows the 
difference between the two model simulations. Note 
that while the tropics and Southern Hemisphere are 
distinctly different, Northern Hemisphere HSCT im- 
pacts are similar in the two cases. 

As noted in Section 3.2.2. 1, observations indicate 
that air in the tropics at 20 km consists of 25 to 40% 
air characteristic of mid-latitudes. Existing 2-D pre- 
dictive models have yet to be diagnosed for this 
quantity, but their tropical NO y /ozone ratio suggest 
that they exaggerate the exchange between mid-lati- 
tudes and the tropics. Agreement with observations 
improves in the AER model when a subtropical bar- 
rier is introduced. 

Simulations using 3-D models also address the ques- 
tion of the subtropical barrier. Weaver et al. [1995] 
carried out a five-year calculation of the spread of 



Month 

Figure 30. Calculation of the column ozone change 
as a function of month and latitude for the 201 5 
Mach 2.4 HSCT fleet with EIno x = 1 5 g N 02 /kg fuel 
from the AER model; a) with "normal" dynamics, 
i.e., no reduction in eddy coefficients in subtropics, 
and b) with reduction of subtropical eddy coefficients 
by a factor of 10 to simulate effects of a subtropical 
barrier (or tropical pipe). 


NO y from a fleet of Mach 2.4 HSCTs using their 3-D 
CTM with winds from the GSFC assimilation. After 
five years, the zonal-average distribution of the NO y 
perturbation compares reasonably with that calcu- 
lated using the GSFC 2-D model (see Figure 31). 
However, the 3-D model does confine the aircraft 
pollutant to mid-latitudes more effectively than the 
2-D model. In addition, as diagnosed with satellite 
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Figure 31. Comparison of 2-D and 3-D calculations 
for June of the added NO y (in ppbv) from a fleet of 
Mach 2.4 HSCTs as a function of latitude and 
altitude (pressure). Contoured values have been 
scaled to an EIno x = 15 g N 02 /kg fuel for 
comparison with other published figures. Dashed 
lines are for the 2-D calculation and solid lines are 
for the 3-D calculation. [Weaver et al., 1995] 

data, the winds used in the 3-D model overmix the 
subtropics with mid-latitudes in the lower strato- 
sphere [Douglass etal., 1995]. This adds more evi- 
dence that the current models incorporate more mid- 
latitude air into the tropics than actually occurs. 

Most 2-D models used in the assessments also do not 
represent the observed isolation of the polar vortices, 
though more dynamically based transport parameter - 
izations can simulate the isolation [Garcia et al., 
1992]. The implications of this model deficiency for 
the HSCT assessment have not been addressed. 

3.2.2.3 Stratosphere-Troposphere Exchange and 
Residence Time of Emissions 

The discussion above focused on the issue of how 
much pollutant gets into the upper stratosphere 
where ozone is more vulnerable to NO x . Another 
important transport issue is stratosphere-troposphere 
exchange. The subject of stratosphere-troposphere 
exchange has a vast literature, as recently reviewed 
by Holton et al. [1995]. Movement of air into the 
troposphere is the primary mechanism for eliminat- 
ing HSCT exhaust products from the stratosphere. 

The simulation of stratosphere-troposphere exchange 
is discussed by Shia et al. [1993]. In 2-D predictive 
models, this exchange is accomplished through ad- 
vection by the residual circulation and diffusive mix- 


ing across the tropopause. It has a tendency to occur 
at high latitudes rather than mid-latitudes, as found in 
3-D models [e.g., Douglass et al., 1993] and analyses 
of atmospheric data. Additionally, tropopause fold- 
ing in some 2-D models is represented as diffusion 
that takes place along isentropes that are in the 
stratosphere at mid-latitudes and in the troposphere 
in the tropics. These models also misplace the 
stratosphere-troposphere exchange when compared 
with 3-D models, as indicated by large amounts of 
stratospheric material being transported into the 
tropical upper troposphere. This transport is dis- 
cussed by Shia et al. [1993]. 

An important aspect of the differences in strato- 
sphere-troposphere exchange for understanding 
HSCT effects is in the overall residence time. 
Jackman etal. [1991] did a sensitivity study of three 
different versions of the residual stream function. 
Their results confirmed that the build-up of pollu- 
tants in the lower stratosphere was directly related to 
the strength of the residual circulation. Comparisons 
of calculations with observations for total ozone and 
tracers N 2 O and CH 4 can be used to determine 
whether a particular circulation is possible. The 
strongest circulation used by Jackman et al. [1991] 
was shown to be impossible within the framework of 
the model, i.e., constituent distributions calculated 
using the model with this circulation were very dif- 
ferent from observations. The circulations that were 
considered possible each produced better qualitative 
comparisons with some observations. There is a 
significant difference in the maximum build-up of 
pollutants allowed using these circulations. Results 
show a significant sensitivity to the circulation, even 
when limited to “possible” by comparisons with 
long-lived species. This sensitivity is damped signif- 
icantly when heterogeneous chemistry is included in 
the calculations. 

In the Weaver et al. [1995] 3-D/2-D intercomparison 
discussed in 3.2.2.2, the calculated change of NO y in 
the 3-D model is in general only slightly less than 
that in the 2-D model. Since the injection rates are 
the same, this indicates that the model residence 
times of the injected NO y are almost equal (~13 
months). Thus, while the stratosphere-troposphere 
exchange mechanisms in 2-D models are very differ- 
ent, the integrated effect on residence time appears to 
be similar. It is not clear how to reconcile this with 
the 2-D and 3-D model differences found in the 14 C 
transport studies described above. 
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3.2.2.4 Longitudinal Asymmetries and Corridor 
Effects 

In addition to assumptions implicit in their transport 
parameterizations, 2-D models must assume zonally 
uniform sources whereas, in reality, emissions will 
be concentrated in localized corridors. This raises 
two issues of transport: how should these localized 
sources be mapped into 2-D, and what is the degree 
and significance of emission concentrations within 
the corridors? 

The sensitivity of the assessment to specification of 
the 2-D source of emissions was addressed by Plumb 
et al. [1995]. In the first method, the 2-D source was 
specified as a simple zonal average of the 3-D source 
distribution (as is done in the assessments). In the 
second method, the potential vorticity/potential tem- 
perature mapping technique of Schoeberl et al. 
[1989] was used to map the 3-D source into pseudo- 
Lagrangian coordinates. The resulting zonally aver- 
aged source then included the effect of short-term 
transport processes that moved emissions away from 
initial locations. Only very small differences were 
found in predicted ozone loss using these 2-D source 
distributions, indicating that this issue is not one of 
significant uncertainty in the assessments. 

The second 3-D issue concerns the possibility that 
some regions in the atmosphere will accumulate 
significantly more exhaust than would be expected 
using a zonal-mean model. In the study of Weaver et 
al. [1995], the largest zonal asymmetry in the pertur- 
bation is found in summer, when light winds fail to 
ventilate the corridors; this is in accord with results 
reported by Sparling et at. [1995], In late summer 
during one of the extreme events, the maximum 
mixing ratio on the 54-mbar surface was only a fac- 
tor of ~2 larger than the zonal mean at the same lati- 
tude. This suggests that the zonal averaging required 
for a 2-D calculation does not produce a large sys- 
tematic error in the assessment calculations. 

3.2.3 PHOTOCHEMISTRY 

Knowledge of the chemical reaction rate coefficients 
has improved significantly in the last 20 years. 
Critical evaluations of the laboratory measurements 
of rate coefficients and estimates of uncertainty are 
made every two years by a NASA panel [see e.g., 
DeMore et al., 1994]. Major changes in the last two 
UNEP/WMO assessments occurred with the recog- 


nition of the importance of heterogeneous chemistry. 
The estimates of uncertainty by the NASA panel can 
be propagated through the 2-D model formulation. 
An example is described in Section 3.2.3.4. 

Uncertainty in the gas-phase, heterogeneous, and 
photolytic processes can also be constrained by ob- 
servations. Measurements of free radicals and 
longer-lived reservoir species made from aircraft, 
satellite, and balloon-borne platforms provide a test 
of the ability of photochemical process models to ac- 
curately partition the chemical families between 
reservoir and radicals. Provided that the photochem- 
istry modules in the 2-D predictive models contain 
the full chemistry explored with the process models, 
knowledge gained in the observations/model inter- 
comparison is transferred into the 2-D predictive 
models. However, such an evaluation does not di- 
rectly test the distribution of free radicals in the 2-D 
predictive models. Since the concentrations of free 
radicals will also depend on the transport of the long- 
lived reservoirs, the distribution of ozone, and the 
temperature. 

3.2.3.1 Middle and Upper Stratosphere 
Photochemistry 

Observations from balloon, shuttle, and satellite- 
borne instrumentation have improved the 
understanding of processes that regulate the 
concentration of hydrogen, nitrogen, and chlorine 
radicals in the middle and upper stratosphere. 
Concentrations of hydrogen radicals in this region 
are primarily sensitive to the abundance of ozone, 
excited atomic oxygen (0( ! D)), H2O, CH4, HNO3, 
and NO. Profiles of OH and HO2 as a function of 
solar illumination have been measured by the 
Smithsonian Astrophysical Observatory Far Infrared 
Spectrometer (FIRS-2), Far Infrared Limb Observing 
Spectrometer (FILOS), and Stratospheric Limb 
Sounder (SLS) balloon-borne remote sensing 
instruments. These observations show good 
agreement (±30%) with theoretical profiles 
constrained by simultaneous measurements of the 
major precursors, as illustrated in Figure 32, for the 
range of altitudes (25 to 40 km) over which balloon- 
borne observations are available [Chance et al., 
1995; Pickett and Peterson, 1995]. A single in situ 
observation of OH and HO2 between -25 and 37 km 
shows broad agreement (±100%) with theory 
[Stimpfle et al . , 1990; Wennberg et al . , 1990], These 
comparisons demonstrate that the primary sources 
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Figure 32. Measured and calculated OH and HO2 
at midday. Measurements were obtained by FIRS-2 
near 34’N, on 26 and 27 September 1989. Error 
bars represent estimates of the 1-a total accuracy of 
measurement, including uncertainties associated 
with molecular parameters, residuals from fits to the 
measured spectra, and intensity calibration of the 
spectra. Theoretical profiles appropriate for the local 
measurement time were found using a 
photochemical model [Salawitch et al., 1994b] 
constrained by profiles of ozone, H2O, and 
temperature measured by FIRS-2; profiles of NOu, 
Cl y> and CH4 estimated from FIRS-2 measurements 
of N2O; and sulfate aerosol surface area measured 
by SAGE II. [Adapted from Chance etal., 1995] 


and sinks of OH and HO 2 should be well represented 
by 2-D predictive models using recommended rates 
and cross sections [DeMore et al., 1994], provided 
fields of precursors (ozone, H 2 O, and CH 4 ) are cal- 
culated in a realistic manner. However, existing ob- 
servations of hydrogen radicals in the middle and 
upper stratosphere are limited to mid-latitudes near 
equinox above altitudes of -25 km. Observation of 
OH and HO 2 in the altitude range 20 to 25 km, for 
latitudes and seasons not yet sampled, are essential 
for testing 2-D predictive models. 

Simultaneous measurements of radicals (NO and/or 
NO 2 ) and major reservoirs (HNO 3 , or total NO y ), 
and ozone (which mediates the conversion efficiency 
of NO x to HNO 3 ) guide our understanding of 
processes that regulate the abundance of NO x gases. 
Concentration profiles of NO y gases NO, NO 2 , 
HNO 3 , CIONO 2 , N 2 O 5 , HNO 4 , and ozone at sunrise 
and sunset have been obtained by ATMOS at a vari- 
ety of latitudes and seasons throughout the middle 
stratosphere. Theoretical profiles of the partitioning 


of NOy gases, constrained by simultaneous mea- 
surements of ozone, H2O, and CH4 from ATMOS 
and the surface area of sulfate aerosol from SAGE II, 
show good agreement with observations, as illus- 
trated in Figure 9. These observations demonstrate 
the NO/NO2 ratio is represented well by models. 
Furthermore, agreement between theory and obser- 
vation for N 2 O 5 suggests that the mechanism for 
conversion of NO x to HNO3 is well simulated by 
photochemical models. Similar levels of agreement 
between theory and observation of NO y gases for 
altitudes between 25 and 40 km are demonstrated by 
balloon-borne measurements of NO, NO2, HNO3, 
and CINO3 by the JPL MkIV interferometer [Sen et 
al., 1995] and measurements of NO2, HNO3, and 
ozone by the Balloon-borne Laser In Situ Sensor 
(BLISS) [Webster etal., 1994]. These observations 
demonstrate that 2-D predictive models should calcu- 
late realistic fields of NO and NO2 for aerosol con- 
ditions appropriate for near-background (i.e., non- 
volcanic) conditions, provided precursor fields (e.g., 
ozone, NOy, Cly, H2O, and CH4) are calculated in a 
realistic manner. 

There are limited measurements of NO y gases in the 
middle stratosphere under conditions highly per- 
turbed by volcanic aerosol. Measurements by 
ATMOS and the JPL MkIV interferometer are lim- 
ited by the inability of their sun trackers to remain 
locked on the solar disk. Measurements of NO 2 by 
SAGE II are compromised by the high extinction as- 
sociated with the volcanic aerosol. Ground-based 
measurements, however, reveal large decreases in 
column NO 2 and increases in column HNO 3 due to 
the Mt. Pinatubo aerosol over Lauder, NZ (47°S) 
[Koike et al., 1994]. The magnitude of these 
changes is underestimated by nearly a factor of 2 by 
2-D models [Koike etal., 1994], suggesting possible 
missing sinks for NO x at high surface area. 
However, balloon-borne in situ measurements within 
volcanic clouds of profiles of NO, NO y , ozone, and 
surface area over Aire sur 1' Adour, France (44°N, 
O’W) on 20 October and 6 November 1992 [Kondo 
etal., 1995], and of NO, NO y , CIO, and ozone over 
Greenland ( 66 °N) during March 1992 [Dessler et al., 
1993] reveal good agreement between observed and 
calculated NO x . These observations suggest the 
major sink for NO x in volcanic clouds is the hetero- 
geneous hydrolysis of N 2 O 5 , and they appear to rule 
out a rapid rate for other sinks at temperatures higher 
than 200 K. Additional measurements of NO y gases 
for volcanically perturbed conditions would be desir- 
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Figure 33. Partitioning of inorganic chlorine species 
in the stratosphere. The triangles represent the 
measured sum of HCI sunset + CINO 3 sunset + CIO 
midmorning; lines labeled "Cly" represent the 
abundance of total inorganic chlorine necessary to 
match the measured sum of these gases. 
Differences at higher altitudes are due primarily to 
HOCI, not measured in this experiment, and CINO 3 , 
not retrieved by ATMOS above ~34 km. Model C 
(solid lines) uses rates for the reactions of Cl + CH 4 
and OH + HCI from a reanalysis of laboratory data 
and allows for a 7% yield of HCI from the reaction of 
CIO + OH. [Adapted from Michelsen etal., 1995] 


able to further test the fidelity of 2-D predictive 
models and verify the chemical mechanism. 

Significant discrepancy between theory and observa- 
tion exists for inorganic chlorine in the middle and 
upper stratosphere. Figure 33 illustrates a compari- 
son of theory and observations for mid-latitude halo- 
gen partitioning [Michelsen et al., 1995]. HCI con- 
centrations are underpredicted while CIONO 2 are 
overpredicted. Measurements of CIO by the 
Microwave Limb Sounder (MLS) are consistently 
lower than predicted in the upper stratosphere. This 
suggests that the production and loss mechanisms of 
HCI in the middle and upper stratosphere are not 
well quantified. The effects of this discrepancy on 
model predictions of HSCT impacts have not yet 
been evaluated. 

Simultaneous measurements of hydrogen, nitrogen, 
and chlorine radicals show NO x catalytic cycles 
dominate removal of ozone for altitudes between 24 
and 38 km, as illustrated in Figure 34. At altitudes 
near 38 km, halogen cycles are less efficient sinks for 
ozone than calculated by models using accepted 


photochemistry [DeMore etal., 1994], resulting in a 
balance between production and loss of ozone [Jucks 
etal., 1995]. 

3.2.3.2 Lower Stratosphere Photochemistry 

The photochemistry of the lower stratosphere is 
somewhat distinct from the higher altitudes because 
the sulfate aerosol layer and PSCs provide sufficient 
surface area for heterogeneous chemistry to play an 
important role in determining the concentrations of 
radicals. Observations from the NASA ER-2 aircraft 
have provided the bulk of the data sets that have been 
used to test the basic photochemistry of the lower 
stratosphere put forward from laboratory studies. 
These observations have motivated the choice of test 
cases used in the 2-D benchmark chemistry test dis- 
cussed earlier in Section 2.3. 
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Figure 34. Production and loss of ozone in the 
middle stratosphere. Loss rates of ozone due to 
hydrogen radicals (closed triangles), nitrogen 
radicals (closed squares), halogen radicals (open 
circles), and the reaction of atomic O and ozone 
(open squares) inferred from calculated O, 
measurements of ozone, HO2, OH, NO, NO2, and 
CIO by the Far Infrared Spectrometer (FIRS-2) and 
the Balloon Microwave Limb Sounder (BMLS) near 
34‘N, on 26 and 27 September 1989. Error bars 
represent the l-o total accuracy of loss rates based 
on uncertainties associated with measurements of 
radicals, and the dotted lines connect individual 
measurements for convenience of interpretation. 
The production rate of ozone from photolysis of O2, 
calculated using the Minschwaner et al. [1993] 
formulation for transmission in the O2 Schumann- 
Runge region, is shown by the solid line. [Adapted 
from Jucks et al., 1995] 
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Figure 35. Diurnal variations of free radicals (OH, 
HO2, NO, NO2, and CIO) measured at 19 km, 37’N 
during SPADE. The lines are calculations using a 
photochemical model. [Adapted from Salawitch ef 
a/., 1994b] 

Detailed analysis of the SPADE data indicate that 
photochemical models using the existing kinetic data 
predict radical partitioning within ±40% [Salawitch 
et al., 1994a, b; Cohen et al . , 1994; Stimpfle et al., 
1994; Jaegl£ et al., 1994]. Figure 35 shows the ob- 
served diurnal variation of the concentrations of free 
radical species from the SPADE data set. Also 
shown is the calculation from a photochemical model 
constrained with measured or directly inferred values 
of aerosol surface area, NO y , ozone, temperature, 
pressure, and the ozone column. The agreement is 
generally better than expected from an error analysis 
of the measurements and models. The level of 
agreement shown in this figure is found consistently 
in analysis of all of the SPADE and most of the 
ASHOE/MAESA flights (except those with recent 
PSC processing or rapid latitudinal excursions) 
which include winter, spring, and fall observations 
over a large range in latitude and aerosol surface area 
loading. The Photochemistry of Ozone Loss in the 
Arctic Region in Summer (POLARIS) aircraft cam- 
paign, scheduled for summer 1997, is intended to fill 


the gap of summer observations when lower strato- 
spheric NO x concentrations are expected to be con- 
siderably higher. 

In general, the photochemical model underpredicts 
the concentration of OH and HO 2 . Most dramati- 
cally, the abrupt rise of OH and HO 2 at sunrise sug- 
gests that a nighttime process results in the accumu- 
lation of an easily photolyzed OH-containing 
molecule [Wennberg etal, 1994]. Recent laboratory 
observations of the rapid hydrolysis of B 1 ONO 2 on 
sulfate aerosols is a possible mechanism [Hanson 
and Ravishankara, 1995; Lary et al., 1995]. 
Inclusion of this process at the rate suggested by 
Hanson and Ravishankara improves the agreement 
between the model and measurement of OH 
throughout the day, particularly at high solar zenith 
angle. However, because this reaction is assumed to 
produce HNO 3 , the predicted NO x concentrations 
fall below the observations [R. Salawitch, personal 
communication]. Because this reaction leads to en- 
hanced concentrations of HO x and depresses NO x , it 
can impact the HSCT assessment. A sensitivity 
study has been performed within the 2-D model for- 
mulation. For a base case scenario (500 aircraft; 3.0 
ppbv Cl y ; EI NOx = 15; UNEP aerosol), the LLNL/UI 
2-D model predicts that global ozone change is re- 
duced from -0.66 to -0.49% when the process is in- 
cluded. Similarly, the CSIRO model predicts the 
change to be reduced from -0.44 to -0.33%. 

Measured noontime OH concentrations are also 
higher than expected, although the differences are 
not beyond the uncertainty in the measurements. 
Nevertheless, it has been suggested in laboratory 
studies, [Ball and Hancock, 1995], atmospheric mea- 
surements [Muller etal., 1995], and modeling work 
[Michelsen et al., 1994] that the DeMore et al. 
[1994] formulation of the O(’D) yields from ozone 
photolysis is not adequate. Addition of a larger yield 
of 0(!D) from ozone in the 310- to 315-nm range 
improves the agreement of the model with OH mea- 
surements [Salawitch et al., 1994a, b]. Although a 
sensitivity study to predict how this affects HSCT 
impact predictions has not been conducted, its inclu- 
sion should affect model predictions. The increase in 
HO x , particularly in the lower stratosphere, will fur- 
ther shift the photochemistry towards a HO ^domi- 
nated lower stratosphere. 

Observations of the NO/NO y ratio in the lower 
stratosphere are generally in good agreement with 
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the photochemical model. However, during the 
SPADE measurements, measured NO2 concentra- 
tions were considerably lower than predicted based 
on simultaneous measurements of NO, CIO, and 
ozone [Jaegld et al., 1994]. NO2 measurements 
during ASHOE/MAESA suggest a much smaller dis- 
crepancy on the order of 15% [Gao et al., 1995]. 
This problem deserves further attention because it 
will impact the HSCT assessment. As mentioned 
earlier, addition of the Br0N02 hydrolysis tends to 
lead to underprediction of the observed NO x /NO y 
ratio. This could arise from an error in measured 
surface area or the formulation of the N2O5 hydroly- 
sis reaction on sulfate aerosol, which dominates the 
heterogeneous conversion of NO x to NO y . 
Variations in the laboratory measurements for this 
process are uncertain up to a factor of 2. 
Additionally, the rate of this reaction is very depen- 
dent on the N2O5 formation rate which depends on 
the concentration of ozone and the temperature - ar- 
eas of known deficiencies in the 2-D models. 

Aircraft in situ observations of HC1 in the Southern 
Hemisphere in 1994 by the Aircraft Laser Infrared 
Absorption Spectrometer (ALIAS) instrument, ob- 
tained for near-background sulfate aerosol loading, 
agree well with model calculations (HCl/Cl y -0.8). 
Measurements of the HCl/Cl y fraction in the 
Northern Hemisphere have evolved upward from 
0.35 in 1991/2 to 0.45 in 1993 to 0.6 in 1994, com- 
pared to model values close to 0.8. For the 1993 and 
1994 data from both hemispheres, the measured 
HCl/Cl y fraction shows a possible dependence on 
aerosol particle surface area, suggesting an unidenti- 
fied heterogeneous loss of HC1 or Cl y that follows 
the evolution with time of the Mt. Pinatubo aerosol. 
From an ATMOS/ER-2 intercomparison in 
November 1994 in the Northern Hemisphere, the 
ATMOS and ALIAS measurements of HCl/Cl y agree 
within their respective uncertainties, although the 
ALIAS measurements are systematically lower by 
-10% at higher HC1 mixing ratios, and up to 40% at 
lower mixing ratios where precision is poor [Chang 
et al., 1995]. UARS measurements in 1992 and 1993 
also indicate that models overpredict the HCl/Cl y ra- 
tio at 550K, but not by as much as indicated by the 
ALIAS observations [Dessler etal., 1995b]. 

Observations of C10/Cl y are in good agreement with 
model calculations [Stimpfle et al., 1994]. The ob- 
served diurnal behavior of CIO indicates that it is 
strongly linked with photolysis of CIONO2; how- 


ever, the inferred fraction C10N02/Cl y = 0.06 to 
0.24 is in apparent conflict with the 1993 aircraft 
measurements of HC1. 

As discussed in Section 3.2.4, the aerosol chemistry 
in the polar regions remains quite uncertain. As 
changes in our knowledge of the aerosol micro- 
physics have occurred, laboratory investigations have 
been made of the rates of heterogeneous processes on 
various surfaces that are identified as potentially im- 
portant. These studies demonstrate that, largely in- 
dependent of the nature of the surface (e.g., liquid 
ternary mixture, NAT, sulfuric acid tetrahydrate 
(SAT)), the reaction of HC1 with CIONO2 becomes 
very efficient at low temperatures found in the winter 
polar regions (<195 K). This suggests that, with re- 
spect to the activation of chlorine, simulation of this 
chemistry in the 2-D models will be much more 
sensitive to the temperature and water than to the 
formulation of the microphysics (see Section 3.2.4). 

Recognition of the various model/observation dis- 
crepancies has led to speculation about possible 
missing chemistry in the models. Missing chemistry 
includes photochemical mechanisms that are not in 
the computer models, and errors in reaction rate 
constants and photolysis rates. These errors may oc- 
cur in the absolute values or in the dependencies on 
pressure, temperature, and the presence of other 
trace-gas constituents. For missing chemistry to be 
significant for this assessment, it must have chemical 
impacts comparable to known gas-phase and hetero- 
geneous chemistry and its inclusion in the models 
must change the calculated effects of stratospheric 
aircraft. Assessing the uncertainties associated with 
missing chemistry is difficult, as evidenced by the 
surprising discovery of the Antarctic ozone hole and 
by the N2O5 hydrolysis on sulfate aerosol that has 
widely altered the predicted impact of aircraft on the 
lower stratosphere. 

In addition to disagreements between model results 
and observations, the search for missing chemistry 
usually is motivated by new laboratory results. As 
laboratory gas-phase measurements have improved 
and become more consistent, much of the emphasis 
has shifted to photolysis rates and heterogeneous 
chemistry. To date, no systematic assessment has 
been made of the uncertainties associated with 
missing chemistry. 
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For a few instances, such as the HO x catalytic reac- 
tions, the possibility of missing chemistry has been 
constrained by observations [Cohen et al., 1994]. In 
situ measurements of the HO 2 /OH ratio appear to 
agree to better than implied by propagation of the 
quoted uncertainties in the measurements and the key 
reaction rate constants. Because the HO x catalytic 
cycle that destroys ozone is responsible for about one 
half of the total ozone loss at mid-latitudes, the con- 
straints on missing HO x chemistry provide important 
bounds on the overall uncertainty. 

Iodine chemistry has been suggested, and has re- 
ceived considerable attention, as a contributor to 
ozone loss in the lower stratosphere [Solomon et al. , 
1994]. However, observations of the iodine monox- 
ide (10) column abundance suggest that the inor- 
ganic iodine concentration in the lower stratosphere 
is lower than the 1 pptv assumed for the calculations 
[P. Wennberg, personal communication]. In addi- 
tion, recent measurements of the reaction rate con- 
stant for the proposed rate-limiting reaction that de- 
stroys ozone, 10 + CIO -* I + Cl + O 2 , is three times 
slower than had been assumed [Gilles et al., 1995] 
and the 10 + BrO -> I + Br + O 2 reaction is about as 
fast as assumed [Laszlo et al., 1995a, b]. Ongoing 
studies will better define iodine’s role in lower 
stratospheric ozone loss. However, if this or any 
other halogen reaction is significant, then any in- 
crease in N0 X from aircraft exhaust very likely will 
lead to a decrease in local ozone destruction due to 
increased control of the halogens by the added N0 X , 
provided that the interaction of the “missing” halo- 
gen oxide family with N0 X is similar to that of C10 X 
and BrO x . 

3.2.3.3 Upper Troposphere Photochemistry 

Changes in tropospheric ozone due to HSCT N0 X 
emissions as predicted by various models are highly 
variable and uncertain. The uncertainty is mostly 
due to the lack of understanding of the distribution 
and budget of N0 X in the troposphere, particularly in 
the upper and middle troposphere. The largest in situ 
source of N0 X in the upper and middle troposphere is 
from lightning. Yet neither the magnitude and alti- 
tude distribution are understood [Fehsenfeld and Liu, 
1993]. In addition, the convective transport of NO x 
emitted near the surface is important to the N0 X 
budget in the free troposphere but is extremely diffi- 
cult to evaluate. These near-surface sources are more 
than two orders of magnitude greater than the 


predicted N0 X emissions from HSCTs in the tropo- 
sphere. A small fraction of the near-surface N0 X 
transported to the upper troposphere can overwhelm 
the contribution of the NO x from HSCTs. 

Because the column density of ozone in the tropo- 
sphere is only -10% of the total column density, 
even with large uncertainty, the predicted changes in 
tropospheric ozone due to HSCT N0 X emissions 
have negligible impact on the total ozone column 
density. The same thing cannot be said about the ef- 
fect of N0 X emissions from subsonic aircraft because 
the emissions are significantly greater than those 
from HSCTs. Their effect on tropospheric ozone and 
climate will be addressed in a separate but coordi- 
nated NASA research program, the Subsonic 
Assessment (SASS) element of the Atmospheric 
Effects of Aviation Project (AEAP). 

3.2.3.4 Propagation of Errors Due to 
Uncertainty in Rate Coefficients 

One aspect of the predictive uncertainty that is 
amenable to a direct numerical estimate involves the 
reaction rate coefficients. Assuming that the chemi- 
cal model is complete, the estimated uncertainties 
from the JPL evaluation [DeMore et al . , 1994] have 
been propagated through a typical 2-D model calcu- 
lation [Stolarski and Considine, 1995]. Some 
preliminary results are shown in Figures 36 and 37. 
Figure 36 shows the calculated 1-a uncertainties on 
the calculated column ozone change due to a Mach 

2.4 HSCT fleet with EIno x = 15 as a function of lati- 
tude and time of year. The results show 0.5 to 0.8% 
uncertainties in the Southern Hemisphere, where the 
perturbation is due mainly to ozone-depleted air 
transported to the upper stratosphere (see Figure 23). 
In the Northern Hemisphere, the uncertainty rises 
rapidly to 1% and then more gradually to 1.5 to 2% 
near the pole. The larger Northern Hemispheric un- 
certainty in the column prediction is the result of 
larger uncertainties on the lower stratospheric pertur- 
bation, as shown in Figure 37. The smaller uncer- 
tainties in the upper stratosphere occur where the 
change in ozone is due primarily to an increase in the 
NO x catalytic cycle. In the lower stratosphere, the 
increased catalysis by NO x is balanced by interfer- 
ence and a subsequent decrease in the catalytic cy- 
cles due to C10 x and HO x . These processes, whose 
absolute uncertainties add together, result in the 
larger overall uncertainty. The next step in this anal- 
ysis is to identify and remove from consideration 
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Figure 36. Calculated 1 -ct uncertainty in column 
ozone change (in percent) for a Mach 2.4 HSCT 
fleet with EIno x = 15 g N02/kg fuel from the GSFC 
model Monte Carlo error analysis as a function of 
season and latitude. 

those model results that predict values of key species 
which are outside the range of measured concentra- 
tions, thereby narrowing the uncertainties. 

3.2.4 Aerosol and Particle m icrophysics 

3.2.4. 1 Aerosol Formation Processes 

The 2-D HSCT predictive models include sulfate 
aerosols prescribed by the UNEP surface area distri- 
bution [WMO, 1992], which represents an average 
derived from satellite observations in 1989, prior to 
the Mt. Pinatubo eruption. Since the variability due 
to volcanic sulfur injection or PSC formation is not 
considered, the evolution and transformations of 
particles are not complete in the models. This omis- 
sion is particularly important for low-temperature 
regions at high latitudes where both particle micro - 
physics and heterogeneous chemical reactivity have 
large temperature dependencies. Separate sensitivity 
studies that include a parameterization for PSCs, 
variations in sulfate aerosols due to volcanoes, and 
HSCT particle formation are discussed in Sections 
3.2.4.5 through 3.2.4.7. 

The composition and abundance of aerosols in the 
stratosphere are controlled by the vapor densities of 
H 2 SO 4 , HNO 3 , and H 2 O that condense into particles 
according to their temperature-dependent ternary 
phase diagram. In most of the low- and mid-latitude 
stratosphere at temperatures above 200 K, super- 
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Figure 37. Calculated 1 -a uncertainty in ozone 
concentration (in percent) for a Mach 2.4 HSCT fleet 
with EInox = 15 g N02/kg fuel from the GSFC model 
Monte Carlo error analysis as a function of altitude 
and latitude for the month of March. 

cooled liquid H 2 SO 4 /H 2 O particles form the ambient 
aerosol layer observed between 10 and 30 km. At 
lower temperatures (<200 K), HNO 3 condensation 
forms ternary particles that become Type I PSCs 
composed largely of HNO 3 and H 2 O at temperatures 
below 192 to 195 K [McCormick etal., 1982; Steele 
et al., 1983]. At still lower temperatures (<190 K), 
H 2 O condensation forms Type II ice clouds. A com- 
plete model of stratospheric aerosols requires a de- 
tailed microphysical formulation of the formation 
and evolution of all these particles [see e.g., Turco et 
al., 1989]. 

Such a model would include sources and sinks of 
gaseous sulfur species. Primary stratospheric 
sources are transport of carbonyl sulfide (OCS) from 
the troposphere and episodic volcanic injections of 
SO 2 . In the stratosphere, all these species are oxi- 
dized to H 2 SO 4 . Microphysical models must ac- 
count for the competition between nucleation of new 
H 2 SO 4 particles and accommodation on preexisting 
particles. In this context, the aircraft exhaust plume 
forms a special case where high acid and water vapor 
levels favor binary nucleation of new particles. 
Models also include the effects of coagulation on 
aerosol concentration and surface. Finally, the pri- 
mary loss mechanism is transport via sedimentation 
to the troposphere. 

In situ and remote observations of PSCs [for exam- 
ple, Dye et al, 1992; Browell et al., 1993; Toon and 
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Tolbert, 1995] show considerable variability that is 
consistent with the complex ternary phase diagram. 
These PSCs are often composed of liquid and solid 
phases that are distinct from the most stable thermo- 
dynamic phase, presumed to be crystalline NAT 
[Hanson and Mauersberger, 1988]. For temperatures 
above 195 K, H 2 O condenses onto particles consis- 
tent with supercooled liquid sulfate composition in 
range of 50-80 weight percent [Steele and Hamill, 
1981]. At temperatures of -192 K, HNO 3 and H 2 O 
can condense to form a ternary liquid solution whose 
composition depends on the overall HNO 3 /H 2 SO 4 
ratio in the air [Molina et al., 1993; Tabazadeh et al., 
1994]. The existence of supercooled liquid particles 
[Dye etal., 1992; Drdla et al., 1994; Carslaw etal., 

1994] is consistent with much laboratory and theoret- 
ical work which predicts significant barriers to freez- 
ing of both binary H 2 SO 4 H 2 O and ternary 
H 2 SO 4 HNO 3 H 2 O liquid particles [Song, 1994; 
Koop et al., 1995; Anthony etal., 1995; MacKenzie 
etal., 1995; Murphy and Gary, 1995]. Moreover, 
other laboratory studies show multiple metastable 
solid phases are possible upon freezing [Worsnop et 
al., 1993; Marti and Mauersberger, 1993; Fox et al., 

1995] . 

Original models that identified Type I PSCs as NAT 
with a temperature threshold of -195 K cannot ex- 
plain Type I PSC appearance between 191 and 192 K 
[Dye et al., 1992; Kawa et al., 1992; Drdla et al., 
1994; Carslaw etal., 1994]. These complexities 
raise questions about the validity of relating tempera- 
ture and global PSC probabilities, such as those cata- 
logued by Poole and Pitts [1994], to heterogeneous 
chemical processes such as those controlling chlorine 
activation which have steep dependence on aerosol 
temperature, composition, and phase [see e.g., 
Hanson and Ravishankara, 1993; Hanson et al., 
1994]. 

Overall, uncertainties in aerosol microphysics in- 
crease at high latitudes, and low temperatures must 
be combined with steeply temperature-dependent 
microphysical and heterogeneous reaction processes 
to predict effects on photochemistry. This is particu- 
larly true for microphysical aerosol processes for 
which there is significant conceptual uncertainty. 
Most importantly, mechanisms of denitrification 
likely depend on transformations between liquid and 
solid metastable phases. More study is needed to 
evaluate the microphysical implications of adding 


HSCT plume sulfate particles to background aerosol 
distributions on these PSC/denitrification processes. 

3.2.4.2 Denitrification and Dehydration 

Field experiments have shown that heterogeneous 
chemical processes during polar winter in both the 
Antarctic and Arctic yield very high levels of active 
chlorine, yet widespread ozone depletion has oc- 
curred to date only over Antarctica. This is thought 
to be due to the fact that the springtime Antarctic 
stratosphere is substantially denitrified, its NO x con- 
tent being irreversibly reduced over large areas by 
some 90% relative to both prewinter and extra- vortex 
values [Crutzen and Arnold, 1986; Solomon et al., 
1986a; McElroy et al., 1986, Fahey et al., 1990]. 
Denitrification has also been observed in the Arctic, 
but only to a much smaller degree and in localized 
“patches.” The springtime Antarctic (but not the 
Arctic) stratosphere is also substantially dehydrated, 
its H 2 O content over large areas being less than half 
that of prewinter and extra-vortex values. 

The presence or absence of NO x is crucial to deter- 
mining the amount of chlorine-catalyzed ozone de- 
struction. If sufficient NO x is present in polar 
springtime, chlorine radicals are converted quickly to 
the reservoir CIONO 2 and more slowly to the reser- 
voir HC1. If there is substantial denitrification, the 
reservoirs remain depleted, and ozone destruction 
can proceed rapidly until nearly all the ozone is de- 
stroyed or until the polar vortex breaks down. 

Denitrification and dehydration are caused by the 
precipitation of large PSC particles that contain con- 
densed HNO 3 and H 2 O, respectively. Microphysical 
models of Type 2 (ice) PSC formation can produce 
particles that are large enough and precipitate fast 
enough to explain observed dehydration levels. 
Widespread dehydration occurs in the Antarctic, but 
not in the Arctic, because Antarctic stratospheric 
temperatures often fall below the frost point for ex- 
tended periods of time [Fahey et al., 1990]. 
However, the mechanism by which denitrification 
occurs remains largely a mystery. It is difficult to 
explain observed levels of denitrification through 
“direct” microphysical models of either Type 1 or 
Type 2 PSC formation [see e.g., Salawitch et al., 
1988; Toon et al., 1990]. Type 1 PSC models gen- 
erally produce particles that contain substantial 
amounts of condensed HNO 3 , but which are too nu- 
merous and hence too small to precipitate at rates re- 
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quired to explain observed denitrification levels. 
Models assuming that Type 2 PSCs nucleate directly 
on preexisting Type 1 PSCs can produce precipitat- 
ing particles, but the particles contain too little HNO3 
to account for observed denitrification levels. Other 
denitrification schemes have been proposed, such as 
the independent scavenging and irreversible removal 
of gas phase HNO3 by falling Type 2 particles 
[Wofsy et al . , 1990], but these are supported by lim- 
ited observations [Deshler et al., 1994; Peter et al., 
1994]. 

A concern with the projected HSCT fleet is whether 
NO x , H2O, and/or particulate emissions from the 
fleet might increase P SC formation in the Arctic. 
PSC formation could produce more widespread deni- 
trification and/or dehydration and hence trigger 
widespread Arctic ozone destruction. This is even 
more of a concern if ambient temperatures in the 
lower stratosphere in the HSCT era are lower than at 
present, as predicted by global climate models. Even 
without HSCT emissions, lower temperatures alone 
would increase the extent of PSC formation in the 
Arctic. There have been attempts to represent PSC 
microphysical processes in 2-D predictive models 
(see Section 3.2.4.6). However, since the basic 
mechanism for denitrification is unknown, one must 
conclude that the impact of HSCT emissions on 
Arctic PSC formation and denitrification cannot be 
predicted confidently. Understanding the mecha- 
nisms controlling denitrification will probably re- 
quire detailed, process-oriented field measurements 
of temperature, HNO3, H2O, and particle size and 
composition. Since it is problematic to predict when 
and where denitrification will occur in the Arctic, the 
Antarctic appears to be the best place to collect the 
required observations. 

3.2.4.3 Plume/Wake Aerosol Processing 

Input for the HSCT fleet emission 2-D predictive 
models is achieved by simple dilution of the aircraft 
plume at the altitude of injection, with no chemical 
changes taking place, into a 1° longitude by 1° lati- 
tude by 1 km altitude grid. Omission of plume ef- 
fects from the model input is supported by 
plume/wake model predictions for NOy gas-phase 
species which show no significant conversion from 
NO x to other NOy species in that regime. However, 
near-field processing of SO x emissions and near- 
field aerosol formation may not be similarly dis- 
missed (see Section 3.1). Model results demonstrate 


that a large number of small particles can form via 
binary homogeneous nucleation of sulfur compounds 
[Miake-Lye et al., 1994; NASA, 1993; Karcher, 
1995; Zhao and Turco, 1995; Brown et al., 1995], 
This has been confirmed by measurements in the 
Concorde plume which show that large numbers of 
small, persistent particles are rapidly formed in the 
engine or in the plume [Fahey et al. , 1995a]. 

Comparison of the particle mass measured in the 
Concorde exhaust with EIso2 inferred from fuel 
analysis [Fahey et al., 1995a] indicates that at least 
12 to 45% of the sulfur emitted from the aircraft is in 
the form of small particles less than one hour after 
emission. This result is consistent with earlier stud- 
ies of plumes behind a stratospheric SR-71 aircraft 
[Hofmann and Rosen, 1978] and subsonic jets 
[Pitchford et al., 1991; Whitefield and Hagen, 1995], 
However, these results conflict with the plume 
chemistry model that predicted only 1% of the sulfur 
being converted to particles, based on OH-limited 
oxidation of SO2 to H2SO4 [Miake-Lye et al. , 1994; 
Brown et al., 1995]. Measurements of OH in the 
Concorde exhaust support the plume model con- 
tention [Fahey et al., 1995a]. Therefore, the 
Concorde particulate measurements apparently indi- 
cate a rapid conversion of SO2 to small sulfate parti- 
cles via some unknown process. Together with the 
observation of small numbers of sulfate particles in 
an ER-2 plume encounter [Fahey et al. , 1995b], these 
results show that current understanding of plume 
chemistry and particle formation processes is not ca- 
pable of predicting particle emission rates for differ- 
ent stratospheric aircraft. 

The potential impact of HSCT-produced sulfate 
aerosols on ozone (see Section 3.2.4. 5) raises the 
possibility of the need to reduce sulfur levels in the 
fuel. The sulfur content of jet fuel is currently 
-0.04% by weight (EIs02 = 0.8), much less than the 
current maximum specification of 0.3%, and is ex- 
pected to diminish in the future [Hadaller and 
Momenthy, 1989, 1993]. One other unknown is 
whether soot plays any role in the particle formation 
process (see Section 3.2.4.7). The Concorde plume 
encounter showed that one third of the particles had 
an in volatile core (larger than 10 nm in diameter) 
that is most likely derived from soot emitted from the 
aircraft. It is clear that prediction of plume particle 
emissions from stratospheric aircraft plumes requires 
further modeling and field study of near-field aircraft 
chemical and microphysical processes. 
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3.2A.4 Sensitivity of HSCT Ozone 

Perturbations to Ambient Sulfate 
Aerosol Loading 

Injection of SO 2 from major volcanic eruptions over 
the past 20 years has caused episodic stratospheric 
sulfate aerosol mass loading increases of more than a 
factor of 10 [McCormick et al., 1995]. These pertur- 
bations remain significant for three or more years; 
post-eruption total stratospheric aerosol mass de- 
creases with a 1 /e-folding time of approximately one 
year. It is difficult to define the true nonvolcanic 
(background) aerosol level because of the frequency 
of major eruptions, their lingering influence 
[Thomason et al., 1995], and the possibility of a 
positive temporal trend linked to increasing subsonic 
aircraft emissions [Hofmann, 1990]. 

Fahey et al. [1993] presented ER-2 data illustrating 
the impact of heterogeneous processing by sulfate 
aerosols on mid-latitude NO y chemistry. The ob- 
served NO x /NO y ratio decreased by a factor of 2 or 
less as aerosol surface areas increased from 0.5-5 
[tm 2 cm -3 to -20 pm 2 cm-3. This result indicates a 
nonlinear dependence of the NOy partitioning into 
reactive and reservoir species on surface area. The 
ratio of NO x to NO y should approach a low, constant 
value at surface areas of 5 to 10 times the nonvol- 
canic background. 

Model calculations of the sensitivity of HSCT-in- 
duced chemical perturbations to ambient sulfate sur- 
face area have been performed by Weisenstein et al. 
[1993], Considine etal. [1995], Tie etal. [1994], and 
Ramaroson and Louisnard [1994]. The first two of 
these studies used the surface area distribution given 
in WMO [1992] to represent ambient conditions, and 
multiplied this distribution by a factor of 4 to repre- 
sent the median aerosol loading observed over the 
past two volcanically active decades. Thomason and 
Poole [1995] show that the WMO distribution is very 
similar to the lowest 10-percentile SAGE II surface 
area distribution over the last decade, while the 
decadal median is, on the average, a factor of 2 to 3 
higher than the 10-percentile distribution. Results of 
the Weisenstein et al. [1993] (Figure 38) and 
Considine et al. [1995] studies are similar, indicating 
that ozone depletion due to HSCT emissions is re- 
duced as the ambient sulfate aerosol surface area is 
increased. It must be noted that these calculations 
assume that the ambient aerosol surface area is not 
affected by the HSCT emissions of sulfate or soot. 


The Tie et al. [1994] model included explicit sulfate 
aerosol microphysics to simulate the global distribu- 
tion of aerosols produced by an eruption such as El 
Chichon. They found that during the post-volcanic 
period, the injection of NO x from the HSCT fleet 
would actually increase the ozone column by inter- 
fering with the catalytic destruction of ozone by 
C10 x . 

3.2.4.5 Sensitivity of HSCT Ozone Perturbations 
to Fleet Sulfur Emissions 

An important issue for a full assessment of the pro- 
posed HSCT fleet is whether the sulfur present in the 
aircraft exhaust could lead to sulfate aerosol surface 
area increases large enough to impact ozone. 
Sensitivity studies addressing this issue, including 
explicit sulfur chemistry and microphysics, have 
been performed by Bekki and Pyle [1993], Pitari et 
al. [1993], and Tie et al. [1994], Using EIso 2 = K 
two of these studies predict surface area increases of 
20 to 50% relative to the nonvolcanic background, 
with a small decrease in calculated ozone depletion 
due to the HSCT emissions. Bekki and Pyle [1993] 
predict a doubling of aerosol surface area in the 
northern lower stratosphere due to HSCT flights, 
probably because they assumed an El for sulfur of 
2 g/kg fuel. As would be expected, they also show a 
more substantial reduction in the sensitivity of ozone 
to the HSCT NO x emissions, accompanied by an in- 
creased sensitivity of ozone to chlorine. These stud- 
ies all assumed that the sulfur was emitted by the en- 
gines as SO 2 so that it would condense onto preexist- 
ing particles. 

A recent study by Weisenstein et al. [1995] used a 
2-D model with coupled sulfur chemistry and aerosol 
microphysics to examine the combined effects of 
NO x , sulfur, and water vapor emitted by HSCTs. 
The results are summarized in Table 5. Of particular 
interest are the results in the last row, where all the 
sulfur is assumed to be emitted as small particles or 
to form small particles rapidly in the aircraft wake. 
The impact of small particles on the calculations de- 
pends in a complicated way on EIno x and back- 
ground chlorine. For EIno x = 5 g/kg fuel, the calcu- 
lation of ozone change is more negative when small 
particles are formed; the perturbation is as large or 
larger than that for the EIno x = 15 g/kg fuel case. 
Inclusion of the reaction of B 1 ONO 2 with water on 
aerosols increased the effect of the small particles. 
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Figure 38. Calculated change in column amount of 
ozone as a function of month and latitude for a Mach 
2.4 HSCT fleet with EIno x = 15 g NC> 2 /kg fuel from 
the AER model; a) aerosols held fixed at 
background as defined in WMO [1992], and b) 
aerosols fixed at four times background (approach 
discussed in Weisenstein etal. [1993]). 

This occurs because the increased surface area 
changes the mid-latitude chlorine partitioning in the 
model. It is important to note that the Weisenstein et 
al. [1995] calculation represents the extreme case. 
This calculation assumed that all of the emitted sul- 
fur formed persistent small particles, thus maximiz- 


ing the perturbation in surface area. (See Section 
3. 2.4.3 for discussion of uncertainty regarding this 
point.) Another assumption was that the perturbation 
was added to a non volcanic aerosol background pro- 
duced by photolysis of tropospheric source gases 
only. If the stratospheric sulfate layer were already 
perturbed by volcanic activity, then saturation of the 
N 2 O 5 hydrolysis reaction should reduce the sensitiv- 
ity to aerosol surface area added by the HSCT. Thus, 
more attention must be given to the role of HSCT 
exhaust products on changing the background 
aerosol surface area in a future stratosphere contain- 
ing lower amounts of chlorine. 

3.2.4.6 Sensitivity of HSCT Ozone Perturbations 
to Inclusion of PSCs 

Heterogeneous reactions on PSCs cause dramatic 
chemical perturbations in the polar regions, and there 
is concern that HSCT exhaust products might lead to 
an increase in PSC occurrence or to changes in their 
microphysical characteristics which would ultimately 
affect ozone. The models used for the basic HSCT 
perturbation calculations do not include PSCs, be- 
cause of microphysical uncertainties and also be- 
cause of fundamental limitations in the 2-D represen- 
tation. For example, data from the Stratospheric 
Aerosol Measurement (SAM II) show that PSC 
sighting probability is not zonally symmetric, espe- 
cially in the Arctic, with higher probabilities occur- 
ring in lower-temperature regions [Poole and Pitts, 
1994], Likewise, it is difficult for 2-D models to 
treat variations in sunlight experienced by a PSC- 
processed air parcel as it moves around the polar 
vortex when the vortex is displaced or elongated. 
Finally, it is difficult for these models to predict the 
transport (or lack thereof) to lower latitudes of air 
parcels that have been chemically processed by PSCs 
in the polar vortices. 

Nonetheless, several groups have used 2-D models 
that include simplified representations of PSC pro- 
cesses to investigate the impact of HSCT emissions 
on stratospheric ozone [Pitari et al., 1993; Considine 
et al., 1994; Tie et al., 1994]. All the models as- 
sumed that Type 1 PSCs form as NAT, but they in- 
corporated different representations of the cloud par- 
ticles and temperature excursions. Pitari etal. [1993] 
included an explicit particle microphysics code with 
PSC nucleation delayed until reasonable levels of 
vapor supersaturation were reached. The tempera- 
ture at each model grid point was treated statistically 
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Table 5. Calculated change in annual average ozone column at 57“N due to emissions of a Mach 2.4 
HSCT fleet including the effects of NO x , H20, and sulfur. Results shown are for EISO 2 = 0.4 9 P er k 9 fuel 
burned and EIno x = 5 and 15. All calculations have been done for background inorganic chlorine amounts 
of 3 ppbv and 2 ppbv. [Adapted from Weisenstein eta!., 1995] 


SO 2 Emission 
(EISO 2 = 0*4) 

Aerosol Surface 
Area Change 

ElNOx = 5 

ElNOx = 15 

3 ppbv 

2 ppbv 

3 ppbv 

2 ppbv 

UNEP aerosols a 

0% 

-0.4%" 

-0.6% 

- 1.1% 

-1.6% 

0% SO2 as particles 

30% 

-0.5% 

-0.6% 

-0.7% 

-1.1% 

10% SO2 as particles 

50% 

-0.7% 

-0.7% 

-0.7% 

-1.0% 

100% SO2 as particles 

150% 

-1.0% 

-0.8% 

-0.8% 

-0.9% 


a UNEP aerosol refers to an atmosphere with a background stratospheric aerosol amount equivalent to a 
relatively clean period between major volcanic eruptions. None of these scenarios included the effects of polar 
stratospheric cloud chemistry. [WMO, 1992] 


by assuming a normal distribution, with a mean tem- 
perature equal to the monthly zonal average value 
and a standard deviation derived from global-atmo- 
spheric circulation statistics from 1958 to 1983. PSC 
condensation probabilities calculated using this ap- 
proach appear to be somewhat lower than the clima- 
tological SAM II data. Considine et al. [1994] repre- 
sented grid-point temperatures by empirical proba- 
bility distributions derived from National 
Meteorological Center (NMC) data for the period 
January 1989 to May 1992. This model also required 
modest supercoolings before PSCs were allowed to 
form. PSC particle surface areas were determined by 
allowing HNO3 and H2O to condense onto log-nor- 
mal distributions with fixed-mode radii of 1 pm for 
Type 1 PSCs and 10 pm for Type 2 PSCs, and a ge- 
ometric standard deviation of 1.8 for both cloud 
types. Cloud formation probabilities calculated by 
Considine et al. [1994] agreed fairly well with the 
SAM II data. 

Both the Pitari et al. [1993] and the Considine et al. 
[1994] studies found that inclusion of PSCs generally 
reduced the sensitivity of model ozone to HSCT 
emissions. For example, Pitari et al. [1993] found 
that HSCT-induced Northern Hemisphere mid-lati- 
tude ozone depletion decreased from 1.5% to 0.5% 
when the PSC formulation was included in the 
model. This decrease is due primarily to enhanced 
NOy removal by the PSCs, which shortens the resi- 
dence time of the NO x emitted by the HSCT fleet. 


However, the NOy removal rates computed in these 
models may be artificially large due to their simpli- 
fied particle microphysics representations. Both 
studies noted an important exception to this reduced 
sensitivity in Antarctic spring, when larger PSC sur- 
face areas caused by the HSCT emissions extend the 
duration of PSCs into late October and November. 
This results in an extended time period in which the 
chlorine catalytic cycle operates efficiently and, thus, 
leads to more ozone loss than in model runs exclud- 
ing PSCs. Other calculations used Type 1 PSC mode 
radii of 0.5 pm and 2 pm (a crude test of NOy re- 
moval rates); these yielded global ozone depletions 
that differed only slightly from that of the standard 
run. The second sensitivity study used both 1980 
and 1990 atmospheric chlorine amounts (-2.0 and 
-3.5 ppbv, respectively). With inclusion of the PSC 
parameterization, the higher chlorine level reduces 
the sensitivity of ozone to HSCT emissions, with the 
exception again of the Southern Hemisphere high- 
latitude spring, where more ozone depletion was 
predicted to occur. 

The 2-D model study by Tie et al. [1994] allowed 
PSCs to form at temperatures several degrees K 
above their thermodynamic thresholds to account for 
the fact that zonal-mean temperatures are somewhat 
higher than those where PSCs are most likely to 
form. This study did not include a microphysics rep- 
resentation, but instead prescribed time constants for 
the pertinent heterogeneous reactions and for NO y 
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and water vapor removal. PSC-formation tempera- 
tures were raised by an additional 1 K to account for 
HSCT-induced increases in HNO3 and H2O vapor. 
Tie et al. [1994] found a substantial (6%) reduction 
in ozone column abundance at high northern lati- 
tudes in late winter and early spring, but this result is 
highly uncertain due to the simplicity of their PSC 
representation parameterization. 

Pitari and Visconti [1994] used their 2-D model to 
evaluate the effect of an HSCT fleet in an atmo- 
sphere where the CO2 concentration was increased 
from 335 to 500 ppmv. The lower temperatures in 
this atmosphere led to a large enhancement in NAT 
surface areas in the Arctic and, hence, to more NO y 
removal and a reduced effect of the added NO x from 
an HSCT fleet. They also used a 3-D circulation 
model to calculate a change in the transport rates in a 
CC>2-perturbed atmosphere. They found that the 
stratospheric residence time of HSCT emissions was 
15% shorter than in an atmosphere with the present 
amount of CO2. This leads to less buildup and, con- 
sequently, smaller perturbations in NO x concentra- 
tions. 

All these sensitivity studies indicate that including 
the effects of PSCs may be important for improving 
the accuracy of the 2-D predictive models. The rep- 
resentation of PSC formation in the 2-D models by 
statistical temperature fluctuations about the zonal - 
mean values reasonably matches observed PSC oc- 
currences, but the photochemical results have not 
been tested against atmospheric observations. 
Fundamentally, the 2-D representation, and the rela- 
tively simple microphysical assumptions included 
within it, contrast with the complex and variable be- 
havior inferred from the laboratory studies and field 
observation of PSCs. Thus, confidence in the ulti- 
mate results of the models - ozone depletion resulting 
from HSCT emissions - would increase with a more 
accurate treatment of PSC microphysics, including 
an assessment of uncertainties in particle composi- 
tion and denitrification mechanisms. 

3.2.4.7 Impact of Soot on Aerosols 

Although an extensive database exists for conven- 
tional aircraft soot emissions, it is comprised of 
“smoke number” data, an industry-standard measure 
that can be related to total mass but not to particle 
size or density. Since the atmospheric effects of par- 
ticles depend on their effective surface area and re- 


activity, the total soot mass is not particularly useful 
and does not effectively constrain atmospheric model 
usage of soot emissions. In addition, the plume evo- 
lution of the soot particles is highly uncertain, in- 
volving competition between nucleation, gas-to-par- 
ticle condensational growth, and coagulation. In 
view of these uncertainties and the lack of laboratory 
data on direct gas/soot interactions, soot is not con- 
sidered explicitly in the assessment model scenarios. 

Recent experimental data have been collected from 
ground-based jet engine tests and in situ plume en- 
counters [Hagen et ai, 1993; Whitefield et aL , 1993; 
Whitefield and Hagen, 1995; Pitchford et ai, 1991; 
Baumgardner and Cooper, 1994; Schumann et ai y 
1995; Fahey et ai , 1995a, b]. The general conclu- 
sions of such tests are that nascent soot emissions 
follow a log-normal size distribution with a peak 
around 0.05 pm. Total number densities are on the 
order of 10 6 cm-3, resulting in soot Els on the order 
of 0.01 g/kg fuel. Latitude surveys of soot abun- 
dance in the lower stratosphere and upper tropo- 
sphere also have been made using wire impactors 
exposed at altitude [Pueschel et aL, 1992; Blake and 
Kato, 1995]. These show that the latitudinal distri- 
bution of soot at 10 to 11 km altitude between 90°N 
and 45 °S covaries with commercial air traffic fuel 
use. Lower stratospheric soot measurements in the 
40°N to 60°N latitude band range from 0.25 to 
3.35 ng nr 3 , averaging ~2 ng nr 3 , while those in the 
Southern Hemisphere range from 0 to 0.6 ng nr 3 , 
with an average of 0.1 ng nr 3 . 

Collectively, these measurements suggest that HSCT 
soot emissions could possibly increase the natural 
stratospheric soot background by a factor of 2 to 3. 
It should be noted that the natural stratospheric soot 
mass loading is only ~1% of that of stratospheric sul- 
fate. However, under the extreme assumption that all 
of the emitted soot particles become activated as CN 
for H2O and H2SO4, the resulting sulfate aerosol sur- 
face area will rival that of the ambient background. 
Potential CN activation is particularly important in 
the near-field region behind the aircraft. It is also 
possible that HSCT-emitted soot particles will serve 
as nuclei for cirrus cloud formation and potentially 
alter cirrus properties or formation mechanisms. 

3.3 Analysis of Uncertainties 

In Chapter 2, a number of 2-D models were used to 
predict the impact of a fleet of HSCTs on ozone. For 
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a Mach 2.4 fleet of 500 airplanes with an assumed 
EInOx = 5 g/kg fuel, the models calculated a pertur- 
bation to the total column amount of ozone averaged 
over the Northern Hemisphere which ranged from 
+0. 1 % to -0.3%. These calculations were made with 
models using the reaction rates and photolysis cross 
sections recommended in the 1994 JPL kinetics re- 
view [DeMore et al . , 1 994]. Results of similar calcu- 
lations with the same models for higher Els of NO x 
(10 and 15 g/kg fuel) and for cruise altitudes of 1 km 
above and 2 km below the nominal Mach 2.4 cruise 
altitude are summarized in Chapter 2. 

The calculations are uncertain for a number of rea- 
sons discussed above. These uncertainties can be 
classified in two ways: those which may be evalu- 
ated directly through parametric studies in an exist- 
ing model (e.g., the impact of uncertainties in input 
reaction rate coefficients), and those for which we do 
not have sufficient information to do this completely 
(e.g., model representation of transport and the effect 
of PSCs). The first class of uncertainties are quan- 
tifiable at a specified confidence level; the second 
class is not yet quantifiable, but still affects the 
overall confidence of the prediction. The results of 
parametric studies, which have been completed as a 
Transport uncertainties have proven to be more diffi- 
cult to quantify. Although there has not been a sin- 
gle comprehensive study of the sensitivity of the 
prediction of HSCT effects to dynamical uncertain- 
ties, there have been systematic efforts to evaluate 
some individual uncertainties as they have been 
identified. 

Stratosphere-troposphere exchange is very important 
for the fate of exhaust emissions in the lower strato- 
sphere, because it determines where and how fast the 
polluted air exits to the upper troposphere. The scale 
and type of dynamical processes which are responsi- 
ble are either too small (e.g., tropopause folds near 
jet streams) or longitudinally dependent (e.g., the 
upward extension of monsoonal circulations) to be 
represented in 2-D models. Such uncertainties will 
be difficult to quantify because even state-of-the-art 
3-D models have not been tested against 
observations in this regard. For the current 
assessment, the five 2-D models show very little 
difference in the total perturbed amount of NO y 
(range from 3.5 to 4.3 x 10 15 molecules cm -2 in the 
global-average column). The 3-D model described 
in Section 3.2.2 calculated perturbed NO y within this 
same range. These calculations are consistent with a 


part of this assessment report, are summarized here. 
Where parametric studies are incomplete or unavail- 
able, the processes which lead to uncertainty in the 
models are identified for future work. 

The uncertainties in gas-phase reaction rate coeffi- 
cients have been propagated through one of the 2-D 
predictive models using the Monte Carlo approach. 
This resulted in an uncertainty of approximately 1% 
(1 o) averaged over the Northern Hemisphere. 
Inclusion of uncertainties in the input cross sections 
for photolysis and heterogeneous reactions would in- 
crease the range of calculated values. Comparing 
the results for many model cases for the current at- 
mosphere with atmospheric measurements of key 
radical concentrations narrows the range for the cal- 
culated ozone change because not all of the model 
cases are consistent with these measurements. 

The effect of B 1 ONO 2 hydrolysis, while absent from 
the suite of 2-D predictive models, has been evalu- 
ated and makes the predicted perturbation less nega- 
tive by a few tenths of a percent. This is consistent 
with our understanding of the role of this reaction in 
converting NO x to HNO 3 in low-temperature strato- 
spheric regions. 

stratospheric residence time for aircraft exhaust of 13 
to 16 months. 

A further test was to use dynamical inputs from the 
GSFC model in the AER model. The resulting per- 
turbation to ozone was less negative than the AER 
value, but more negative than the GSFC value. This 
indicates that a significant portion of the difference 
between the models is due to the representation of 
dynamics. However, the differences were not due to 
the total steady-state accumulation of HSCT pollu- 
tants. Rather, the dynamics in the separate 2-D 
models create differences in constituents such as 
ozone and NO y that are affected by transport in the 
model base atmosphere and differences in transport 
of pollutants to higher altitudes. 

The question of the effect of the zonal asymmetry of 
the pollutant source has been investigated using a 
3-D calculation. The maximum accumulations in 
flight corridors differ from the zonal mean by only a 
factor of 2. Since the photochemistry is expected to 
be linear over this range, the uncertainty in the 2-D 
calculation introduced by the “corridor effect” is 
small compared with other uncertainties. 
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The transport between the tropics and mid-latitudes 
is poorly represented in current 2-D models. The 
importance of the tropical pipe concept was tested 
through arbitrary restriction of subtropical diffusion 
coefficients in the AER model. This produced a cal- 
culated perturbation which was slightly more posi- 
tive than the result for AER in Chapter 2 (see Table 
4). This is qualitatively consistent with the idea that 
the barrier restricts the transport of aircraft effluent 
from mid-latitudes to the tropics, allowing efficient 
transport to high altitudes where NO x catalysis of 
ozone destruction is rapid. However, the assumed 
altitude dependence of the barrier suggests that this 
study does not represent the extreme case because of 
competing effects. The existence of a barrier to 
transport between the tropics and mid-latitudes re- 
duces the amount of exhaust gases which can be 
transported to high altitudes. The worst-case sce- 
nario for ozone depletion would be a leaky barrier in 
the lower stratosphere, allowing exhaust gases from 
the mid-latitudes to reach the tropics, and a stronger 
barrier above, effectively confining the gases and 
allowing them to reach high altitudes. This case 
should be evaluated. 

As discussed in Section 3. 2.4.6, the 2-D models used 
in this assessment do not include the effects of PSCs. 
The important processes are inherently difficult to 
represent in a 2-D framework. Several groups have 
used models which attempt to overcome these diffi- 
culties, such as the nonuniform distribution of winter 
temperature with longitude. Over most of the globe, 
the model response is generally less negative because 
models with PSCs have more chlorine-catalyzed 
ozone destruction which is damped by addition of 
NO x . At high-latitude spring, the model response is 
more negative due to increased persistence of PSCs 
as a result of more HNO 3 and H 2 O. The formation 
of PSCs is a strong function of temperature, H 2 O, 
and HNO 3 , and the chemical response is nonlinear 
both as a result of the threshold process and model 
sensitivity to sedimentation of PSC particles, which 
permanently remove NO y . Calculations to date do 
not represent the range of possible answers for the 
impact of PSCs on the HSCT perturbation. 
However, it should be noted that the importance of 
chlorine-catalyzed ozone destruction resulting from 
PSC processes will decrease as chlorine from man- 
made gases decreases during the next century. 

A recent study [Weisenstein et al., 1995] considered 
the impact of sulfur in the exhaust on the strato- 


spheric sulfate particle surface area. The impact de- 
pends critically on the form of the sulfur in the ex- 
haust. If small particles are formed, as indicated in 
the ER-2 measurements of the Concorde plume 
[Fahey et al . , 1995a], then the increase in sulfate sur- 
face area is maximized. Under the extreme condition 
of small particles forming from all of the available 
sulfur in the fuel and being emitted into a non vol- 
canically perturbed stratosphere, the calculated ozone 
change for the HSCT perturbation is — 1 %. The 
impact of the production of small particles on 
calculated ozone change depended on both EIno x 
and background chlorine amount. 

If the range of the parametric sensitivity tests is taken 
as an estimate of uncertainty and these results are 
combined, the overall uncertainty due to the pro- 
cesses not included in the basic scenarios is equal to 
or greater than the reaction-rate uncertainty. For all 
of the uncertainties, the model results must be con- 
strained by observations. For example, the extreme 
condition of the tropical pipe must produce tracer 
distributions which are compatible with UARS and 
ER-2 observations. This may reduce the estimates of 
uncertainties, and is an approach which should be 
pursued in the future. 

The Monte Carlo uncertainty propagation showed 
that the maximum uncertainty is in the high-latitude 
lower stratosphere. This is true without any consid- 
eration of PSCs. The uncertainty at high latitudes is 
significantly larger when considering that condensi- 
bles from aircraft exhaust could change the PSC 
formation probability. Of particular concern is the 
possibility that the increase in HNO 3 and H 2 O could 
change the conditions in the Arctic to be more like 
the Antarctic, in which the total NO y levels are sig- 
nificantly reduced during Antarctic spring through 
processes associated with PSC formation that are not 
well understood. Since this loss of NO y is responsi- 
ble for allowing the chlorine-catalyzed ozone de- 
struction to proceed unchecked during the Antarctic 
spring, it is important to improve understanding of 
the possible role of the aircraft effluent in producing 
conditions that would lead to large-scale permanent 
removal of NO y from the Arctic lower stratosphere. 

The evaluations described above were all for a pro- 
jected 2015 atmosphere with low levels of sulfate 
area and two different values of chlorine concentra- 
tion (3 ppbv and 2 ppbv). The 2015 atmosphere it- 
self is not certain. The sulfate surface area may be 
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affected by future (unpredictable) volcanic eruptions. 
Atmospheric CO2 is increasing, which may lead to 
changes in the atmospheric temperature structure and 
to the atmospheric circulation. The effects of these 
conditions may be evaluated with sensitivity studies. 

Not all of the processes are included in this uncer- 
tainty estimate, and including better estimates for 
processes such as the effect of PSCs may lead to a 
higher estimate of overall uncertainty. However, 
there will be a higher level of confidence in the un- 
certainty estimate when all identified processes are 
included realistically. 
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4 SUMMARY 


This report is an assessment of the potential atmo- 
spheric effects of a proposed fleet of HSCT aircraft 
based on the current scientific understanding of the 
atmosphere. Recent advances in that understanding 
allow us to make more definitive statements about 
potential HSCT effects than were possible in past as- 
sessments. The assessment consists of three parts. 
The first is a description of how recent advances 
have improved our understanding of the stratosphere 
and how it works. Second is a set of calculations of 
the potential impacts of an HSCT fleet. These calcu- 
lations were made by 2-D atmospheric models that 
attempt to include all of the latest advances in under- 
standing for the important processes controlling 
ozone. The last is a discussion of the shortcomings 
in our understanding and in the model representa- 
tions of atmospheric processes. These shortcomings 
lead to a discussion of the uncertainties inherent in 
our predictions. These uncertainties can be only 
partially quantified. 

4.1 Plume Measurements 

The evolution of the exhaust plume in the atmo- 
sphere and conversions of the exhaust components 
have been addressed with new in situ aircraft exhaust 
sampling experiments and modeling of engine pro- 
cesses and the dilution of exhaust gases in the atmo- 
sphere. Measurements of NO y species, CN, and CO 2 
were recently made in the exhaust of the Concorde 
supersonic aircraft and the NASA ER-2 subsonic 
aircraft while in flight. The Els were calculated 
independently of plume dynamics because 
measurements of CO 2 provide a measure of fuel 
burned. The measured El of NO x , expressed as 
grams of equivalent NO 2 per kilogram of fuel 
burned, is in excellent agreement with earlier results 
deduced from altitude chamber measurements of 
exhaust from the Concorde engine as part of CIAP. 
These results suggest that the methodology 
developed to calculate the El of NO x at cruise from 
altitude chamber data will be appropriate to evaluate 
new HSCT engines. 

Observations show that NO x is the most abundant 
NO y species in both supersonic and subsonic aircraft 
plumes. The presence of only a small percentage of 
HNO 3 and HONO is in agreement with observations 
of HO x in the plume and theoretical model 
calculations. 


The large number of detectable particles in the 
Concorde plume conflicts with the results of plume 
models that predict the formation of a large number 
of much smaller sulfate particles and a smaller frac- 
tional conversion of SO x to condensed sulfate. 
These results indicate that current understanding of 
plume chemistry and particle formation processes is 
not complete. If new HSCT aircraft produce parti- 
cles at a rate comparable to the Concorde, significant 
increases in particle number and surface area could 
occur in the lower stratosphere of the Northern 
Hemisphere. Including these increases in predictive 
models could change calculated ozone and climate 
impacts from an HSCT fleet. 

4.2 Predicted Changes 

Five 2-D photochemical models were used to 
calculate the impact of an HSCT fleet for a variety of 
cases including sensitivity tests for the El of NO x , 
cruise altitude, background atmospheric chlorine 
amount, and fleet size. Individual models were used 
to test sensitivity to a variety of effects not included 
in these basic predictive calculations. Conditions 
include a baseline fleet of 500 Mach 2.4 aircraft 
burning 8.2 x 10 10 kg of fuel per year in an 
atmosphere with a background stratospheric aerosol 
amount equivalent to a relatively clean period 
between major volcanic eruptions. The models do 
not include the formation or reactivity of PSC 
particles. 

For a fleet of 500 Mach 2.4 HSCTs with EIno x = 5 g 
NC> 2 /kg fuel, the model-calculated effect of the 
emission of NO x and H 2 O on the total column 
amount of ozone in the Northern Hemisphere ranged 
from -0.3% to +0.1%. The calculated total ozone 
change generally increases with increasing emission 
index for NO x . At EInOx = 15 g NC> 2 /kg fuel, the 
range of average values from the five models was 
-1.0% to -0.02%. The calculated total ozone change 
depends on the cruise altitude of the aircraft. For the 
case of 500 HSCTs with EIno x = 5 g N 02 /kg fuel, 
increasing the assumed fleet altitude by 1 km results 
in a model range for the perturbation of -0.5% to 
+0.02% averaged over the Northern Hemisphere. 
Decreasing the cruise altitude by 2 km 
(approximately equivalent to a Mach 2.0 fleet) re- 
sulted in a calculated range of -0.06% to +0. 1 %. 
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For an assumed fleet size doubling to approximately 
1000 aircraft, the calculated perturbation of the aver- 
age total ozone over the Northern Hemisphere for 
EInOx = 5 ranged from -0.7% to +0.03%. At EInOx= 
15 g NC> 2 /kg fuel, the range was from -2.7% to 
- 0 . 6 %. 

Reduction of the assumed background chlorine 
amount from 3 ppbv (projected 2015 concentration) 
to 2 ppbv (projected 2050 concentration) makes the 
calculated HSCT perturbation of Northern 
Hemispheric total ozone somewhat more negative; a 
range of -0.4% to +0.02% for EInOx = 5 g NC> 2 /kg 
fuel. 

The calculated changes in ozone result from a de- 
crease in the middle and upper stratosphere and an 
increase in the lower stratosphere and troposphere. 
The time required for photochemistry and transport 
to effect these ozone changes precludes large ozone 
changes near the assumed HSCT flight corridors. 
The EIjmOx values reflect those expected for future 
HSCT combustors (5 g N02/kg fuel) and those 
closer to current engines (15 g NO 2 /kg fuel). 

Four of the five models used in the calculations were 
tested against two benchmarks, one for photolysis 
rates and one for chemistry. These indicate mostly 
minor differences, implying that the range in model 
results follows, at least in part, from differences in 
the 2-D representation of dynamics. The chemistry 
benchmark tested the model chemistry integrators 
against a photochemical steady-state model that has 
been used in detailed comparisons with atmospheric 
measurements of key radicals involved in the chem- 
istry of stratospheric ozone. 

When the background aerosol surface area used in 
the models is increased by a factor of 4 (roughly 
equivalent to the median amount in the stratosphere 
over the last 2 decades), the calculated perturbation 
to total ozone is less negative by a few tenths of a 
percent. Variations in future stratospheric aerosol 
levels, primarily due to the frequency of large vol- 
canic eruptions, are a source of uncertainty in HSCT 
ozone impact predictions. 

Sulfur emitted by a fleet of HSCTs can increase the 
sulfate aerosol surface area throughout the lower 
stratosphere. Although the increase in surface area 
decreases the impact of NO x emissions, it also in- 
creases ozone loss due to the existing chlorine and 


hydrogen radicals in the stratosphere. If the sulfur 
emitted by HSCTs immediately forms small particles 
in the plume, then the global increase in sulfate sur- 
face area is maximized. 

A model that included aerosol processes was used to 
calculate ozone changes when 100% of the sulfur 
(EIso 2 = 0.4 g sulfur/kg fuel) from an HSCT fleet 
was emitted as small particles into a volcanically un- 
perturbed stratosphere. The change in total column 
ozone from the injection of small particles depended 
in a complex way on EInOx and background chlorine 
amount. For EIno x = 5 g/kg fuel, the calculated 
ozone depletion increased when small particles were 
formed; the overall perturbation became as large or 
larger than that for the EInOx = 15 g/kg fuel case. 
Inclusion of the reaction of B 1 ONO 2 with water on 
aerosols increased the effect of the small particles. 

The model results for ozone change in the basic sce- 
narios vary more widely when heterogeneous reac- 
tions on PSCs are included. This results from uncer- 
tainty in PSC microphysical processes and from 
ambiguity in how to best represent such processes in 
a 2-D model. 

A GCM was used to calculate the global surface air 
temperature and stratospheric temperature responses 
to predicted changes in ozone and water vapor from 
a fleet of HSCTs with EIno x = 15 g NC> 2 /kg fuel. 
Since the resulting changes in temperature are 
smaller than the climatic fluctuations in the model, 
they likely would not be detectable in the 
atmosphere. 

Climate effects due to sulfur and soot emissions, 
while likely to be small, have yet to be evaluated. 
Formation of persistent contrails from HSCTs is ex- 
pected to be infrequent. 

4.3 Model/Measurement Comparisons 

Atmospheric observations can provide tests of our 
understanding of the photochemical and transport 
processes responsible for controlling the distribution 
of ozone and the potential impact of HSCTs. 
Laboratory measurements define the formulation of 
chemistry in models. Subsets of this formulation are 
tested by comparing the model predictions against a 
sufficiently complete set of atmospheric observa- 
tions. Although constraints on transport formula- 
tions in 2-D models are less direct, observations of 
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the distribution of and the correlations between long- 
lived gases provide a test of these formulations. An 
extensive series of 2-D model tests, conducted in 
1992 [Prather and Remsberg, 1993], led to improve- 
ments in the models and an understanding of the dif- 
ferences among them. Comparisons to observations 
made since that time are ongoing, as are further 
model improvements. 

Models reproduce the basic features of the ozone 
column and vertical profile. Consistent with obser- 
vations, the models suggest that ozone concentra- 
tions should have decreased in the past decade. 
However, the predicted changes are smaller than 
those measured by a factor of 1 .5 to 2, particularly in 
the lower stratosphere. Inclusion of a realistic repre- 
sentation of the temporal changes in stratospheric 
sulfate aerosol in a model improves the agreement 
with the observed trends, but the predicted changes 
are still smaller than observed. These differences 
suggest that important uncertainties remain in the 
ozone changes calculated for an HSCT fleet. 

Model simulations reproduce some, but not all, of the 
basic features of the distribution of long-lived trace 
gases in the stratosphere. Since trace gases are the 
sources of the free radicals that catalyze ozone re- 
moval, model-measurement differences reflect inac- 
curacies in the model transport formulation and af- 
fect comparisons of reactive species calculations 
with measurements. 

Important regions for simulation of atmospheric 
transport are those with high spatial gradients in 
long-lived species such as N 2 O, NO y , or ozone. 
Examples are the regions in the stratosphere between 
the tropics and subtropics and at the edge of the polar 
vortices in winter. Aircraft and satellite observations 
have increased the level of understanding of the gra- 
dients and the barriers to transport that create them. 
In general, the 2-D models used to predict HSCT ef- 
fects exhibit weak gradients, suggesting that they 
have too much mixing across these barriers. 

Observations of HO x , NO x , C10 x , and BrO x in the 
stratosphere allow direct evaluation of the photo- 
chemical destruction cycles of ozone. Above 25 km 
altitude, observations of HO x and NO x free radicals 
agree with photochemical models to within -40% 
when the models are constrained by observed abun- 
dances of long-lived chemical species, temperature, 
aerosol surface area, and overhead ozone column. 


Observations of CIO at these altitudes are up to a 
factor of 2 lower than model predictions. Below 20 
km, midday observations of HO x , NO x , C10 x , and 
BiO x also generally agree with similarly constrained 
models to within -40%. However, HO x observations 
exceed model calculations by more than a factor of 3 
near sunrise and sunset, times when HO x concentra- 
tions are small. Incorporation of the reaction of 
B 1 ONO 2 with H 2 O on sulfate aerosols reduces this 
discrepancy but increases the discrepancy in other 
species. Between about 20 and 25 km, observations 
of HO x and C10 x radicals are not available to pro- 
vide a rigorous test of models. 

The observations and modeling of NO x , HO x , and 
C10 x free radicals in the lower stratosphere agree in 
showing that the concentrations of HO x and C10 x 
depend on NO x radicals. When the NO x concentra- 
tion is low, such as observed in the lower strato- 
sphere of the Northern Hemisphere during the fall of 
1992 and spring of 1993, the HO x and C10 x /BrO x 
reactions dominate the local ozone loss rates. Thus, 
according to model calculations, when NO x from 
HSCTs is added to these regions, reductions in HO x 
and C10 x /BrO x abundances generally result in reduc- 
tions in local ozone loss rates. 

Observations show that PSCs occur during the 
Antarctic winter and spring and Arctic winter. 
Heterogeneous reactions on these PSC surfaces acti- 
vate large fractions of the available chlorine that can 
then destroy ozone in the presence of sunlight. 
Widespread ozone destruction does occur in the 
Antarctic because sedimenting PSC particles also 
denitrify the stratosphere and alter the buffering of 
C10 x by NO x . The proper inclusion of PSC pro- 
cesses in the 2-D predictive models would reduce an 
important uncertainty in the HSCT impact results. 

4.4 Uncertainties 

Several classes of uncertainty are associated with 
predictions of HSCT impacts made with 2-D models 
as shown in Table 1. Uncertainties in some model 
parameters can be readily quantified and evaluated 
directly with the models. Other sources of uncer- 
tainty have been identified but not yet quantified. 
Finally, uncertainty from the possibility of missing 
processes or major errors in understanding cannot be 
quantified, but only discussed in terms of the overall 
confidence in the assessment. 
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The only quantified uncertainty in this assessment is 
that due to uncertainties in gas-phase reaction rate 
coefficients. These have been propagated through a 
model to give an estimated uncertainty of ~1% (1-a) 
in the calculated ozone column perturbation in the 
Northern Hemisphere. Uncertainties in heteroge- 
neous chemical reactions and photolysis rates will 
increase the overall chemical uncertainty. 
Constraining these results with atmospheric observa- 
tions will reduce the uncertainty. 

A number of parametric sensitivity tests were per- 
formed for processes not included in the basic as- 
sessment scenarios. If the combined range of these 
tests is taken as an estimate of uncertainty, the over- 
all uncertainty due to the processes not included in 
the basic scenarios is equal to or greater than the re- 
action-rate uncertainty. 

Uncertainties in microphysics include both the for- 
mation of numerous small particles in the HSCT ex- 
haust plume and the possibility of increased denitri- 
fication in polar regions due to increases in conden- 
sible gases from the HSCT exhaust. Although the 
numerous small particles in the HSCT plume can 
indirectly affect ozone loss rates, the particle forma- 
tion process is poorly understood and the uncertain- 
ties are large. Although the likelihood of a signifi- 
cant increase in Arctic denitrification and in polar 
ozone loss through chlorine and bromine chemistry 
is small, the implications of such processes for this 
assessment are large. 

Uncertainties in model transport remain important 
and difficult to quantify. Observations of seasonally 
varying and short- and long-lived tracers provide di- 
agnoses of the degree of mixing and rates of trans- 
port between tropical and mid-latitude regions of the 
stratosphere, as well as providing information on res- 
idence times for stratospheric trace species. 

The impact of ozone and water vapor from HSCTs 
on climate was calculated to be within the climatic 
fluctuations of the model used. Thus, the expected 
perturbation and its uncertainty are difficult to quan- 
tify. 


The possibility of missing processes in the models 
creates an unquantifiable uncertainty in the HSCT 
assessment. However, the probability that important 
processes have been overlooked decreases as we 
continue to explore HSCT perturbations with models 
and continue to challenge these models with atmo- 
spheric measurements. 

4.5 Reducing Uncertainties 

Future efforts to reduce the overall assessment uncer- 
tainty are expected to focus on the major unquanti- 
fied uncertainties and attempts to quantify them. The 
primary methodology for improving our predictive 
capability continues to involve testing model repre- 
sentations of the atmosphere through comparison 
with observational data. Substantial progress in that 
direction will be realized in the near term by system- 
atic comparisons with recent data from UARS, 
ATLAS on the Space Shuttle, and ER-2 aircraft, as 
well as a host of balloon platforms. Specific 
additional strategies for addressing the major 
uncertainties identified in this assessment are listed 
below in no specific order of priority: 

• Extension of the observational database for long- 
lived tracers and radicals to include a number of 
seasonal cycles and undersampled regions of the 
atmosphere, notably the summer high latitudes 
and tropics. 

• Systematic development of 3-D chemical-trans- 
port models for testing transport formulations and 
zonal-symmetry assumptions in 2-D predictive 
models. 

• Inclusion of standardized descriptions of particle 
formation and evolution in 2-D predictive models 
and comparison with global aerosol distributions 
derived from existing satellite and aircraft data. 

• Further characterization, through laboratory and 
atmospheric observations, of the composition and 
temporal evolution of aircraft exhaust aerosols in 
the wake/plume regime, and of the composition 
and microphysics of PSCs. 
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Appendix A 


Changes to the 2-D Predictive Models Since the 1993 NASA 
Interim Assessment of the Atmospheric Effects of 
Stratospheric Aircraft 


aer 

Malcolm Ko, Debra Weisenstein, Jose Rodriguez, and N. D. Sze 
Atmospheric and Environmental Research, Inc. 

1) The vertical resolution has been increased by a factor of 3, from 3.5 km to 1 .2 km. 

2) Temperatures are now the National Meteorological Center (NMC) 8 -year climatology which was 
used as the standard in the 1993 Models and Measurements Workshop [Prather and Remsberg, 
1993]. Tropopause heights are computed monthly based on temperature gradients. 

3) The AER code for calculation of photolysis rates has been replaced by the code of Michael Prather, 
which was used as the standard in the 1995 intercomparison exercise. 

4) The AER model now calculates water (H 2 O) concentrations in the stratosphere. Tropospheric 
values of H 2 O are based on assumed relative humidities. 

5) The reactions rates have been updated to DeMore et al. [1994], The species nitrous acid (HONO) 
and bromine chloride (BrCl) have been added to the model, along with a number of reactions. 
Sticking coefficients from Hanson and Ravishankara [1994] are now used for the reactions 

HC1 + CIONO 2 -» HNO 3 + Cl 2 , 

CIONO 2 + H 2 0 -» HNO 3 + HOC1, and 
HOC1 + HC1 — > H 2 O + Cl 2 on sulfate aerosols. 


CAMED 

J. S. Kinnersley, R. S. Harwood, Y. V. Roberts, and J. A. Pyle 
University of Cambridge/University of Edinburgh 

The model has been improved since Kinnersley and Harwood [1993] in the following ways, as 
documented in Kinnersley [1996]. A spectrum of gravity waves is emitted from the tropopause which 
help to produce realistic stratospheric and mesospheric winds and temperatures; in particular, a 
tropical-winter mesospheric jet and a good semiannual oscillation. The three planetary waves now 
interact with each other, leading to improved behavior [Kinnersley, 1995], and the troposphere has been 
made more interactive. The chemistry has been updated with the JPL-94 rates [DeMore et al., 1994] 
and some more reactions have been included. Heterogeneous chemistry now incorporates Hanson and 
Mauersberger [1988] Type 1 polar stratospheric cloud (PSC) onset and surface area calculations with 
Tabazadeh and Turco [1993] reaction probabilities in the calculation of PSC reaction rates. The 
radiation scheme now includes a treatment of multiple scattering from air molecules. 
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CSIRO 

Keith Ryan, Lakshman Randeniya, Peter Vohralik, and Ian Plumb 
Division of Applied Physics 

Commonwealth Scientific and Industrial Research Organization, Australia 

This model did not participate in the 1993 Atmospheric Effects of Stratospheric Aircraft (AESA) 
Interim Assessment or the 1992 Models and Measurement Workshop. A description of the model is 
provided below. 

Introduction 

The model domain extends from 90°S to 90°N and from the ground to 80 km. Model grid cells are 2 
km by 5 degrees and are centered at 87.5°S, 82.5'S, ..., 87.5°N latitude and 1, 3, ..., 79 km altitude (log- 
pressure vertical coordinates, z* = 16 • log io( 1000/p)), where z* is in km and p is in mbar. Diabatic 
residual circulation together with appropriate diffusion terms provide the dynamics transport. The 
chemistry module is family-based and employs full diurnal averaging methods. Long-lived species are 
transported over 10-day intervals using a transport time step of 0.25 days. After each 10-day transport 
interval, solar fluxes are updated, diurnal averaging is done to calculate diurnal coefficients, and 24- 
hour-average concentrations of equilibrium species are determined (including the partitioning of 
transported family totals). The CSIRO model uses a 360-day year, with updating of solar fluxes, 
diurnal coefficients, and equilibrium species carried out every 10 days. 

Transport 

Derivation of Residual Circulation 

Residual circulation velocities for the model are calculated from heating rates and observed 
temperatures using a full solution of the energy conservation equation as described by Solomon et al. 
[1986b] together with continuity. Temperatures from the ground to 60 km are 8-year zonal- and 
monthly-mean NMC values obtained from the 1992 Models and Measurements database [Prather and 
Remsberg, 1993]. From 60 to 80 km, temperatures are from the Middle Atmosphere Program (MAP) 
database [Barnett and Corney, 1985]. Radiative heating rates from the ground to 60 km are as 
calculated by Rosenfield [1987] and obtained from the 1992 Models and Measurements database 
[Prather and Remsberg, 1993]. Above 60 km, because of the breakdown of local thermodynamic 
equilibrium (LTE) assumptions, cooling rates are calculated using a Newtonian cooling approximation 
with coefficients from Holton and Wehrbein [1980] and heating rates are calculated using 
parameterizations developed by Strobel [1978] and Schoeberl and Strobel [1978]. In the troposphere, 
latent and sensible heating rates are added to the radiative heating rates to obtain net heating rates. 
Sensible and latent heating rates are obtained from Dopplick [1974]. 

Advection 

For a given model run, vertical velocities are calculated for the 12 months of the year at the beginning 
of the run. Both vertical velocities and temperatures are interpolated to the beginning of each 10-day 
transport interval. Horizontal velocities are then computed from the vertical velocities using continuity, 
small mass-conservation correction factors are applied to correct for departures from mass conservation 
(<1%) due to the process-splitting algorithm used to integrate the continuity equation, and horizontal 
eddy diffusion coefficients (K yy ) are derived for the current set of horizontal velocities. 
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Advection is implemented using the variable-order Bott scheme [Bott, 1989]. This method, which is a 
high-order implementation of the upwind scheme, is conservative and positive-definite and produces 
only small numerical diffusion. The 1995 AES A assessment results were obtained using the second- 
order version of the scheme. 

Diffusion 

Stratospheric eddy diffusion coefficients are computed self-consistently with the residual circulation 
using the method of Jackman etal. [1988], namely: 

Kyy = 0.25 • R e - v • tan(latitude) 

where R e denotes the radius of the Earth. Where the calculated values fall outside the range 2 x 10 5 < 
Kyy < 30 x 10 5 m 2 s _1 , they are reset to these limiting values. Stratospheric K yy values between the 
poles and 80°N or S are set to 2 x 10 5 m 2 s _1 . Tropospheric K yy values are set to 10.0 x 10 5 m 2 s _1 . 
Off-diagonal components of the diffusion tensor are included by assuming that diffusion occurs 
predominantly along isentropes [e.g., Jackman etal . , 1988]. 

Vertical eddy diffusion coefficient (K zz ) values are set to 4 m 2 s _1 in the troposphere, 0.25 m 2 s _1 
between 21 and 42 km, and 0.1 m 2 s' 1 throughout the remainder of the stratosphere. 

Tropopause Height 

The CSIRO model uses an annually averaged zonal-mean tropopause height distribution obtained by 
averaging daily NMC tropopause data for the period 1983-1993. Tropopause heights used by the 
CSIRO model for latitudes 87.5°S, 82.5°S, ... , 82.5°N, 87.5°N, using the log-pressure coordinates 
defined above (i.e., z* values in km) are: 


Latitude 

z* 

Latitude 

z* 

Latitude 

z* 

87.5‘S 

9.57 

27.5°S 

14.64 

32.5°N 

12.94 

82.5‘S 

9.54 

22.5°S 

15.24 

37.5°N 

11.68 

77.5'S 

9.50 

17.5°S 

15.48 

42.5°N 

10.70 

72.5‘S 

9.49 

12.5°S 

15.60 

47.5°N 

10.00 

67.5‘S 

9.44 

7.5°S 

15.65 

52.5°N 

9.53 

62.5°S 

9.39 

2.5°S 

15.66 

57.5°N 

9.23 

57.5‘S 

9.47 

2.5°N 

15.64 

62.5°N 

9.02 

52.5'S 

9.78 

7.5°N 

15.63 

67.5’N 

8.82 

47.5’S 

10.32 

12.5°N 

15.59 

72.5°N 

8.63 

42.5°S 

11.06 

17.5°N 

15.47 

77.5°N 

8.47 

37.5°S 

12.07 

22.5°N 

15.14 

82.5°N 

8.37 

32.5°S 

13.43 

27.5°N 

14.29 

87.5°N 

8.33 























































Treatment offyO 

Stratospheric H2O is determined by transporting H2O in the stratosphere with a lower boundary 
condition of 3 parts per million by volume (ppmv) at the tropopause. Tropospheric H2O is set to values 
determined by the Oort climatology where available [Oort, 1983] and the Stratospheric Aerosol and Gas 
Experiment (SAGE II) otherwise [Prather and Remsberg, 1993]. 

Emission Scenarios 

Subsonic and high-speed civil transport (HSCT) emissions for nitrogen oxides (NO x ) (as odd nitrogen, 
NO y ), carbon monoxide (CO), hydrocarbon (HC) (as methane, CH4) and H2O (using EIH2O = 1230 gm 
H20/kg fuel) are calculated using the data files supplied by the National Aeronautics and Space 
Administration (NASA) for the CSIRO grid. 

Additional Tropospheric NOy Source 

An upper tropospheric source of NO y in included in the model. This NO y source is included between 
30°S and 30*N and from 4 to 12 km (model levels at 5, 7, 9, and 1 1 km), and between 25°S and 25°N 
from 12 to 14 km (model level at 13 km). The NO y source strength used is 1100 cnr 3 s* 1 , 
corresponding to a total source of approximately 2 Mt N yr 1 . 

Rainout and Surface Deposition 

Rainout rates are included for transported family totals NO y , inorganic chlorine (Cl y ), and inorganic 
bromine (Br y ), and for equilibrium species H2O2, formaldehyde (CH2O), and methyl peroxide 
(CH3OOH) for model grid cells below 10 km which are also more than 3 km below the tropopause. 
Rainout rates for individual members of NO y , Cl y , and Br y are not explicitly included when partitioning 
these family totals. Rainout rates for model levels at z = 1, 3, 5, 7, and 9 km are 3.74 x 10' 6 , 3.98 x 
10* 6 , 2.96 x 10* 6 , 1.38 x 10* 6 , and 0.42 xlO' 6 s* 1 , respectively. In addition to these rainout processes, 
the following surface deposition rates are included: 6.7 x 10' 6 s* 1 for NO y , hydroxyl radical (OH), 
hydroperoxy radical; (HO2), hydrogen peroxide (H2O2), formyl radical (HCO), methyl radical (CH3); 
and 3.3 x 10* 8 s' 1 for CH2O. (Surface deposition is not applied to odd oxygen (O x ), Cl y , Br y , atomic 
hydrogen (H), and CO as these species are held fixed at the ground by the model.) 

Surface Boundary Conditions 

All surface boundary conditions are independent of latitude and time. At the surface (CSIRO model 
level at z* = 1 km), the boundary conditions are: 

• Molecular hydrogen (H2) = 0.515 ppmv 

• CO = 0.1 ppmv 

• O x = 20 ppbv 

• NO y = flux of 6 x 10 9 cm' 2 s' 1 

• Cl y = 80 pptv 

• Br y = 0.5 pptv 

• CFC-114 = CFC-115 = 

• Other surface mixing ratios as specified by NASA for the Cl y = 2 ppbv and Cl y = 3 ppbv 
atmospheres: 
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Species 

a y = 3 ppbv 

a y = 2 ppbv 

C0 2 

412 ppmv 

412 ppmv 

CH 4 

1994 ppbv 

1994 ppbv 

n 2 o 

330 ppbv 

330 ppbv 

cfci 3 

220 pptv 

120 pptv 

C F 2 C1 2 

470 pptv 

350 pptv 

CFC-1 13 

80 pptv 

60 pptv 

CC1 4 

70 pptv 

35 pptv 

HCFC-22 

250 pptv 

15 pptv 

CH 3 CC1 3 * 

1 1 pptv 

0 

halon-1301 

1.4 pptv 

0.9 pptv 

halon -121 1 

1.1 pptv 

0.2 pptv 

CH 3 C1 

600 pptv 

600 pptv 

CH 3 Br 

15 pptv 

15 pptv 


♦(including HCFC-141b) 


Solar Radiation 

The solar radiation module includes fully interactive, multiple scattering calculations and allows for the 
curvature of the Earth. The method of Meier et al. [1982] is employed to account for the effects of 
multiple scattering, while the spherical atmosphere treatment uses ray-tracing techniques. The model 
uses 158 wavelengths from 175 to 850 nm, as defined in WMO [1986]. Photodissociation cross 
sections were taken from DeMore et al. [1994], except for J(HOBr), which was taken from Orlando and 
Burkholder [1995], 

Chemistry 

Species and Family Definitions 
Transported families are defined as follows: 

O x =0 + 0( 1 D) + 0 3 

NO y = N + NO + N0 2 + N0 3 + HONO + HN0 3 + HNO 4 + 2N 2 0 5 + C10N0 2 + BrON0 2 
Cl y = Cl + 2C1 2 + CIO + C10 2 + OCIO + 2C1 2 0 2 + HOC1 + HC1 + C10N0 2 
Br y = Br + BrO + HOBr + HBr + BrON0 2 

Other species included in the CSIRO model are: 

Fixed Species: N 2 , 0 2 , C0 2 

Transported Species: N 2 0, CH 4 , CO, H 2 0, H 2 , CFC1 3 , CF 2 C1 2 , CFC-1 13, CFC-1 14, CFC-1 15, 
CC1 4 , HCFC-22, CH 3 CC1 3 , CH 3 C1, halon-1301, halon-121 1 , CH 3 Br 

Equilibrium species: H, OH, H0 2 , H 2 0 2 , HCO, CH 2 0, CH 3 , CH 3 0, CH 3 0 2 , CH 3 OOH 




Rate Coefficients 

Cross sections and other reaction-rate parameters are taken from DeMore et al. [1994], except for 
J(HOBr) (as noted above). 

The sulfate aerosol loading is fixed at background levels as defined in Chapter 8 of WMO [1992]. 
Sulfate aerosol reaction rates are calculated using 7^205 + H2O) = 0.1 and y(C10N02 + H2O), as 
defined by Hanson et al. [1994]. For the additional runs which include bromine nitrate (BrONC>2) 
hydrolysis, it was assumed that v(BrON02 + H2O) = 0.4. 

Solution Method 

A typical simulation is carried out using a number of consecutive five-year model runs. During the first 
year, diurnal averaging is carried out every 10 days for all model grid cells which are not in polar night 
to determine diurnal coefficients for all reactions. The same diurnal coefficients are then used for the 
remaining four years of the five-year run when partitioning the transported families and calculating 
other equilibrium species concentrations. Transported species are held constant during the diurnal 
integration, and the 24-hour variation of all equilibrium species is obtained as a by-product. 

Equilibrium Chemistry 

24-hour-average rate coefficients are obtained using the diurnal coefficients and calculated noon photo- 
rate coefficients. These 24-hour-average rate coefficients are used to solve the coupled set of steady - 
state equations dC/dt = 0 for all 24-hour-average equilibrium species concentrations (including 
transported family members). For each transported family, one of the steady-state equations dC/dt = 0 
is replaced by an equation to ensure that the family total concentration equals the sum of its member's 
concentrations. An iterative Newton-Raphson algorithm is used to solve the coupled system of 
equations. 


GSFC 

Charles H. Jackman, David B. Considine, and Eric L. Fleming 
NASA/Goddard Space Flight Center 

1) A look-up table is now used for the photolytic source term (PST) and J(C>2) generated by R. Kawa 
(GSFC) using a radiation code developed by D. Anderson and co-workers at the Johns Hopkins 
Applied Physics Laboratory [Anderson and Meier, 1979; Anderson and Lloyd, 1990]. The PST 
table is a function of wavelength, solar zenith angle, and column ozone and is given on the pressure 
grid of the GSFC two-dimensional (2-D) model. The PST is computed for the particular 
wavelength, solar zenith angle, column ozone, and pressure of interest and is multiplied by the solar 
irradiance at the top of the model to compute the reduced flux at the point of interest. 

2) Sticking coefficients from Hanson and Ravishankara [1994] are now used for the reactions 

HC1 + CIONO2 -> HNO3 + Cl 2 , 

CIONO2 + H 2 0 -> HNO3 + HOC1, and 
HOC1 + HC1 — > H2O + Cl 2 on sulfate aerosols. 
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3) H2O is now transported above 400 mbar rather than fixed over the entire model domain. H2O is set 
to the Oort [1983] climatology below 400 mbar. A detailed description of our treatment of H2O is 
given in Fleming etal. [1995]. 

4) The reactions rates have been updated to DeMore et al. [1994]. 

5) Several minor modifications have been made in the model treatment of the photochemistry of the 
NO y , Cly, and Br y families. 

6) A parameterization and treatment of Type 1 and 2 PSCs was formulated and used in the GSFC 2-D 
model. This PSC treatment is described in detail in Considine et al. [1994]. PSCs were not 
included in the base supersonic transport plane assessments; however, sensitivity studies were run 
for some of these assessments with PSCs included. 


LLNL/UI 

D. E. Kinnison, P.S. Connell, K. E. Grant, D. A. Rotman, and L. Li 
Lawrence Livermore National Laboratory 
K. O. Patten and D. J. Wuebbles 
University of Illinois 

Since the NASA Models and Measurement Workshop [Prather and Remsberg, 1993], the LLNL 2-D 
chemical radiative transport (CRT) model has been improved and refined, with combined support of 
AESA, the NASA Atmospheric Chemistry Modeling and Analysis Program (ACMAP), and the 
Department of Energy. The grid resolution of the model has been increased from 3 km (pressure 
altitude) by ~10 degrees in latitude to 1.5 km and 5 degrees in latitude. The vertical extent of the model 
is extended from 60 to ~85 km to allow for the breaking of planetary and gravity waves, which is 
important to the interactive calculation of the dynamics of the residual circulation. 

The treatment of dynamical processes has been revised to take the 2-D model a major step closer to 
having fully interactive, complete derivation of dynamics and temperature fields. In past studies, an 
interim approach to planetary-wave parameterization using climatological temperature fields, based on 
the Newman method [Newman et al., 1988], was used in the 2-D model to calculate K yy . A significant 
additional step has been the development of a first-principle treatment of planetary waves (both wave 
number 1 and 2), based on the groundbreaking work of Garcia [1991] and Garcia et al. [1992]. This 
improvement, along with gravity-wave effects following Lindzen [1981], have been implemented and 
are currently being used in the LLNL 2-D CRT assessment. These additions have greatly improved the 
dynamical representation of long-lived tracers (e.g., N2O) when compared to satellite data (e.g., the 
Upper Atmosphere Research Satellite). 

We have updated our treatment of reaction on sulfate aerosols as presented in Hanson et al. [1994]. 
Hydrolysis of chlorine nitrate (CIONO2) and dinitrogen pentoxide (N2O5) on the surface of 
stratospheric sulfate aerosol is included. In addition, BrON02 hydrolysis is included as a sensitivity 
study. The LLNL 2-D CRT model can also represent heterogeneous chemical processes to account for 
PSCs. Here, polar heterogeneous chemistry is represented with climatological fields of PSC 
frequencies. This first-order, linearized approach cannot capture the actual situation of processing, 
recovery, and reprocessing that must occur, but is probably adequate to represent the effect of polar 
processing on the global zonally averaged stratosphere. Several sensitivity studies were completed with 
this version of the LLNL 2-D CRT. 
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